PL-TR-97-1143 


PL-TR- 

97-1143 


DYNAMICS  OF  SEMICONDUCTOR  LASERS  UNDER 
EXTERNAL  OPTICAL  INJECTION 

Thomas  B.  Simpson  and  Jia-ming  Liu 


Jaycor,  Inc 

P.O.  Box  85154 

San  Diego,  CA  92186-5154 


August  1997 


Final  Report 


[unc  QUALHY  me?SCTE]D  S 

PHILLIPS  LABORATORY 

Lasers  and  Imaging  Directorate 

AIR  FORCE  MATERIEL  COMMAND 

KIRTLAND  AIR  FORCE  BASE,  NM  87117-5776 


PL-TR-97-1143 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for  any 
purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S. 
Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data,  does  not  license  the  holder  or  any  other  person  or  corporation;  or 
convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that  may 
relate  to  them. 

This  report  has  been  reviewed  by  the  Public  Affairs  Office  and  is  releasable  to  the  National 
Technical  Information  Service  (NTIS).  At  NTIS,  it  will  be  available  to  the  general  public, 
including  foreign  nationals. 

If  you  change  your  address,  wish  to  be  removed  from  this  mailing  list,  or  your  organization  no 
longer  employs  the  addressee,  please  notify  PL/LIMS,  3550  Aberdeen  Ave  SE,  Kirtland  AFB, 
NM  871 17-5776. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notice  on  a 
specific  document  requires  its  return. 

This  report  has  been  approved  for  publication. 


CHRISTOPHER  M.  CLAYTON 
Project  Manager 


FOR  THE  COMMANDER 


J  lUSL^ 

WILLIAlvfT  WELKER,  LtCol  USAF 
Chief,  Optical  Sensing  Division 


ROBERT  A.  DURYEA,  Cbl,  USAF 
Director,  Lasers  and  Imaging  Directorate 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  074-0188 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  this  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information,  including 
suggestions  for  reducing  this  burden  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports.  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington.  VA  22202-4302, 
and  to  the  Office  of  Management  and  Budget,  Papenftrork  Reduction  Project  (0704-01681  Washington.  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  blank) 


2.  REPORT  DATE 

August  1997 


3.  REPORT  TYPE  AND  DATES  COVERED 

Final;  Sep  94  -  Aug  97 _ 


4.  TITLE  AND  SUBTITLE 

Dynamics  of  Semiconductor  Lasers  Under  External  Optical 
Injection 


6.  AUTHOR(S) 

Thomas  B.  Simpson 
Jia-ming  Liu 


5.  FUNDING  NUMBERS 

C:  F29601-94-C-0166 

PE:  62601F 

PR:  3326 

TA:  LH 
WU:  AK 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Jaycor,  Inc 

P.O.  Box85154 

San  Diego,  CA  92186-5154 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

J21 1-97-0071/7211 


9.  SPONSORING  /  MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Phillips  Laboratory/LIMS 
3550  Aberdeen  Ave  SE 
Kirtland  AFB  NM  87117-5776 


10.  SPONSORING  /  MONITORING 
AGENCY  REPORT  NUMBER 

PL-TR-97-1143 


11.  SUPPLEMENTARY  NOTES 

12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited 

12b.  DISTRIBUTION  CODE 

1 3.  ABSTRACT  (Maximum  200  Words) 

Injection  of  an  external  optical  signal  into  the  oscillating  mode  of  a  semiconductor  laser  can  modify  the  coupling  between  the 
circulating  optical  field  and  the  carrier  density  (gain  medium).  The  external  optical  signal  modifies  both  the  gain  and  the  refractive 
index  of  the  semiconductor  gain  medium  by  changing  the  ffee-carrier  density.  Depending  on  the  strength  of  the  optical  injection  and  the 
frequency  offset  between  the  injected  signal  and  the  free-running  oscillation  frequency,  this  can  lead  to  modulated  output  including 
multiwave  mixing,  unstable  dynamics,  and  deterministic  chaos,  or  to  significantly  improved  modulation  and  noise  characteristics 
relative  to  the  free-running  case.  Experimental  investigations  were  coupled  with  numerical  and  analytical  modeling  to  investigate  the 
modified  output  in  a  commercially  available  laser  diode  and  a  vertical-cavity  surface-emitting  laser  (VCSEL).  The  modeling  used  a 
conventional  coupled-equation  model  which,  we  found,  identifies  the  key  phenomena  that  we  observe  in  our  data. 


14.  SUBJECT  TERMS 

semiconductor  lasers,  nonlinear  dynamics,  injection  locking, 
bandwidth  enhancement,  noise  reduction 

15.  NUMBER  OF  PAGES 

134 

16.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 
OF  REPORT 

Unclassified 

18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 

Unclassified 

19.  SECURITY  CLASSIFICATION 

OF  ABSTRACT 

Unclassified 

20.  LIMITATION  OF  ABSTRACT 

Unlimited 

NSN  754(M)1-280-5500 


Standard  Form  298  (Rev.  2-89) 

Prescftbed  by  ANSI  Std.  Z39-1B 


CONTENTS 


Section  Page 

HGURES  iv 

ACKNOWLEDGEMENTS  vii 

1.0  INTRODUCTION  1 

2.0  COUPLED  EQUATION  MODEL  6 

2.1  SMALL-SIGNAL  MODULATION  ANALYSIS  7 

2.1.1  Theoretical  Model  8 

2.1.2  Linearized  Analysis  16 

2.1.3  Modulation  Spectra  26 

2.2  LINEARIZED  ANALYSIS  INCLUDING  FEEDBACK  32 

2.3  NUMERICAL  SIMULATION  OF  INJECTION- 

LOCKED  OPERATION  38 

2.3.1  Resonance  Frequency  Enhancement  and  Noise  Reduction  43 

2.3.2  Dynamic  Range  47 

3.0  EXPERIMENTAL  APPARATUS  51 

4.0  SEMICONDUCTOR  LASER  ABOVE  THRESHOLD  60 

4.1  NONLINEAR  DYNAMICS  63 

4.2  ENHANCED  MODULATION  BANDWIDTH  AND  REDUCED 

NOISE  I  INJECTION-LOCKED  SEMICONDUCTOR  LASERS  75 

4.3  FEEDBACK  EFFECTS  ON  MODULATION  SPECTRA  IN 

INJECTION  LOCKED  SEMICONDUCTOR  LASERS  8 1 

4.4  VCSEL  ABOVE  THRESHOLD  86 

4.5  VCSEL  POLARIZATION  AND  NOISE  CHARACTERISTICS  92 

5.0  VCSEL  NEAR  THRESHOLD  101 

5.1  FREE-RUNNING  OPERATION  101 

5.2  OPERATION  UNDER  EXTERNAL  OPTICAL  INJECTION  1 08 

6.0  CONCLUSIONS  119 


iii 


FIGURES 


Figure  Page 

1  Transceiver  concept  using  injection-locked  semiconductor  lasers  3 

2  Mapping  of  the  locking  and  stability  characteristics  of  a  semiconductor  laser  14 

3  Frequency  offsets  and  damping  rates  as  a  function  of  the  injection  parameter 

along  lines  of  constant  phase  1 5 

4  Frequency  offsets  and  damping  rates  as  a  function  of  the  locked  phase  for 

constant  values  of  the  injection  parameter  21 

5  Spectrum  of  the  regenerative  amplification  sideband  characteristic  due  to  a 

weak  optical  modulation  29 

6  Spectrum  of  the  four-wave  mixing  sideband  characteristic  due  to  a  weak 

optical  modulation  30 

7  Spectrum  of  the  amplitude  modulation  characteristic  due  to  a  weak 

optical  modulation  3 1 

8  Spectrum  of  the  optical  field  positive  sideband  characteristic  due  to  a 

weak  current  modulation  33 

9  Spectrum  of  the  amplitude  modulation  characteristic  due  to  a  weak 

current  modulation  34 

10  Mapping  of  different  characteristic  operation  regions  of  the  laser  under 

injection  42 

1 1  Resonance  frequency  of  the  laser  under  injection  locking  as  a  function  of 

the  injection  parameter  44 

12  Total  noise  power  of  the  laser  under  injection  locking  46 

1 3  Saturation  and  nonlinear  frequency  generation  due  to  a  large  modulation 

current  49 

14  Schematic  diagram  of  the  experimental  apparatus  52 

1 5  Spectra  of  the  regenerative  amplification  sideband  due  to  a  weak  optical 

probe  showing  the  effects  of  a  weak  optical  cavity  for  the  probe  56 

1 6  Output  power  of  the  VCSEL  under  study  at  bias  currents  where  the  output 

power  is  in  one  spatial  mode  with  two  orthogonal  polarization  modes  58 

1 7  Mapping  of  the  locking  and  stability  characteristics  of  the  conventional,  edge 

emitting  semiconductor  laser  61 

1 8  Calculated  eigenfrequency  and  measured  frequency  of  the  optical  spectrum 

sideband  peak  for  the  conventional,  edge-emitting  semiconductor  laser  62 


IV 


FIGURES  (Continued) 


24 

25 

26 


27 


28 

29 

30 

31 

32 

33 

34 


Pa2e 


Measured  optical  spectra  of  the  conventional,  edge-emitting  semiconductor 
laser  under  optical  injection  at  the  free-running  optical  frequency  64 

Numerically  calculated  bifurcation  diagram  of  the  extreme  of  the  normalized 
optical  field  amplitude  for  injection  is  at  the  free-running  frequency  67 

Computed  optical  spectra  showing  the  effects  of  the  spontaneous-emission 
noise  at  three  levels  of  injection  69 

Measured  optical  spectra  of  the  conventional,  edge-emitting  semiconductor 
laser  under  optical  injection  at  a  frequency  offset  ~  5.5  GHz  70 

Mapping  of  the  experimentally  observed  dynamic  regions  for  a  conventional 
edge-emitting  semiconductor  laser  72 

Frequency  pushing  of  the  initially  free-running  spectral  feature  as  a  function 
of  injection  with  the  master  laser  frequency  offset  by  5.5  GHz  74 

Spectra  of  the  regenerative  amplification  sideband  due  to  a  weak  optical 

probe  under  four  different  operating  conditions  78 

Calculated  amplitude  modulation  spectra  due  to  a  weak  current  modulation 
corresponding  to  the  regenerative  amplification  spectra  plotted  in  Figure  25  80 

Changes  in  the  amplitude  noise  spectrum  of  a  quantum-well  laser  diode,  where 
the  dominant  noise  source  is  spontaneous  emission,  due  to  stable  injection 
locking  by  an  external  laser  82 

Spectra  of  the  regenerative  amplification  sideband  due  to  a  weak  optical 

probe  showing  the  effects  of  feedback  84 

Spectra  of  the  amplitude  modulation  due  to  a  weak  optical  probe  showing  the 
effects  of  feedback  85 

Spectra  of  the  regenerative  amplification  sideband  characteristic  due  to  a  weak 
optical  modulation,  showing  the  effects  of  optical  feedback  87 

Spectra  of  the  amplitude  modulation  characteristic  due  to  a  weak  current 
modulation  88 

Measured  optical  spectra  of  the  VCSEL  under  optical  injection  at  the  free 
running  frequency  89 

Variation  of  the  resonance  frequency  of  the  VCSEL  under  external  optical 
injection  91 

Measured  optical  spectra  of  the  VCSEL  with  as  the  injection  frequency  offset 
is  shifted  93 


V 


! 


FIGURES  (Continued) 

Figure  Page 


35  Low-frequency  power  spectrum  showing  the  resonance  due  to  the  antimode 

coupling  between  the  two  orthogonally  polarized  modes  95 

36  Optical  spectra  showing  the  effect  of  external  optical  injection  on  the  two  • 

polarization  components  97 

37  Optical  spectra  showing  the  effects  of  optical  injection  when  the  two 

polarization  components  are  of  similar  magnitude  98 

38  Power  spectrum  at  the  beat  frequency  between  the  injection-locked 
optical  signal  component  and  the  shifted  Fabry-Perot  cavity  resonance 
feature  for  the  laser  under  optical  injection  where  the  optical  spectrum 

showed  two  strong  optical  components  99 

39  Measured  optical  spectra  of  a  VCSEL  as  the  bias  current  is  varied  102 

40  Calculated  optical  spectra  and  histograms  of  the  Fabry-Perot  resonance 

frequency  variation  based  on  fluctuations  in  the  carrier  density  105 

41  Calculated  optical  spectra  and  histograms  of  the  Fabry-Perot  resonance 
frequency  variation  based  on  fluctuations  in  the  carrier  density  using 

modified  values  of  the  experimentally  determined  VCSEL  parameters  106 

42  The  optical  spectrum  of  the  free-running  VCSEL  biased  well  above 
threshold  and  the  regeneratively  amplified  spectrum  due  to  the  injection 

of  a  weak  optical  probe  1 10 

43  Representative  optical  spectra  of  the  free-running  and  optically-injected 

VCSEL  near  threshold  1 1 1 

44  The  dependence  of  the  injection-locked  power  and  the  shift  of  the  resonance 

from  the  free-running  frequency  for  the  VCSEL  biased  near  threshold  1 1 3 

45  Measured  optical  spectra  of  the  VCSEL  biased  at  4. 1  mA  and  subjected  to 

external  optical  injection  1 15 

46  Calculated  optical  spectra  of  the  VCSEL  biased  at  4. 1  mA  and  subjected  to 

external  optical  injection  116 

47  Spectra  of  the  VCSEL  under  strong  optical  injection  as  the  detuning  is 

varied  117 

48  Measured  pushing  of  the  Fabry-Perot  resonance  as  a  function  of  the  * 

detuning  of  the  injected  signal  118 


VI 


ACKNOWLEDGEMENTS 


This  report  represents  the  conclusion  of  the  formal  collaborative  effort  between  the  au¬ 
thors  and  the  scientists  of  the  Nonlinear  Optics  Center  of  Phillips  Laboratory.  While  the 
report  concentrates  on  the  experimental  meeisurements,  analytical  modelling  and  numerical 
calculations  performed  by  the  authors,  many  of  these  results  build  upon  the  collaborative 
papers  written  with  our  colleagues  at  Phillips  Laboratory  which  are  cited  at  the  appropriate 
points  in  the  text.  We  especially  want  to  express  our  appreciation  to  Dr.  Christopher  M. 
Clayton  and  Dr.  Athanasios  Gavrielides  who  not  only  served  as  technical  monitors  but  as  key 
collaborators  with  us  in  this  research.  They  provided  key  technical  guidance  which  helped 
to  shape  the  direction  of  this  research  while,  simultaneously,  they  gave  us  the  freedom  to  in¬ 
vestigate  unexpected,  exciting  new  findings  that  appeared  along  the  way.  Dr.  Gavrielides  is 
also  coauthor  of  the  early  paper  describing  the  bandwidth  enhancement  phenomena.  Based 
on  the  findings  in  this  paper,  a  patent  is  being  filed  on  behalf  of  the  Air  Force.  We  consider 
ourselves  to  have  been  remarkably  fortunate  to  have  had  the  opportunity  to  pursue  this  re¬ 
search  effort  under  their  oversight  and  with  their  participation.  We  would  also  like  to  thank 
Dr.  Vassilios  Kovanis  for  many  discussions  and  his  active  collaboration  in  the  comparison 
of  experimental  data  and  numerical  calculations.  He  prepared  the  bifurcation  diagram  in 
Figure  20,  performed  the  numerical  calculations  shown  in  Figure  21,  and,  more  generally, 
performed  many  numerical  calculations  which  showed  the  ability  of  the  model  to  recover 
experimentally  generated  spectra.  Also,  his  enthusiasm  was  infectious  and  enlivened  many 
discussions.  Finally,  some  of  the  work  described  in  this  report  was  also  supported,  in  part, 
by  Phillips  Laboratory  under  contract  F29601-96-C-0047  and  by  the  Army  Research  Office 
under  contract  DAAH04-96-C-0038. 


Vll 


viii 


1.0  INTRODUCTION 


Semiconductor  lasers  have  been  the  subject  of  considerable  research  and  development  activity 
over  the  past  30  years  and  now  are  used  in  a  variety  of  applications  including  compact 
disc  players  and  telecommunications.  They  are  also  important  as  a  test  system  for  the 
understanding  of  many  key  features  of  laser  operation.  Due  to  the  small  size  of  the  optical 
cavities  that  can  be  constructed  with  semiconductor  material,  and  the  high  gain  of  that 
material,  features  of  the  optical  spectrum  can  be  related  back  to  fundamental  quantum 
noise  sources  of  the  system  [1].  Semiconductor  lasers  differ  from  many  other  laser  systems 
in  that  the  free  carriers  which  control  the  optical  gain  of  the  system  also  have  a  strong 
effect  on  the  refractive  index,  giving  it  the  characteristics  of  a  detuned  oscillator  [2].  This 
feature  leads  to  an  enhancement  of  the  linewidth  of  a  semiconductor  laser  [3]  and  also  has 
been  exploited  to  induce  unstable  dynamics  and  deterministic  chaos  in  semiconductor  lasers 
by  the  addition  of  an  extra  external  modulation  [4]  such  as  a  current  modulation  [5],  weak 
optical  feedback  [6]  or  external  optical  injection  [7]. 

As  a  test  system  for  studying  nonlinear  dynamics  in  semiconductor  lasers,  the  laser  undergo¬ 
ing  external  optical  injection  has  certain  advantages.  The  use  of  an  independent  leiser,  rather 
than  feedback  from  the  laser  under  investigation,  adds  the  freedom  to  control  the  amplitude 
and  offset  frequency  of  the  perturbing  optical  field.  Optical  injection,  as  opposed  to  current 
modulation,  does  not  require  high-speed  electronics  and  bypasses  the  electronic  parasitic  ef¬ 
fects  that  can  obscure  the  intrinsic  laser  properties.  Experimental  investigations  of  nonlinear 
dynamics  in  a  conventional  edge-emitting  Fabry-Perot  laser  diode  biased  well  above  thresh¬ 
old  have  shown  quantitative  agreement  with  a  coupled  equation  model  for  single-mode  laser 
operation  [7,  8,  9,  10].  Here,  we  summarize  that  work  and  present  the  results  in  the  form 
of  a  mapping  of  the  observed  dynamics  of  the  laser  output  as  a  function  of  the  amplitude 
and  detuning  of  the  external  optical  injection.  Many  features  of  the  observed  dynamics  are 
expected  to  appear  in  a  variety  of  laser  structures,  while  some  relate  to  characteristics  of 
the  specific  laser  structure.  We  also  present  results  on  optical  injection  of  a  vertical  cavity 
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surface  emitting  laser  (VCSEL).  This  laser  displayed,  in  its  free-running  condition  when  bi¬ 
ased  near  threshold,  an  optical  output  with  a  near-Gaussian,  as  opposed  to  a  Lorentzian, 
spectrum.  When  the  external  optical  injection  destabilizes  the  dynamics  of  this  system,  the 
optical  spectrum  shows  spectral  spikes  and  holes  in  the  Gaussian  envelope.  As  before,  all  of 
the  results  can  be  reproduced  with  the  coupled  equation  model. 

In  recent  work,  we  have  proposed  that  the  modulation  bandwidth  and  broadband  noise 
characteristics  of  a  conventional,  edge-emitting  laser  diode  are  significantly  improved  under 
appropriate  injection  locking  to  a  stable  master  laser  [11 , 12].  By  strong,  we  mean  an  injection 
power  on  the  order  of  10"^  to  2-3  x  10"^  of  the  output  power  of  the  free-running  slave  laser. 
Calculations  have  indicated  that  the  modulation  bandwidth  of  the  injection-locked  laser 
could  be  improved  beyond  the  theoretical  limit  of  the  free-running  laser.  This  bandwidth 
improvement  can  be  accompanied  by  significant  reduction  in  the  noise  power  of  the  laser. 
Significantly,  the  injected  optical  power  required  to  induce  these  changes  can  be  only  a  few 
percent  of  the  free-running  output  power  of  the  laser.  Injection  locking  of  a  semiconductor 
laser  also  allows  one  to  tailor  the  response  characteristics  to  emphasize  a  specific  feature 
such  as  modulation  bandwidth,  flatness  of  response,  or  gain  at  a  particular  modulation 
frequency.  This  capability  can  be  used  to  generate  a  wide  variety  of  output  characteristics 
with  controlled  amplitude  and  phase.  A  block  diagram  of  a  generic  transceiver  element 
built  around  an  injection-locked  laser  diode  is  shown  in  Figure  1.  The  master  laser  is  a 
stable,  narrowband  laser  oscillator  which  may  have  very  poor  high-frequency  modulation 
characteristics.  Between  the  master  and  slave  laser,  or  multiple  slave  lasers,  are  optical 
elements  which  provide  the  proper  coupling  between  the  lasers  and  external  signals  and 
minimize  unwanted  feedback.  The  transceiver  can  be  operated  in  the  transmit  mode  where 
an  electronic  signal  can  be  added  to  the  bias  current  of  the  slave  laser  to  produce  an  analog 
or  digital  modulation  signal.  Because  the  slave  laser  is  operated  as  an  injection-locked 
oscillator,  it  can  have  relatively  poor  free- running  noise  characteristics  to  emphasize  high- 
frequency  response.  In  the  receive  mode,  the  slave  laser  becomes  a  bandwidth  tunable  filter 
and  preamplifier  for  modulated  optical  input  to  the  photodetector.  The  master  laser  controls 
the  central  wavelength  of  the  filter.  By  adjusting  the  bias  of  the  slave  laser,  the  bandwidth 
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Figure  1:  Transceiver  concept  using  injection-locked  semiconductor  lasers. 
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of  the  filter  can  be  controlled.  Because  the  master  laser  only  has  to  inject  a  locking  signal 
of  a  few  percent  of  the  power  of  the  slave  laser,  a  single  master  laser  can  control  multiple 
transceiver  elements. 

In  this  report,  we  summarize  our  findings  on  the  nonlinear  dynamics  and  bandwidth  en¬ 
hancement  which  can  be  induced  in  a  semiconductor  laser  under  external  optical  injection. 
In  the  next  section,  we  review  the  model  of  a  semiconductor  laser  subject  to  a  strong  optical 
signal,  first  describing  the  full  nonlinear  coupled  equations  and  the  proceeding  to  a  pertur¬ 
bation  analysis  of  the  modulation  characteristics  of  an  injection-locked  semiconductor  laser. 
This  is  followed  by  descriptions  of  our  experimental  apparatus  and  some  general  features  of 
the  specific  semiconductor  lasers  under  investigation.  Because  nonlinear  dynamics  plays  a 
critical  role  in  the  understanding  of  the  semiconductor  laser  under  strong  optical  injection, 
we  then  discuss  a  detailed  mapping  of  the  kinds  of  output  dynamics  that  we  observe  under 
external  optical  injection.  The  mapping  is  generated  by  comparing  optical  spectra  generated 
from  calculated  time-series  data  with  and  without  noise  terms  with  the  spectra  we  observe. 
Removing  the  noise  terms  allows  us  to  identify  the  underlying  deterministic  dynamics  in  a 
way  which  cannot  be  duplicated  experimentally.  We  then  discuss  some  interesting  features  of 
the  VCSEL  which  were  not  observed  in  our  previous  studies  of  a  conventional,  edge-emitting 
laser  diode,  with  emphasis  on  VCSEL  characteristics  near  threshold.  There  has  been  consid¬ 
erable  controversy  about  the  description  of  microcavity  lasers,  like  VCSELs,  near  threshold. 
However,  we  show  how  our  model  recovers  all  of  the  novel  features  that  we  observe.  Next, 
we  discuss  the  output  characteristics  of  an  injection-locked  semiconductor  laser  under  strong 
injection  where  stable  dynamics  are  expected  to  be  observed.  We  confirm  the  capability 
of  the  model  to  predict  the  enhancement  of  the  characteristic  resonance  frequencies  of  the 
injection-locked  laser  and  the  noise  reduction.  The  model  presented  in  past  publications 
[7,  12,  13]  is  also  augmented  by  the  inclusion  of  feedback  from  optical  elements.  As  we 
describe  below,  we  found  that  feedback  effects  could  be  observed  in  our  data.  However, 
the  coupled-equation  model  correctly  shows  the  influence  of  feedback.  We  then  describe 
the  measurements  which  showed  that  injection  locking  could  be  used  to  achieve  modulation 
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bandwidths  beyond  the  theoretical  limit  of  the  free-running  laser  diode 
back,  the  enhanced  operating  characteristics  could  still  be  observed, 
with  a  brief  discussion  of  possible  future  research. 


.  Even  with  the  feed- 
Finally,  we  conclude 
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2.0  COUPLED  EQUATION  MODEL 


Semiconductor  lasers  have  been  the  subject  of  extensive  modelling.  Like  a  variety  of  other 
laser  systems,  the  polarization  of  the  semiconductor  gain  medium  has  a  relaxation  rate  which 
dominates  all  other  dynamic  frequencies  in  the  system,  except  the  optical  emission  frequency. 
In  these  lasers,  called  Class  B  lasers,  the  differential  equation  describing  the  polarization  of 

V 

the  medium,  P,  can  be  adiabatically  eliminated  by  writing  P  explicitly  as  a  function  of  other 
parameters  of  the  laser: 

P  =  -eo—{h  +  i)gAe-'^^^  (1) 

Here,  the  gain,  and  the  linewidth  enhancement  factor,  b,  relate  the  complex  intracavity 
field,  A,  and  the  polarization  induced  by  the  free  carriers  which  make  up  the  gain  medium.  In 
Equation  1,  Co  is  the  permittivity  of  free  space,  n  is  the  refractive  index  of  the  semiconductor 
medium,  c  is  the  speed  of  light  and  ujq  is  the  optical  frequency  of  the  oscillating  mode.  The 
gain  is  assumed  to  obey  a  linear  dependence  on  changes  in  the  carrier  density,  with  coefficient 
and  changes  in  the  photon  density,  with  coefficient  g.p,  about  an  appropriately  chosen 
operating  point  [13]. 


An  isolated,  single-mode  laser  of  this  type  that  is  biased  well  above  threshold  is  sufficiently 
described  by  only  two  rate  equations:  one  for  the  photon  density  and  another  for  the  pop¬ 
ulation  density  of  the  gain  medium  [2,  14).  When  the  laser  is  coupled  to  the  outside  world 
and  is  subject  to  the  perturbation  of  external  sources,  such  as  the  vacuum  field  or  external 
optical  signals,  an  equation  for  the  complex,  oscillating  electric  field  must  replace  the  simpler 
equation  for  the  photon  density.  The  two  equations,  one  for  the  complex  oscillating  field 
within  the  laser  cavity  and  one  for  the  carrier  density,  which  describe  the  dynamics  of  the 
single-mode  laser  under  external  optical  injection  are  [7,  13,  15]: 


dt 


+  i{uQ  -  WcjA  -I-  -(1  -  ih)gA  +  F  +  rjAi  exp{-iQt), 


dt 


2con^ 

Ku>q 


9\A?- 


(2) 

(3) 
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Here,  7c  is  the  cavity  decay  rate,  is  the  longitudinal  mode  frequency  of  the  cold  laser 
cavity,  and  F  is  the  confinement  factor  which  gives  the  spatial  overlap  between  the  active 
gain  volume  and  the  optical  mode  volume.  TV  is  the  carrier  density,  J  is  the  injection  current 
density,  e  is  the  electronic  charge,  d  is  the  active  layer  thickness,  7^  is  the  spontaneous  carrier 
decay  rate,  and  n  is  the  refractive  index  of  the  semiconductor  medium.  Two  source  terms 
in  the  field  equation  are  due  to:  1)  external  optical  injection  with  rj  being  the  coupling 
parameter,  Ai  the  amplitude  of  the  external  field  and  0  the  offset  frequency  with  respect  to 
the  free-running  frequency  of  the  injected  laser,  and  2)  noise.  F  is  the  complex  Langevin 
source  term  for  the  field  noise  due  to  spontaneous  emission  into  the  mode  and  external 
coupling  out  of  the  mode  [3,  16]. 


{Fit))  =  0 

(4) 

{F(t)F‘it'))  =  R.,Sit  -  !') 

(5) 

{F{t)F(t'))  =  0 

(6) 

where  R,p  gives  the  strength  of  the  field  source  term.  The  noise  source  terms  in  the  carrier 
equation  are  omitted.  They  are  relatively  unimportant  in  determining  the  optical  field 
spectra  which  we  will  emphasize.  These  equations  form  the  basis  for  a  full,  nonlinear  analysis 
of  the  semiconductor  laser  under  external  optical  injection.  The  full  nonlinear  form  must 
be  retained  when  the  laser  output  departs  from  stable  operation  under  external  optical 
injection.  The  equations  have  proven  to  be  quite  robust  and  can  be  applied  for  arbitrary 
values  of  the  injected  field  both  inside  and  outside  the  range  of  injection-locked  operation. 
While  variations  in  spatial  characteristics  can  be  retained  in  this  model  [17],  generally  it  is 
assumed  that  spatial  inhomogeneities  can  be  averaged  over  the  cavity  and  a  lumped  element 
approach  can  be  used. 

2.1  SMALL-SIGNAL  MODULATION  ANALYSIS 

In  previous  publications,  we  have  described  the  weak  injection  limit  and  pointed  out  how  key 
dynamic  parameters  of  a  semiconductor  laser  can  be  determined  from  the  optical  modulation 
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spectra  [13,  25].  We  have  also  described  similarities  between  the  noise  spectra  and  the 
weak  external-field  modulation  spectra  of  a  free-running  semiconductor  laser  [15].  In  these 
analyses,  we  used  the  lumped-element  model  of  the  semiconductor  laser,  emphasizing  certain 
key  features.  The  modulation  characteristics  are  determined  from  a  perturbation  analysis 
about  the  injection-modified  operating  point.  We  next  make  a  more  complete  explication  of 
the  small  signal  analysis  for  a  semiconductor  laser  subject  to  a  strong,  externally  injected 
field  which  satisfies  the  conditions  for  injection  locking.  This  small-signal  analysis  uses  both 
a  weak  injection  current  and  a  second,  weak  optical  field  as  the  modulation  sources.  Both 
optical  and  current  modulation  characteristics  are  of  interest  in  a  variety  of  applications.  We 
have  attempted  to  use  notation  consistent  with  our  earlier  work  on  weak  optical  injection  in 
semiconductor  lasers  [13].  The  goals  will  be  to  further  explain  the  similarities  and  differences 
between  the  modulation  characteristics  of  the  free-running  and  the  injection-locked  lasers 
and  to  show  the  effects  of  a  strong  external  field  on  the  modulation  characteristics  of  the 
laser.  Specifically,  we  emphasize  the  strong  optical  injection  regime  where  enhancement  of 
the  modulation  bandwidth  and  a  broad  tuning  range  for  stable  injection  locking  can  be 
achieved.  Past  work  on  the  modulation  characteristics  of  injection-locked  semiconductor 
lasers  has  not  emphasized  this  regime  [18,  19,  20,  21,  22].  We  will  also  compare  the  results 
of  our  analysis  of  modulation  in  an  optical  cavity  with  past  analyses  of  four-wave  mixing 
in  semiconductor  travelling- wave  amplifier  configurations  [23,  24]  to  show  that  the  optical 
cavity  plays  a  key  role  in  the  modified  modulation  characteristics.  We  will  be  able  to  use 
some  of  the  results  described  here  to  analyze  the  noise  spectra  of  a  semiconductor  laser 
subject  to  strong  optical  injection.  A  more  complete  analysis  of  the  noise  characteristics 
requires  a  full  nonlinear  treatment  and  will  be  described  later  in  this  report. 

2.1.1  Theoretical  Model 

In  our  modulation  analysis,  we  will  consider  a  single-mode  semiconductor  laser  subject  to 
external  optical  injection  and/or  current  modulation  [7,  13,  25].  When  free  running,  the 
laser  is  assumed  to  be  biased  above  the  threshold  for  oscillation.  We  replace  the  noise  source 
term  in  the  field  equation  with  a  second  external  optical  injection  field.  This  field  is  assumed 
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to  be  sufRciently  weak  so  that  it  will  act  as  a  linear  probe.  A  modulation  current  will  also 
be  added  to  the  injection  current.  The  modulation  of  the  injection  current  will  also  be  small 
and  be  used  as  a  linear  probe.  The  noise  terms  are  dropped  for  the  modulation  analysis. 
The  interaction  can  be  described  by  two  coupled  equations: 

^  =  —'^A  +  i{u>L  —  uJc)A  +  —{l—ib)gA  +  T]{Ai-‘rAie  (7) 

at  2  2 

^  ^  J  +  J^cos(m)  (8) 

dt  ed  *  huL 

where  A  is  the  total  complex  intracavity  field  amplitude  at  the  locked  oscillating  frequency 
ul,  7c  is  the  cavity  decay  rate,  iVc  is  the  longitudinal-mode  frequency  of  the  cold  laser  cavity, 
r  is  the  confinement  factor,  A;  and  A,-  are  the  amplitudes  of  the  strong  injection  signal  at 
the  locking  frequency  and  the  weak  probe  signal  at  ul  +  respectively,  t]  is  the  coupling 
rate,  N  is  the  carrier  density,  J  is  the  steady-state  injection  current  density,  Jm  is  the  weak 
current  modulation  at  a  frequency  f],  e  is  the  electronic  charge,  d  is  the  active  layer  thickness 
and  is  the  spontaneous  carrier  decay  rate. 

The  exact  functional  dependence  of  the  gain  on  parameters  such  as  the  carrier  density  and 
the  optical-field  intensity  is  unknown,  but  it  is  important  to  consider  the  gain  coefficient 
g{NL,SL)  as  a  function  of  the  carrier  density,  A^l,  and  the  intracavity  photon  density.  Sl 
[13].  The  subscript  L  indicates  the  values  of  the  parameters  in  the  injection-locked  condition. 
Therefore,  we  make  the  assumption  that  the  gain  can  be  modeled  perturbatively  about  a 
steady-state  operating  point  determined  by  the  injection  current,  J,  the  carrier  density,  .V^. 
the  field  amplitude,  |Al|,  and  the  phase  relative  to  the  strong  injection  field,  <I>l.  The  photon 
density  is  related  to  the  field  amplitude  through: 

5i  ^  (9) 

nuJL 

The  steady-state  operating  point  is  found  by  setting  the  time  derivatives  and  the  weak- 
modulation  terms  to  zero  in  Equations  7  and  8.  At  the  operating  point,  this  procedure 
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yields: 


I^lI 


2‘q\ Ai\  cos  <t>L 
7c  -  r^L 


sm 


-1 


’{uL  -u)c-  hTgLl2)\Ai\ 

v\H 


Nl 


-i- 

Ka  \cd 


(10) 

(11) 

(12) 


The  dependence  of  the  gain  on  the  carrier  density  and  the  photon  density  is  included  in 
the  analysis  through  two  new  relaxation  rates,  'jnL  and  'fpL.  While  7^  and  7c  are  assumed 
to  be  constant,  7„l  and  7pL  are  proportional  to  the  intracavity  photon  density  through  the 
following  relations: 


InL  =  OnSh  = 


dg 

dN 


Sl  and  7pL 


Nl 


=  -r  I 


Sl, 


(13) 


Nl 


where  and  gp  are  the  differential  and  nonlinear  gain  parameters  defined  as  the  derivatives  of 
the  gain  coefficient,  g,  with  respect  to  the  carrier  density  and  the  photon  density,  respectively, 
evaluated  at  the  operating  point. 


These  definitions  can  be  related  back  to  our  previous  definitions  used  for  weak  injection  in  a 
free-running,  single-mode  laser  by  setting  Ai  =  0,  returning  the  operating  point  to  its  free- 
running  value  [13].  This  requires  ul  wo,  the  free-running  optical  frequency,  Al  ^  Aq,  the 
free-running  average  field  amplitude,  and  Nl  =>  N,  the  free-running  average  carrier  density, 
as  previously  defined.  Because  the  dynamic  parameters  of  free-running  laser  diodes  can  be 
experimentally  determined  [13,  26],  we  can  use  these  dynamic  parameters  in  an  analysis  of 
the  injection-locked  laser.  First,  we  must  relate  the  oscillation  points  of  the  free-running  and 
the  injection-locked  lasers.  Within  the  injection-locked  region,  the  frequency  and  phase  of 
the  slave  laser  are  shifted  relative  to  the  free-running  values.  The  frequency  and  phase  shifts 
are  related  by  [18,  19]: 


u>L  —  4*  V. 


(14) 
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Here,  U  =  t]\Ai/Al\  cos  <I>l  is  the  in-phase  injection  term  and  V  =  T]\Ai/AL\s\n(l>L  is  the 
quadrature  term.  From  Equations  10  and  11,  the  gain  defect  at  the  injection-locked  operating 
point  is  2U  and  the  frequency  offset,  uil  —  u>c  —  bTgi/2,  is  V.  Equation  14  places  an  upper 
bound  on  the  detuning  range  for  injection-locked  operation:  \u)i  —  wqI  <  vl  +  b'^rj\Ai/ Al\. 
The  injection-locked  range  is  further  limited  by  the  requirement  that  the  injection-locked 
solution  has  a  lower  gain  than  the  free-running  solution  [19].  The  injection-locked  solution 
is  always  in  competition  with  an  unlocked  solution  with  unbounded  phase.  From  Equation 
10  this  requires  that  |<^£,|  <  5r/2,  or  — \/l  +  b'^rjlAi/Ail  <  u>l  —  u)o  <  rjlAi/Ail- 

The  relation  between  \Al\  and  the  free-running  value,  |/lo|,  can  be  obtained  through  the 
cubic  equation: 

al-Yal-aL-Y^  =  0  (15) 

7n 

where  ai,  =  \Ai,IAo\^  7„  is  the  field-induced  decay  rate  due  to  the  differential  gain  term  for 
the  free-running  laser  [13]  and  the  factor,  F,  is  given  by: 

F  =  2^^^^  cos((?i>L).  (16) 

Here,  Or  is  the  radial  frequency  of  the  carrier-field  resonance,  0^  =  7c7n  +  7»7p?  where  7p 
is  the  field-induced  decay  rate  due  to  the  nonlinear  gain  term  for  the  free-running  laser 
[13],  and  ^  =  7;|>l;|/7c|i4o|  is  the  injection  parameter  [8].  The  injected  signal  also  modifies 
the  threshold  for  oscillation.  The  pump  parameter,  Jl  —  iJl^d)f'ysNL,  —  1,  of  the  locked 
oscillator  is  related  to  the  free-running  pump  parameter,  J  =  {Jfed)f%N  -  1  [8],  by: 

j  ^ _ [1  -2^cos((^L)/gL]qi>/ _ 

^  1  -  7sJ[2^  cos{(l>L)laL  -  7p{al  -  l)/7c]/7n 
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Alternatively,  |Ao|  and  \Ai\  are  related  through  an  equation  based  on  the  injection- locked 
parameters: 


2  2U ')nL 

ni-2U^s 

where  the  factor  =  (7c  -  2U)-(nL  +  7*7pL  is  used  for  convenience. 


(18) 


We  will  illustrate  the  effects  of  injection  locking  in  semiconductor  lasers  using  the  exper¬ 
imentally  measured  parameters  of  a  commercially  available,  index-guided,  GaAs/AlGaAs 
quantum-well  laser  diode.  The  relevant  parameters  for  this  laser,  along  with  a  reference  list 
of  key  symbols,  are  summarized  in  Table  1.  For  the  calculations  described  in  this  paper 
we  will  set  J  =  2/3  where  experimental  measurements  have  been  consistent  with  the  model 
[11].  Figure  2  is  a  mapping  of  the  operating  point  characteristics  in  the  detuning-frequency 
versus  injection-level  plane.  The  injection-locking  range  is  bounded  by  the  lines  of  constant 
phase  at  (t>L  =  ?r/2,  in  the  positive  frequency  range,  and  <t>L  =  —  arctan  at  negative  offset 
frequencies.  Shown  in  Figure  3(a)  are  lines  of  constant  phase  at  <l>i,  =  arctan  6,  the  zero  de¬ 
tuning  axis,  <i>L  =  0.67  «  0.5  arctan  6,  and  <f>L  =  0.  The  =  0  line  is  barely  distinguishable 
from  the  lower  injection-locking  boundary.  The  upper  bound  on  the  locking  range  is  also 
a  line  of  constant  amplitude,  ol  =  1.  When  the  injection  fields  are  small  compared  to  the 
output  field,  i  ~  0.05,  the  lines  of  constant  <j)L  are  nearly  straight  because  the  circulating 
field  is  essentially  unchanged  by  the  injection.  At  higher  injection  levels,  the  lines  of  constant 
phase  at  =  7r/2  and  at  <f>L  =  arctan  b  remain  straight  while  the  lines  at  =  0.67,  <i>L  =  0, 
and  4>l  =■  —  arctan  bend  over  as  \Al\  begins  to  depend  on  |A/|.  For  the  laser  under  con¬ 
sideration,  at  J  =  2/3  the  linearized  treatment  of  g  begins  to  break  down  for  ai  >  2,  when 
7nL  and  ')pL  become  a  sizeable  fraction  of  7c.  In  addition,  the  applicability  of  a  single-mode 
model  must  be  questioned  for  large  detunings,  on  the  order  of  the  longitudinal  mode  spacing. 
The  laser  with  the  parameters  listed  in  Table  1  has  a  mode  spacing  of  «  86  GHz. 


Before  continuing,  we  wish  to  consider  one  limiting  case  of  a  semiconductor  laser  under 
strong  optical  injection.  It  is  often  true  in  semiconductor  lasers  that  7^  7c  while  7p£,  is 
somewhat  larger  but  of  the  same  order  as  7„£,.  When  the  free-running  laser  is  biased  just 
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Parameter 

Symbol 

Reference/ Value* 

Polarization  of  gain  medium 

p 

Equation  1 

Gain 

9 

Equation  1 

Linewidth  enhancement  factor 

b 

4 

Differential  gain  parameter 

9n 

Reference  [13] 

Nonlinear  gain  parameter 

Sp 

Reference  [13] 

Intracavity  field 

A 

Equation  2 

Cavity  decay  rate 

Ic 

2.4  X  10“  s-i 

External  modulation  Field 

Ai 

Equation  2 

Injection  coupling  rate 

V 

Equation  2 

Carrier  Density 

N 

Equation  3 

Spontaneous  carrier  relaxation  rate 

Is 

1.458  X  10®  s-i 

Steady-state  injection  current  density 

J/ed 

Equation  3 

Modulation  offset  frequency 

n 

Equations  2  and  3 

Free-running  oscillation  frequency 

UJQ 

Equation  2 

Free-running  steady-state  field 

Aq 

Reference  [13] 

Free-running  steady-state  carrier  density 

N 

Reference  [13] 

Free-running  pump  parameter 

J 

Reference  [11] 

Field-induced  carrier  decay  rate 

In 

2J  X  10®  s-i 

Nonlinear  gain  decay  rate 

7p 

3.6  J  X  10®  s-i 

Relaxation  resonance  angular  frequency 

Qr 

(7c7n  +  Islpf'^ 

Total  carrier  relaxation  rate 

Ir 

75  “b  7n  "b  7p 

External  locking  field 

Ai 

Equation  7 

Locked  oscillation  frequency 

<^L 

Equation  7 

Modulation  current  amplitude 

Jm 

Equation  8 

Locked  steady-state  field  amplitude 

|AlI 

Equation  9 

Locked  Photon  Density 

Sl 

Equation  9 

Locked  steady-state  phase 

4>l 

Equation  11 

Locked  steady-state  carrier  density 

Nl 

Equation  12 

In-phase  locking  parameter 

U 

r)(\Ai\l\AL\)c.os(t>L 

Quadrature  locking  parameter 

V 

■n{\Ai\l\AL,\)  sin 4>l 

Normalized  steady-state  amplitude 

O-L 

\Al\I\Ao\ 

Injection  parameter 

e 

■q\Ai\hc\AQ\ 

Locked  Pump  parameter 

Jl 

Equation  17 

Normalized  modulation  field  amplitude 

a 

\A\I\Al\-1 

Phase  of  the  modulation  field 

Equation  23 

Normalized  modulation  carrier  density 

ft 

N/Nl  -  1 

Normalized  modulation  current 

dm  L 

Locked  field-induced  carrier  decay  rate 

InL 

Inal 

Locked  nonlinear  gain  decay  rate 

IpL 

Ipol 

Reference  angular  frequency 

Q.L 

[(7c  -  2[/)7„l  +  IslpLf^'^ 

Locked  decay  rate 

ItL 

Is  +  7nL  +  7pL  + 

TABLE  1.  Relevant  symbols  and  experimentally  determined  parameters  of  a  conventional, 
edge-emitting  laser  diode. 
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Injection  Parameter 


Figure  2:  Mapping  of  the  locking  and  stability  characteristics  of  a  semiconductor  laser  as 
a  function  of  the  frequency  offset  between  the  locking  and  free-running  laser  fields  and 
the  injection  parameter,  The  thick  solid  line  is  the  boundary  of  the  injection-locked 
region.  The  low-frequency  locking  boundary  is  a  constant  phase  line,  (f>i  =  —  arctan  6“^,  and 
the  high-frequency  boundary  is  the  line  where  (?i>L  =  7r/2  and  the  relative  circulating  field 
amplitude  is  ai,  =  1.  Also  shown  in  the  interior  of  the  locking  region  is  the  Hopf  bifurcation 
line  (thick  dashed  line)  separating  the  stable  and  unstable  operating  regimes  and  the  lines 
of  constant  amplitude  at  ai  =  1.1  (dotted  line),  =  1.5  (dot-dashed  line),  and  a^,  =  2 
(double  dot-dashed  line). 
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Injection  Parameter 


Figure  3(a).  Frequency  offsets. 


Injection  Parameter 

Figure  3(b).  Damping  rates. 


Figure  3:  Frequency  offsets  and  damping  rates  as  a  function  of  the  injection  parameter, 
along  lines  of  constant  phase.  In  (a),  the  frequency  offset  of  the  locking  frequency  and  eigen- 
frequency  of  the  non-zero  frequency  resonance  are  shown:  =  0  offset  (double  dot-dashed 

line)  and  eigenfrequency  (dotted  line),  <f)L  =  0.67  offset  (dot-dashed  line)  and  eigenfrequency 
(dashed  line),  and  <f>L  =  arctan  b  offset  (light  solid  line)  and  eigenfrequency  (dark  solid  line). 
In  (b),  the  damping  rates  for  the  eigenfrequencies  are  shown:  =  0  zero  (double  dot-dashed 

line)  and  non-zero  (dotted  line)  eigenfrequency  damping  rate,  =  0.67  zero  (dot-dashed 
line)  and  non-zero  (dashed  line),  and  <j)i  =  arctan  b  zero  (light  solid  line)  and  non-zero  (dark 
solid  line). 
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above  threshold  and  is  subjected  to  a  strong  injection  field,  so  that  ~  0,  then. 

IclnL 


u 


V 


2(7s  +7nL)’ 


(j^L  —  CJo  +  ^ 


^c'^nL 


2(7,  +  7nL)’ 


(19) 

(20) 


arctan 


2(75  +  7n£,)(^L  —  <*>0)  +  tndnL 
IclnL 


(21) 


We  will  use  these  relations  to  connect  our  results  on  the  modulation  characteristics  of  laser 
oscillators  with  the  results  of  optical  modulation  in  semiconductor  travelling-wave  amplifiers. 


All  modelling  in  the  following  sections  is  derived  from  this  basic  model,  these  definitions,  and 
the  assumption  that  both  the  modulation  field  and  the  modulation  current  can  be  modeled 
perturbatively  through  a  linearization  of  the  coupled  equations  about  the  operating  point. 
Here,  we  conclude  by  noting  that  the  weak  external  optical  probe  and  the  modulation  current 
are  assumed  to  be  independent  sources  which  perturb  the  operating  point.  In  the  linearized 
analysis  this  leads  to  independent  contributions  to  the  spectra  from  uncorrelated  sources. 


2.1.2  Linearized  Analysis 

With  these  relations  we  can  recast  the  coupled  equations  into  a  form  which  emphasizes 
the  different  dynamic  parameters.  The  full  nonlinear  equations  have  been  used  to  describe 
the  chaotic  dynamics  of  a  semiconductor  laser  under  optical  injection  [7].  We  perform  a 
perturbation  analysis  of  the  effects  of  the  weak  optical  probe  and  the  current  modulation  by 
assuming  that  the  normalized  field  magnitude,  a  =  |A|/|24/,|  —  1,  the  phase  variation  about 
^Li  the  normalized  carrier  density,  n  =  N/Nl  —  1,  are  all  much  less  than  unity.  The 

coupled  equations  are  then  linearized  as: 

—  =  -{U  + 'ypL)a  -  V ^  -k  77a,  cos  nt,  (22) 

dt  V  'P  /  27,  Jl 

^  =  {V  +  lrfpL)a  -  U ^  -  b— — -  T/a.  sinHt,  (23) 

dt  27,  Jl 
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^  =  27,X>cosfif  -  27,Jl(1 - -  (74  +  7„i)n.  (24) 

Here,  Jm  =  Jml^^d'ysNL  is  the  normalized  modulation  current  magnitude  and  we  have 
assumed  that  a,-  =  {Ail is  a  real  quantity  so  that  we  do  not  have  to  retain  an 
extraneous  phase  factor  through  the  calculation. 

Considerable  insight  about  the  modifications  in  laser  characteristics  can  be  obtained  by  per¬ 
forming  an  eigenvalue,  or  stability,  analysis,  of  the  coupled  equations  without  the  perturbing 
modulation  and  noise  terms.  The  resonance  frequency  and  bandwidth  can  be  approximated 
by  finding  the  complex  eigenvalues  of  the  coupled  equations  [7].  For  the  free-running  laser, 
the  eigenvalues  are  determined  by  the  equation: 

-in{-n^  +  nl-in-tr)  =  o,  (25) 

where  7r  =  la  +  7n  +  7p  is  the  field-enhanced  decay  rate.  There  are  two  non-zero  eigen- 
frequencies  given  by  ±(1)?  “  associated  damping  of  7r/2.  The  third 

eigenfrequency  of  the  free-running  laser  is  zero,  with  zero  damping,  an  unphysical  result 
that  arises  from  the  linearization  procedure’s  handling  of  the  phase  diffusion  of  the  oscillat¬ 
ing  field.  Retaining  the  full  nonlinearity  of  the  coupled  equations  yields  the  proper  damping 
behavior. 

The  injected  field  modifies  the  eigenfrequencies  and  damping,  as  given  by  the  zeros  of  the 
determinant, 

D  =  {—iCl  -|-  U)[—fl^  +  (Ij  +  ')nL)U  — 

-f (fl^  —  -f  (—if)  -+•  7s  -f  (26) 

where  7rL  =  7s  +  7nL  +  7pL  +  U.  For  the  injection- locked  laser,  the  zero  offset  eigenfrequency 
is  now  with  respect  to  the  locking  optical  frequency,  not  the  free-running  frequency,  and  has 
positive  damping.  The  non-zero  eigenfrequencies,  and  associated  damping,  are  also  modified. 
In  all  cases  the  modifications  follow  a  complex  dependence  on  the  strength  of  the  locking  field 
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and  the  detuning  beteen  the  master  and  slave  lasers.  The  eigenvalue  solution  is  relatively 
simple  for  the  case  where  (I>l  =  0.  In  this  case,  V  =  0  and  the  zero  frequency  eigenvalue  has 
an  associated  damping  of  U.  The  relaxation  resonance  frequency  is  now  modified  so  that 

=  T'cTni  +  7s7pL  +  {is  -  IulW 


and  the  field-enhanced  decay  rate  is  now  7rL.  The  non-zero  eigenfrequencies  and  associated 
damping  follow  with  the  obvious  substitutions  into  the  free-running  formulas. 

The  damping  associated  with  the  non-zero  eigenfrequencies  can  become  negative,  indicating 
that  a  Hopf  bifurcation  to  unstable  dynamics  has  occured.  The  condition  for  the  occurance 
of  a  Hopf  bifurcation  is  the  solution  of  Equation  26  with  a  real  non-zero  eigenfrequency.  By 
separately  setting  the  resulting  real  and  imaginary  parts  of  Equation  26  equal  to  zero,  this 
requires  that: 

-  bV)Ql  +  {irL  +  U)[lrLU  -f  V{ln,L  +  V)]  +  {is  +  lnL){lrLU  -  V^)  =  0.  (28) 


Depending  upon  the  detuning  between  the  master  and  slave  oscillators,  and  the  operating 
point  of  the  slave  laser.  Equation  28  will  be  satisfied  for  zero,  one  or  two  injection  field  levels. 

Figure  2  shows  the  line  denoting  the  Hopf  bifurcation  in  the  detuning  versus  injection- 
parameter  plane.  At  relatively  low  injection  levels,  0.01  <  ^  <  0.1,  the  region  of  instability 
occupies  a  sizeable  fraction  of  the  injection-locked  range.  At  higher  injection  levels,  the  Hopf 
bifurcation  line  asymptotically  approaches  the  <f>L  =  7r/2,  ai  =  I  boundary,  and  most  of  the 
injection-locked  range  ha^  stable  dynamics.  For  this  laser  at  this  operating  point,  stability  is 
assured  when  ai,  >  1.1.  Increasing  the  value  of  b  relatively  expands  the  unstable  operating 
range  while  increasing  J  reduces  it. 

To  give  some  insight  into  these  solutions,  we  can  make  some  approximations  to  simplify  the 
formula.  If  7r  <  Dr,  as  is  often  the  case  in  operating  semiconductor  lasers,  and  U  and  V  are 
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of  the  same  order  cis  7r,  then  Ai,  ^  Aq  and  Equation  28  will  have  a  solution  when  bV—U  ~  7r. 
At  this  bifurcation,  there  will  be  a  resonant  frequency,  of  ~  fir-  At  high  injection 

levels,  the  laser  will  undergo  a  reverse  Hopf  bifurcation  back  to  stable  operating  conditions. 
This  occurs  when  2U{U^  +  V^)/ibV  -  f/)  «  assuming  that  \Al\  «  |Aol,  and  the  resonant 
frequency,  now  satisfies  ^1^52  ^  d"  d"  these  simplified  formulas,  it 

follows  that  when  <^l  <  arctanfe"^  the  laser  remains  injection  locked  with  stable  operating 
characteristics.  Numerically,  for  the  laser  under  study  we  find  that  the  range  of  dynamic 
instability  is  bounded  by  (f>i,  w  0.67  >  arctanft  Injection  at  the  free-running  frequency, 
where  the  phase  offset  is  =  arctan  6,  induces  strong  modifications  in  both  the  frequency 
and  damping  with  increasing  injection  strength.  The  first  Hopf  bifurcation  point  occurs 
when  r}\Ai\/\Ao\  byr/{b^  —  ^)  and  the  reverse  bifurcation  when  7/|A;|/|Ao|  «  y^{b'^  —  l)/2Qr- 
Further, 


o  /^'d-l  1A,| 


(29) 


showing  that  for  strong  injection  fields  the  eigenvalue  solution  is  proportional  to  the  injection 
field.  If  the  semiconductor  laser  has  a  free- running  bias  point  such  that  7r/flr  >  byV2. 
then  an  injection  field  can  destabilize  the  dynamics  in  only  a  relatively  narrow  range  near 
the  =  7r/2  boundary.  Because  of  the  relations  between  fir  and  7r,  this  operating  point 
usually  demands  that  7„  >  7c.  In  practice,  this  requires  the  use  of  an  external  cavity,  where 
the  assumption  of  single-mode  operation  of  this  analysis  is  generally  more  difficult  to  justify. 

Figure  3(a)  plots  the  non-zero  positive  eigenfrequency  and  the  associated  frequency  offset, 
i'L  ~  of  the  master  laser  with  respect  to  the  free-running  slave  laser  as  a  function  of 
(  along  the  constant  phase  lines  =  0,  0.67  and  arctan  6.  Along  the  <t>L  =  0  line,  the 
non-zero  eigenfrequency  vanishes  for  ^  >  0.24,  indicating  a  transition  from  underdamped  to 
overdamped  dynamics.  Damping  rates  associated  with  the  non-zero  eigenfrequencies  and  the 
injection-locked  oscillation  frequency  are  shown  in  Figure  3(b).  Negative  damping  indicates 
unstable  dynamics.  The  4>l  -  0.67  line  approximately  bounds  the  unstable  region.  When 
<f>L  =  0  and  the  dynamics  are  underdamped,  the  damping  at  the  locking  frequency  is  U . 
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Figure  4  shows  similar  plots  at  constant  values  of  ^  as  the  injection-locked  phase  of  the 
field  is  varied.  Again,  the  non-zero  eigenfrequency  and  the  master-slave  frequency  offset  are 
shown  in  Figure  4(a)  and  the  associated  damping  rates  in  Figure  4(b).  Note  the  zero  crossing 
of  the  offset  at  (I>l  =  arctan  b.  At  high  injection  levels,  the  damping  remains  positive  over 
essentially  the  entire  injection-locking  range.  Also,  though  the  non-zero  eigenfrequency  is 
smaller  than  the  free-running  value  for  a  range  of  <t>L,  it  vanishes  only  very  close  to  =  0 
at  the  higher  injection  levels.  The  shift  to  true  overdamped  operation  occurs  for  only  a  very 
narrow  range  of  offset  frequency  and  injection  level. 


When  the  free-running  laser  is  biased  close  to  threshold,  then  U  and  V  can  dominate  the  other 
dynamic  rates,  except  7^  at  high,  but  accessible,  injection  powers.  Using  Equations  19  to  21, 
and  considering  injection  at  the  free-running  frequency,  ujq,  the  damping  rate  associated  with 
the  injection-locked  oscillation  frequency  is  approximately  874  -|-7nL-  Note  the  similarity  to 
the  field-enhanced  carrier  decay  rate,  75-l-7„i,.  The  non-zero  eigenfrequency  is  approximately 
U,  with  an  associated  damping  rate  of  U.  Both  U  and  V  depend  on  the  ratio  jclnL/l^ils  + 
7„l)]  in  this  limiting  case. 


When  the  optical  and  current  modulation  sources  are  added,  the  spectral  characteristics  are 
determined  by  the  Fourier  Transform  of  Equations  22  to  24. 


jF(a) 


(— ifl  -I-  74  +  7nL)[~^(^  —  V)  +  U]  r]\Ai\ 
D  2|Al| 


+  c.c. 


(-in  +  U  +  bV){'fc  -  2U)  7nL  Jm  , 

^ - = - h  C.C. 

D  2Jr, 


(30) 


+  7.  +  7nL)[-(fi  -V  -  brjpL)  -  HU  +  ~jpL)] 
+{b  -  l)(7c  -  2t/  -  7pt)7n£,}||^  +  C.C. 


+4(>i'n  -w+  v){i,  -  2U)'^ 

Lf  2Jl 


+  C.C., 


(31) 


20 


Locked  Phase  (radians) 

Figure  4(b).  Damping  rates. 

Figure  4:  Frequency  offsets  and  damping  rates  as  a  function  of  the  locked  phase,  <f>L,  for 
constant  values  of  the  injection  parameter.  In  (a),  the  frequency  offset  of  the  locking  fre¬ 
quency  and  eigenfrequency  of  the  non-zero  frequency  resonance  are  shown:  ^  =  1.0  offset 
(double  dot-dashed  line)  and  eigenfrequency  (dotted  line),  ^  =  0.5  offset  (dot-dashed  line) 
and  eigenfrequency  (dashed  line),  and  if  =  0.2  offset  (light  solid  line)  and  eigenfrequency 
(dark  solid  line).  In  (b),  the  damping  rates  for  the  eigenfrequencies  are  shown:  f  =  1.0  zero 
(double  dot-dashed  line)  and  non-zero  (dotted  line)  eigenfrequency  damping  rate,  f  =  0.5 
zero  (dot-dashed  line)  and  non-zero  (dashed  line),  and  f  =  0.2  zero  (light  solid  line)  and 
non-zero  (dark  solid  line). 
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(32) 


:r(fi)  =  \(-isi  +  u){-m  +  u +  ^rL)  +  v(v■^■ln^L)\^  +  c■c■ 

-(-,,  -2U-  -V)  +  f/]^  +  c.^. 

Here,  ^(x)  is  the  Fourier  Transform  of  x,  D  is  the  determinant  given  by  Equation  26  and  c.c. 
is  the  complex  conjugate  which  corresponds  to  a  frequency  component  equally  and  oppositely 
shifted  with  respect  to  dl  for  the  optical  frequency  spectra,  and  with  respect  to  zero  for  the 
carrier  density  and  optical  power  spectra. 

For  optical  modulation,  the  linearized  analysis  predicts  that  the  output  field  of  the  injection- 
locked  laser  consists  of  spectral  components  at  the  locked  oscillating  frequency,  the 
injection  frequency,  ul  +  and  the  four-wave  mixing  parametric  frequency,  ujl  —  From 
the  view  point  of  regenerative  amplification  and  four-wave  mixing,  it  can  be  mathematically 
described  by 

E{t)  =  +  c.c.  (33) 


with 


A{t)  =  Al  + Are~'^^  A  (34) 

where  Ay  and  Aj  are  the  complex  amplitudes  of  the  regeneratively  amplified  and  four-wave 
mixing  fields,  respectively.  Alternatively,  it  can  be  described  by 

E{t)  =  El[\  -t-  crcos(nt  -  <(3i)]  cosfwit  +  6cos{Q.t  -  (^2)],  (35) 

as  a  field  undergoing  both  amplitude  and  phase  modulations,  where 


&  =  cre"^^  and  6  = 


(36) 


are  the  complex  amplitudes  of  the  amplitude  and  phase  modulations,  respectively.  The  two 
pictures  are  connected  through 


..  4-  A'jAl 


and  S  =  i 


AtA'^  —  AjAjj 


(37) 


or 


A 

Al 


and 


a  -\riS 
2 


(38) 


Thus  both  Ar  and  A j  have  contributions  from  both  the  amplitude  and  phase  modulations  of 
the  laser  field.  The  situation  is  similar  for  the  case  of  current  modulation,  where  Ar  can  be 
replaced  by  Ap  and  Aj  hy  An  to  represent  the  positive  and  negative  frequency  components 
of  the  sinusoidal  modulation,  and  the  same  formalism  can  be  used. 


As  is  the  case  with  the  free-running  laser,  the  key  component  of  the  four-wave  mixing  process 
is  the  carrier  pulsation  in  the  semiconductor  gain  medium  caused  by  the  beating  between 
the  mixing  optical  fields.  In  a  single-mode  semiconductor  laser,  the  spatial  gratings  are  not 
important  and  the  process  is  solely  contributed  by  the  temporal  pulsation  of  the  carrier 
density  at  the  beat  frequency  [25,  23].  In  the  linearized  analysis,  the  carrier  density  can  be 
expressed  as 


N{t)  =  NlA  Me"’”'  + 


(39) 


where  N\  is  the  amplitude  of  the  carrier  pulsation. 


To  reveal  the  physical  mechanisms  involved  in  the  process,  we  first  use  Equations  35  and  39 
with  Equations  30  to  32  to  obtain 


iVi  = 


7c  -  2f/  -  7p  Sl 

Is  +  7nL  -  r 


(40) 
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Clearly,  this  indicates  that  the  carrier  pulsation  is  caused  by  the  amplitude  fluctuation  of 
the  total  optical  field  in  the  laser  cavity.  When  — >■  0,  and  a  are  out  of  phase  by  tt. 
This  means  that  the  induced  carrier  pulsation  acts  to  reduce  the  amplitude  fluctuation  in 
the  laser  field  caused  by  the  optical  modulation. 


Using  Equations  34  and  40  with  Equations  30  to  32,  we  have 

Ar  i  f  ^  ~  —  e'07pL  r M 

H  ^  n  +  V  +  iU  V  2  7c  -  2C/  -  7pL  Sl  ^  Al)' 


A*j  i  /l+ib  9,1-  i9'^pL  TNi\ 

Al  ""  9-V  +  iU\  2  'fc-2U  --fpL  Sl  J  ' 


From  Equation  37,  we  further  obtain 

i{9  -V)-U 


i9  +  iU  +  ibV){9l-i9'fpL)TNi  rjAj 
Q2  _  (f/2  ^  y2)  +  2iU9  [(D- V  +  iU)(7c-2C/-7pL)  Sl  ^  Al  , 


8  = 


ii9  -V)-U 


Q2  _  (JJ2  +  ^  2iU9 


ib9  +  ibU  -  iV)i9l  -  in-fpi)  TNi  .rjAi 

(fi  -  V  +  ft/)(7c  -2U-  -fpL)  Sl  '  Al 


(42) 


(43) 

(44) 


These  relations  show  the  various  physical  components  of  contribution  to  the  regeneratively 
amplified  and  four-wave  mixing  signals,  as  well  as  those  to  the  amplitude  and  phase  modu¬ 
lations  of  the  field.  The  factor 


G{9)  = 


91  -  iQ'fpL  TNi 
7c  —  2t/  -  7pL  Sl 


(45) 


appears  in  all  four  relations.  It  is  the  gain  caused  by  the  carrier  pulsation.  This  can  be 
easily  seen  from  the  fact  that  when  the  nonlinear  gain  coefficient  is  negligible  so  that  7p  =  0, 
G  =  g„Ni.  In  a  semiconductor  medium  in  general,  and  a  free- running  laser  in  particular,  a 
change  of  gain  by  G  is  always  accompanied  by  a  change  of  refractive  index  by  bG.  Thus  an 
amplification  of  the  amplitude  modulation  by  G  is  accompanied  by  an  amplification  of  the 
phase  modulation  by  bG.  This  is  seen  in  Equations  43  and  44  when  U  =  V  =  0  and,  more 
clearly  in  Reference  [13].  Under  injection  locking,  the  quadrature  injection  term,  V,  changes 
the  amplitude  and  phase  coupling  characteristics  among  the  carriers,  the  circulating  optical 
field  and  the  injected  optical  field. 
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Both  amplitude  and  phase  modulations  also  have  contributions  from  the  externally  injected 
field.  In  contrast,  the  externally  injected  field  directly  contributes  only  to  the  regeneratively 
amplified  signal  at  the  same  frequency,  but  not  to  the  four- wave  mixing  signal  at  a  different 
frequency.  It  is  only  through  the  carrier  pulsation  that  the  external  modulation  field  induces 
the  four-wave  mixing  signal.  This  is  seen  by  using  Equations  32  and  39  to  obtain 


Ni  = 


(7c  -2U-  ^j,l)[U  -  ^•(^  -  V)]  Sl  yAi 
D  T  Al' 


(46) 


In  the  free-running  laser,  the  modulation  on  the  laser  field  at  ujq  by  the  gain  fluctuation  G 
and  the  accompanying  refractive-index  fluctuation  bG  at  generates  sidebands  at  wo  + 
and  Wo with  equal  efficiency  although  different  phases.  Thus  the  regeneratively  amplified 
signal  and  the  four-wave  mixing  signal  both  have  contributions  from  G  and  bG  with  equal 
magnitude  but  different  phases.  This  is  no  longer  true  in  the  injection-locked  laser  due  to  the 
different  resonant  denominators  in  Equations  41  and  42.  The  resonant  denominator  for  the 
four-wave  mixing  sideband  is  cancelled  by  a  similar  term  in  the  numerator  of  Equation  46. 
There  is  no  such  cancellation  in  the  regeneratively  amplified  sideband  so  that  it  will  contain 
a  resonance  term  at  negative  offset  frequencies  and  a  quenching  term  at  positive  offset 
frequencies. 

The  connection  between  optical  modulation  and  carrier  pulsation  is  direct  in  the  case  of 
current  modulation  of  the  carrier  density.  Now,  from  Equations  32  and  39,  we  have 


Ni  = 


(-in  -h  u){in  +  U  +  ipl)  +  v{v  +  inpi)  Jm 

D  2ed 


(47) 


The  solutions  for  the  positive  and  negative  optical  frequency  components,  Ap  and  A„,  re¬ 
spectively,  can  be  obtained  in  the  same  way  as  Ar  and  Aj: 


Al 


K 


i  +  b 


i-n  +  V-iUhr^L 


TNi 


2  i-m  U){-in  +  U  +  Jpl)  +  V{V  +  67pl)  Sl  ’ 

-i  +  b _ {n  +  v-  iu)jnL _ 
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Similarly,  we  can  obtain  the  solutions  for  the  amplitude  and  phase  modulation  coefficients: 

. _ Hn  +  f/  +  6y)7nL _ w  (50) 

^  {—iO,  +  U){—iCl  +  U  +  7pl)  -\-V{V  b^ph)  Sl 

.  _  _ {-m^hU  +V)'ir.L _ m  (51) 

^  "  (-if]  +  t/)(-in  +  t/  +  7pL)  +  y(v  +  hpL)  ■ 

In  Equations  48  to  51  the  combination  of  parameters  '^tiL^NxISl  =  is  common  and 

is  the  gain  due  to  the  carrier  pulsation.  Because  the  current  modulates  the  carrier  density 
directly,  the  gain  term  here  is  simpler  than  G  derived  for  the  case  of  optical  modulation.  As 
was  the  case  with  optical  modulation,  the  quadrature  injection-locking  term,  V ,  removes  the 
simple  ratio  of  b  between  the  phase  and  amplitude  modulations  that  exist  in  the  case  of  the 
free-running  laser. 

2.1.3  Modulation  Spectra 

The  spectral  characterization  of  semiconductor  lasers  uses  a  variety  of  techniques.  Techniques 
which  use  a  narrowband  optical  filter,  such  as  an  optical  spectrum  analyzer  or  an  optical 
heterodyne  configuration,  can  be  used  to  map  out  the  squared  amplitude  of  the  modulation 
sidebands,  |Arp,  |A/P,  |App  and  |A„p.  The  sidebands  contain  differing  combinations  of 
the  amplitude  and  phase  modulation  characteristics.  Spectral  techniques  which  analyze  the 
spectrum  of  the  signal  generated  by  a  fast  photodiode  when  illuminated  by  the  semiconductor 
laser  contain  only  information  about  the  amplitude  modulation  characteristics.  When  the 
detected  signal  from  a  fast  photodiode  is  fed  into  an  RF/microwave  spectrum  analyzer, 
the  spectrum  of  the  amplitude  modulation  coefficient,  |d-p,  can  be  measured.  The  spectra 
of  other  quantities  of  potential  interest,  such  as  |^p  and  |A^iP,  can  be  determined  from 
measurements  of  the  optical  sidebands  and  the  amplitude  modulation  coefficient. 

Here,  we  present  representative  spectra  for  optical  and  current  modulation  of  an  injection- 
locked  laser  operating  in  the  region  of  stable  dynamics.  The  operating  parameters  of  the  laser 
are  listed  in  Table  1  and  J  =  2/3  for  these  calculations.  Comparisons  will  be  made  for  two 
cases.  First,  the  master  laser  has  the  same  optical  frequency  as  the  free-running  slave  laser. 
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corresponding  to  =  arctani,  and  ^  is  varied:  0.005,  0.2,  0.5  and  1.0.  The  first  two  values 
of  ^  roughly  bracket  the  region  of  unstable  operation.  Second,  ^  =  0.2  and  <f>L  is  varied: 
0,  0.67  and  arctan  b.  At  this  injection  strength,  when  <?!>l  =  0  there  is  a  frequency  offset  of 
—22.4  GHz  between  the  locked  and  free-running  oscillation  frequencies.  For  (f>L  =  0.67  the 
offset  is  -14.9  GHz,  and  for  <I>l  =  arctan  b  the  offset  is  0.  Also  shown  in  this  series  are  the 
spectra  for  the  free-running  laser.  Note  that  the  ^  =  0.2,  <j)L  =  arctan  b  curves  are  common 
to  both  cases  to  assist  in  comparisons  betweem  them.  Spectra  of  optical  field  modulation 
characteristics  will  use  the  free-running  slave-laser  frequency  as  the  zero  frequency  reference 
for  the  calculated  spectra  that  follow. 


Regenerative  amplification  and  four-wave  mixing  sidebands  generated  by  weak  optical  mod¬ 
ulation  of  an  injection-locked  semiconductor  laser  will  have  the  following  spectral  character¬ 
istics: 
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The  regenerative  amplification  sideband  spectrum  is  asymmetric  with  respect  to  detuning 
about  the  injection-locked  operating  frequency.  This  asymmetry  is  extreme  under  a  broad 
range  of  operating  conditions,  with  a  strong  resonance  at  negative  offset  frequencies  and  a 
quenching  at  positive  offset  frequencies.  The  four-wave  mixing  sideband  spectrum,  on  the 
other  hand,  is  symmetric  with  relatively  weaker  resonances  equally  and  oppositely  offset 
from  the  injection-locked  frequency.  The  pole  at  zero  detuning  for  the  free-running  laser 
is  removed  under  injection  locking,  though  the  spectral  characteristics  around  zero  will  be 
obscured,  to  some  degree,  by  the  spectrum  of  the  master  laser.  Away  from  zero  detuning,  the 
resonant  carrier-field  coupling  may  be  destabilized,  leading  to  oscillation  at  the  negative,  and 
sometimes  the  positive,  eigenfrequency  over  a  range  of  operating  conditions.  Equations  52 
and  53  are  certainly  not  a  good  description  of  experimentally  observed  spectra  in  the  unstable 
operating  regime.  In  that  case  the  full  nonlinear  equations  must  be  retained. 
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Figure  5  plots  the  regenerative  amplification  spectra  for  different  values  of  ^  and  the  signal 
value  is  (7c|/lr|/27/|A,|)^.  The  frequency  axis  is  taken  relative  to  the  free- running  operating 
point.  A  weak  asymmetry  in  favor  of  negative  offset  frequencies  in  the  free- running  laser 
becomes  quite  strong  as  the  injection  level  is  increased. 

Figure  6  plots  the  four- wave  mixing  spectra  for  different  values  of  ^  and  the  signal  value 
is  (7c|A/|/27/|A,|)^.  Again,  the  frequency  axis  is  relative  to  the  free-running  operating  point. 
The  frequency-dependent  term  in  the  numerator  is  small  compared  to  the  7c7nL  term  so 
that  the  shape  of  the  spectrum  is  determined  by  the  |Dp  term  in  the  denominator.  Here, 
the  three  resonances,  or  eigenfrequencies,  remain  clearly  developed  in  a  symmetric  spectrum 
when  the  injection  is  at  the  free-running  frequency,  Figure  6(a).  When  the  injection  is 
frequency  offset.  Figure  6(b),  the  increased  damping  obscures  the  resonances  offset  from  the 
locking  frequency. 

The  amplitude  modulation  spectrum  is,  in  general,  asymmetric  about  the  injection-locked 
frequency: 


vAi  ^  [ii^+hs+inifm-vr+u^] 

Al  w 


(54) 


Figure  7  plots  the  spectra  for  different  values  of  (  and  4>l'i  the  signal  value  is 
\^^{lc\Ao\l2r}\Ai\Y.  Only  in  the  free-running  and  <^i,  =  0  cases  are  the  amplitude  mod¬ 
ulation  spectra  symmetric.  The  importance  of  V  in  this  asymmetry  is  clear  in  Equation  54. 
Also  note  the  strong  quenching  of  the  amplitude  modulations  near  the  locking  frequency. 
The  spectral  width  of  this  quenching  region  is  primarily  determined  by  the  field-enhanced 
carrier  decay  rate,  7s  -f  7ni,. 

Because  |Ap(-H)P  =  |An(fl)P,  the  optical  spectra  of  the  current  modulation  characteristic 
are  given  by 

|ApP  _  1  +  9  [(n  -  Vf  +  U^]{lc  -  mWnL  .55^ 

|A,P”  4  w  j/ 
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Figure  5;  Spectrum  of  the  regenerative  amplification  sideband  characteristic  due  to  a  weak 
optical  modulation,  the  signal  value  is  (7c|.4r|/27?|/l,|)^.  The  frequency  offset  is  with  respect 
to  the  free-running  optical  frequency,  (a)  Injection-locking  field  at  the  free-running  optical 
frequency  with  ^  =  0.005  (dashed),  ^  =  0.2  (solid),  .f  =  0.5  (dotted),  and  ^  =  1.0  (dot- 
dashed).  (b)  Injection-locking  field  at  =  0.0  (dashed)  and  at  (f  =  0.2  with  Ol  =  0.0 
(dot-dashed).  6l  =  0.67  (dotted),  and  01  =  arctan  5  (solid). 
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Figure  6:  Spectrum  of  the  four-wave  mixing  sideband  characteristic  due  to  a  weak  optical 
modulation,  the  signal  value  is  (7c|.4/|/27;|.4i|)^.  The  frequency  offset  is  with  respect  to 
the  free-running  optical  frequency,  (a)  Injection-locking  field  at  the  free-running  optical 
frequency  with  ^  =  0.005  (dashed),  ^  =  0.2  (solid),  ^  =  0.5  (dotted),  and  ^  =  1.0  (dot- 
dashed).  (b)  Injection-locking  field  at  ^  =  0.0  (dashed)  and  at  =  0.2  with  ol  =  0.0 
(dot-dashed),  Ol  =  0.67  (dotted),  and  Si  =  arctanft  (solid). 
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Figure  7:  Spectrum  of  the  amplitude  modulation  characteristic  due  to  a  weak  optical  modu¬ 
lation.  The  signal  value  is  |dp(7c|/lo|/‘-^Mi|)^-  (a)  Injection-locking  field  at  the  free-running 
optical  frequency  with  =  0.005  (dashed),  ^  =  0.2  (solid),  ^  =  0.5  (dotted),  and  ^  =  1.0 
(dot-dashed),  (b)  Injection-locking  field  at  =  0.0  (dashed)  and  at  ^  =  0.2  with  pi  =  0.0 
(dot-dashed).  6i  =  0.67  (dotted),  and  61  =  arctan  6  (solid). 


Except  for  the  cases  of  the  free-running  and  <f>L  =  0  operation,  the  spectra  are  asymmetric 
with  respect  to  the  oscillation  frequency  and,  like  the  regenerative  amplification  spectra,  have 
a  stronger  resonance  at  negative  offset  frequencies.  Figure  8  plots  the  spectra  for  different 
values  of  ^  and  <I>l]  the  signal  value  is  \Ap/Al\'^ /{J^/2jy.  The  frequency  axis  is  relative  to 
the  free-running  operating  point.  Injection  locking  breaks  the  spectral  symmetry  exhibited 
by  the  free-running  spectrum.  Like  the  regenerative  amplification  spectrum,  the  modulation 
is  enhanced  for  negative  offset  frequencies. 

The  amplitude  modulation  spectra  are  symmetric  with 


+  (U  +  bvn^,  -  2U)WnL  Jm' 


(56) 


Figure  9  plots  the  amplitude  modulation  under  current  modulation  for  different  values  of  ^ 
and  (t>L‘,  the  signal  value  is  |(7p/( Only  positive  frequencies  are  shown  because  of 
the  inherent  spectral  symmetry.  Strong  optical  injection  increases  the  non-zero  resonance 
frequencies  and  quenches  the  resonances.  In  Figure  9(b)  we  see  how  either  modulation 
bandwidth  or  flatness  of  response  can  be  emphasized  depending  on  the  injection  conditions. 
The  bandwidth  can  be  more  than  doubled  by  injection  signals  which  are  weak  relative  to 
the  laser  output. 


2.2  LINEARIZED  ANALYSIS  INCLUDING  FEEDBACK 


Here,  we  augment  the  analysis  described  above  by  including  the  effects  of  optical  feedback  on 
the  injection-locked  operating  characteristics.  Optical  feedback  can  be  included  in  Equation  2 
by  adding  a  term  rA{t  —  T)  where  r  is  the  coupling  parameter  for  the  strength  of  the  feedback 
and  T  is  the  round-trip  time  delay  for  the  back-reflected  output.  It  is  well  known  that  external 
feedback  can  induce  nonlinear  dynamics  in  a  laser  diode  [27].  Therefore,  this  analysis  will 
only  be  appropriate  for  feedback  which  is  weak  compared  to  the  external  optical  injection 
and  where  the  injection  is  not  inducing  an  operating  point  close  to  the  Hopf  bifurcation 
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Figure  8:  Spectrum  of  the  optical  field  positive  sideband  characteristic  due  to  a  weak  current 
modulation.  The  signal  value  is  \Apl Al\^ ■  The  frequency  offset  is  with  respect 
to  the  free-running  optical  frequency,  (a)  Injection-locking  field  at  the  free-running  optical 
frequency  with  i  =  0.005  (dashed),  ^  =  0.2  (solid),  {  =  0.5  (dotted),  and  i  =  1.0  (dot- 
dashed).  (b)  Injection-locking  field  at  ^  =  0.0  (dashed)  and  at  4  =  0.2  with  <i>L  =  0.0 
(dot-dashed),  (t>L  =  0.67  (dotted),  and  4>l  =  arctan6  (solid). 
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Figure  9:  Spectrum  of  the  amplitude  modulation  characteristic  due  to  a  weak  current  mod¬ 
ulation.  The  signal  value  is  \a\'^ l{Jjnl'2JY .  (a)  Injection-locking  field  at  the  free-running 
optical  frequency  with  ^  =  0.005  (dashed),  ^  =  0.2  (solid),  ^  =  0.5  (dotted),  and  ^  —  1.0 
(dot-dashed),  (b)  Injection-locking  field  at  ^  =  0.0  (dashed)  and  at  =  0.2  with  Sl  =  0.0 
(dot-dashed),  <i>L  =  0.67  (dotted),  and  6l  =  arctan6  (solid). 
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separating  the  regions  of  stable  and  unstable  dyanmics.  We  proceed  following  the  same  line 
of  analysis  as  used  above.  The  notation  will  be  kept  as  consistent  as  possible.  The  addition 
of  the  feedback  further  modifies  the  operating  point  of  the  laser  so  that  now: 


UL  -uo  =  -hU  +  V  -  br  cos{ijjlt)  +  r  sin(u;£,T).  (57) 

Again,  U  =  'r}\Ail Ai\  cos  <t>i  is  the  in-phase  injection  term  and  V  =  sin<?!>L  is  the 

quadrature  term.  In  the  absence  of  feedback,  the  gain  defect  at  the  injection-locked  operating 
point  is  2U  and  the  frequency  offset,  ul  —  Uc  —  bTgif2,  where  gi,  is  the  steady-state  gain 
at  the  injection-locked  operating  point,  is  V,  as  described  above.  The  feedback  causes  an 
additional  shift  such  that  the  gain  defect  is  now  2U  +  2r  cos{ult)  =  2U'  and  the  locked  phase 
is  modified  so  that  the  frequency  shift  is  now  V  -|- sin(a;i:,r)  =  V .  The  relation  between  \Ai\ 
and  the  free- running  value,  |Ao|,  can  be  obtained  through  the  cubic  equation: 


al  -  Y'al  _  ai  -  r  ^  =  0 

^  ^  In 

where  ai  =  IAl/AqI,  the  factor,  Y\  is  given  by: 


(58) 


Y'  =  2 


7c7n 


^  cos(^l)  H - cos{ult) 

7c 


(59) 


The  injected  signal  and  feedback  also  modify  the  threshold  for  oscillation.  The  pump  pa¬ 
rameter  of  the  locked  oscillator  is  related  to  the  free-running  pump  parameter  by: 


_ (1  -  2U')alj _ 

\-l,J[2U'-i,{al-l)hc]hn 


(60) 


Alternatively,  |Ao|  and  \Ai\  are  related  through  an  equation  based  on  the  injection- locked 
parameters: 


_  flj  -H  2f/^7nL 
-  2t/'7 
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where,  now,  the  factor  =  (7c  —  2/7')7„£,  +  ')a'1pL-  The  coupled  equations,  including  weak 
feedback,  are  then  linearized  as: 


—  =  — (t/' +  7p£,)a  +  rcos(w£,T)a(t  —  r)  —  yV  —  rsin(a;i:,T)^(i  —  r) 

dt 


(7c  -  2C/')7nL  -  .  r.. 

-n  +  i}o,i  cos  Sit, 


(62) 


dt 


=  {V  +  lrjpL)o-  +  f  sm{uLT)a{t  —  r)  —  U'^  +  r  cos{u}it )^{t  —  t ) 

hc-2U')'fnL.  (63) 
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-h  —  r]ai  sin  fit. 


7pL 


(64) 


lc-2U' 

Here,  X.  =  Jml’^.ed^aNi  is  the  normalized  modulation  current  magnitude  and  we  have 
assumed  that  ai  =  {Ai/AL)e~"^‘‘  is  a  real  quantity  so  that  we  do  not  have  to  retain  an 
extraneous  phase  factor  through  the  calculation,  as  before. 


The  spectral  characteristics  are  determined  by  the  Fourier  Transform  of  Equations  62  to  64. 

’2\Al\ 


_  .  _  i-in  +  i3  +  inL)[-i{n-vn  +  u"]v\M  ,  ^ 


= 


D 


{-in  +  U"  +  bV"){ya-2U')'ynLJm  ,  _ 

4-^ —  ~  ’1“  C.C., 

D  2J'r 


^{(-ifl  +  75  +  7nL)[-(fl  -  V"  -  IrfpL)  -  i{U"  +  7pl)] 

■n\A\ 

+{b  -  0(7c  -  2t7'  -  7pi)7nL}^i^  +  c.c. 


+i(ii.f!  -  bU"  +  V"){1,  -  2U'y!^  +  C.C. 


(65) 


2J'r 


(66) 


IsJv 


J^(n)  =  |(-tn  +  £/")(-if2+(/"  +  7,L)  +  V"(V"'  +  hn))^  +  cc 
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-h.-2£/'-7,i)l-i(n-n  +  C'"l;^  +  c‘:-  (6^) 

Here,  T{x)  is  the  Fourier  Transform  of  x,  U"  =  f/'-e*“‘rcos(a>LT),  V"  =  F'  +  e'^VsinKr), 
D  is  the  determinant  given  by: 

D  =  +  +  9i^{-is+lnL)U" 

+  (f)^  —  iCtjpL)bV"  +  (  — jfi  +  7a  +  InhW'^i  (6^) 

where  7rL  =  7a  +  7nL  +  7pL  +  U",  and  c.c.  is  the  complex  conjugate  which  corresponds  to 
a  frequency  component  equally  and  oppositely  shifted  with  respect  to  ul  for  the  optical 
frequency  spectra,  and  with  respect  to  zero  for  the  carrier  density  and  optical  power  spectra. 
The  equations  above  are  given  in  the  same  format  as  was  used  above  and  in  Reference  [12]. 
and  reduce  to  those  equations  when  r  =  0. 


Regenerative  amplification  and  four-wave  mixing  sidebands  generated  by  weak  optical  mod¬ 
ulation  of  an  injection-locked  semiconductor  laser  will  have  the  following  spectral  character¬ 
istics: 
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As  before,  the  regenerative  amplification  sideband  spectrum  is  asymmetric  with  respect  to 
detuning  about  the  injection-locked  operating  frequency.  This  asymmetry  is  extreme  under 
a  broad  range  of  operating  conditions,  with  a  strong  resonance  at  negative  offset  frequencies 
and  a  quenching  at  positive  offset  frequencies.  The  four-wave  mixing  sideband  spectrum,  on 
the  other  hand,  remains  symmetric  with  relatively  weaker  resonances  equally  and  oppositely 
offset  from  the  injection-locked  frequency.  The  pole  at  zero  detuning  for  the  free-running 
laser  is  removed  under  injection  locking,  though  the  spectral  characteristics  around  zero 
will  be  obscured,  to  some  degree,  by  the  spectrum  of  the  master  laser.  Away  from  zero 
detuning,  the  resonant  carrier-field  coupling  may  be  destabilized,  leading  to  oscillation  at 
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the  negative,  and  sometimes  the  positive,  eigenfrequency  over  a  range  of  operating  conditions. 
Equations  69  and  70  are  certainly  not  a  good  description  of  experimentally  observed  spectra 
in  the  unstable  operating  regime.  In  that  case  the  full  nonlinear  equations  must  be  retained. 
The  amplitude  modulation  spectrum  is,  in  general,  asymmetric  about  the  injection-locked 
frequency: 


{n^  +  i%+ir.L?m-vy  +  u"^] 

w 


(71) 


The  optical  spectra  of  the  current  modulation  characteristic  are  modified  from  the  feedback- 


free  case  to: 


_  1  +  p  f(n  -  vr  -h  Jm" 

4  w  JC 


(72) 


As  before,  except  for  the  cases  of  the  free-running  and  <i>L  =  ^  operation,  the  spectra  are 
asymmetric  with  respect  to  the  oscillation  frequency  and,  like  the  regenerative  amplification 
spectra,  have  a  stronger  resonance  at  negative  offset  frequencies.  The  amplitude  modulation 
spectra  are  symmetric  with 


|dP  = 


+  jU"  +  fcl/^0^](7c  -  mWnL  Jn 


|Z>P 


r 

'^L 


(73) 


Again,  the  spectra  revert  to  the  formulas  above  and  in  Reference  [12]  when  r  —  0. 


2.3  NUMERICAL  SIMULATION  OF  INJECTION-LOCKED  OPERATION 


Because  the  coupled  field  and  carrier  density  equations  are  nonlinear,  in  general,  they  must 
be  solved  numerically.  Here,  we  describe  the  techniques  used  to  solve  the  full  nonlinear 
equations.  By  comparison  of  the  numerical  solutions  of  the  full  nonlinear  equations  with  the 
analytical  solutions  of  the  linearized  approximation,  we  can  show  where  the  linearization  is 
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appropriate.  Some  care  must  be  taken  in  the  numerical  integration  of  the  coupled  equations. 
Equation  2  is  a  stochastic  nonlinear  differential  equation  because  it  contains  the  stochastic 
noise  source  term  as  well  as  the  deterministic  external  field  source  term.  The  numerical 
integration  of  differential  equations  is  typically  done  using  a  Taylor-series  type  expansion 

[28] .  One  wants  to  determine  the  quantity  X{t  +  h),  at  the  time  t  h,  knowing  X{t).  The 
Taylor-series  expansion  parameter  is  h.  There  are  several  standard  solution  algorithms  for 
solving  deterministic  differential  equations.  For  example,  the  standard  4th  order  Runge- 
Kutta  integrations  scheme  is  available  in  commercial  software  packages  such  as  Mathcad. 
It  is  called  4th  order  because  the  expansion  is  accurate  to  4th  order  in  the  parameter  h. 
However,  it  can,  at  best,  only  be  made  accurate  to  1st  order  in  h  for  a  stochastic  source 

[29] .  Therefore,  the  stepsize,  h,  must  be  much  smaller  with  a  stochastic  term  in  the  equation 
than  for  a  purely  deterministic  equation.  Different  algorithms  are  required  for  the  numerical 
integration  of  stochastic  differential  equations  with  accuracy  higher  than  first  order  in  h 
[30,  31].  An  algorithm  accurate  to  2nd  order  in  h  was  implemented  in  Fortran  for  comparison 
with  the  modified  Runge-Kutta  algorithm  which  could  be  implemented  with  Mathcad.  Based 
on  this  comparison,  both  algorithms  gave  the  same  solutions  over  the  range  of  parameters 
discussed  in  this  report.  We  have  found  that  as  long  as  the  stepsize  is  sufficiently  small, 
typically  77^  is  appropriate,  and  the  integration  time  is  sufficiently  long,  typically  on  the 
order  of  10®  steps,  the  numerical  integration  produced  consistent  results.  These  results 
agreed  with  the  analytical  solutions  of  the  linearized  analysis  in  the  appropriate  limit. 


For  numerical  simulation,  the  following  coupled  equations  [7]  are  used: 


da 

dt 


7c7n- 


IsJ 


-  7p(2a  -1-  a^) 
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(74) 

(75) 


^  =  — 7,n  —  7„(1  +  a)^n  —  7i  j(2a -h  a^) -h  j(2a -h  a^)(l -H  a)^,  (76) 

dt  7c 

where  7^  is  the  spontaneous  carrier  relaxation  rate,  7„  is  the  differential  carrier  relaxation 
rate,  7p  is  the  nonlinear  carrier  relaxation  rate,  b  is  the  linewidth  enhancement  factor,  and 
J  =  {Jfed--fsNo)hsNo  with  e  being  the  electronic  charge  and  d  the  active  layer  thickness  of 
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the  laser.  The  normalized  noise-source  parameters  Fa  and  are  related  to  the  spontaneous 
emission  rate,  i2jp,  of  the  laser  through  the  following  relations: 

(F.{t)F,(f))  =  {Fi(i)Ft(t'))  =  ^^S(t-t'),  (77) 

{F.(t)Ft(t'))  =  (Ft{i)F.(t'))=(S.  (78) 

As  we  discuss  below,  the  carrier-field  resonance  is  key  to  the  issues  of  stability,  bandwidth 
enhancement  and  noise  reduction.  Because  field  noise  generally  dominates  carrier  noise 
where  there  is  resonance  enhancement,  we  do  not  include  the  carrier-noise  source  terms. 

Many  characteristics  of  a  semiconductor  laser  can  be  modified  by  injection  locking.  Under 
weak  injection  locking,  single-mode  oscillation,  frequency-chirp  reduction  [32],  and  linewidth 
improvement  [33]  in  a  laser  can  be  realized.  Under  strong  injection,  the  dynamic  character¬ 
istics  of  a  laser  can  be  substantially  changed.  Our  recent  analysis  [11,  12]  has  revealed  that 
the  modulation  characteristics  of  a  semiconductor  laser  can  be  significantly  improved  with 
strong  injection  locking  under  proper  conditions.  Using  a  tunable  optical  signal  to  probe 
a  strongly  injection-locked  semiconductor  laser,  we  have  observed  a  parasitic-free  intrinsic 
bandwidth  that  is  beyond  the  -factor  limit  and  is  more  than  three  times  the  bandwidth  of 
the  laser  in  the  free- running  condition  [34].  Besides  bandwidth  enhancement,  our  previous 
analysis  also  indicated  that  strong  injection  locking  can  lead  to  a  broadband  reduction  of 
the  laser  noise  as  well  [11].  However,  an  enhancement  in  the  modulation  bandwidth  is  not 
necessarily  accompanied  by  a  reduction  in  the  laser  noise. 

Referring  back  to  the  bandwidth  enhanced  transceiver  in  Figure  1,  the  specific  application 
will  determine  what  properties  of  the  semiconductor  laser  under  optical  injection  need  to  be 
emphasized.  To  capitalize  on  the  potential  benefit  of  strong  injection  locking  of  a  semicon¬ 
ductor  laser  for  its  practical  application  as  a  large-bandwidth,  low-noise  device,  it  is  clearly 
necessary  that  conditions  exist  for  bandwidth  enhancement,  noise  reduction,  and  stable  lock¬ 
ing  to  be  simultaneously  achieved.  Here,  we  quantify  such  conditions.  It  is  well  known  that 
except  for  very  weak  injection,  stable  locking  occurs  predominantly  with  the  frequency  of  the 
injection  signal  negatively  detuned  from  the  free-running  oscillation  frequency  of  the  laser 
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[18,  19].  However,  even  in  the  stable  locking  region,  a  locking-unlocking  bistability  has  been 
observed  above  a  certain  injection  level  [35].  We  use  the  conventional,  edge-emitting  laser  as 
our  model  system.  Because  its  parameters  have  been  accurately  determined,  the  simulation 
results  agree  exceedingly  well  with  any  experimental  data  that  we  have  collected.  In  the 
cases  where  the  linearized  analytical  model  [12]  applies,  the  analytical  results  are  also  shown 
for  comparison. 

We  fix  the  current  injection  level  at  J  =  2/3,  which  corresponds  to  an  injection  current  level 
of  40  mA  and  a  free-running  output  power  of  9  mW.  At  this  current  level,  the  corresponding 
laser  parameters  that  are  used  in  our  numerical  and  analytical  calculations  for  the  laser 
characteristics  are  7c  =  2.4  x  10^^  s“^,  7,  =  1.458  x  10^  s  7n  =  1.34  x  10®  s  7p  “ 
2.41  X  10®  s“^,  6  =  4,  and  Rgp  =  4.7  x  10^®  V^m~^s“^  Meanwhile,  at  j  =  2/3,  the  resonance 
frequency  of  this  laser  in  the  free-running  condition  is  fr  =  2.93  GHz.  All  of  these  parameters 
were  determined  experimentally  [13]. 

At  a  given  injection  current  level,  the  characteristics  of  a  semiconductor  laser  subject  to  an 
external  optical  injection  depend  on  two  key  parameters  of  the  injection  signal:  the  frequency 
detuning,  H/27r,  of  the  injection  signal  with  respect  to  the  free-running  frequency  of  the  laser 
and  the  injection  parameter,  that  indicates  the  strength  of  the  injection  signal  received 
by  the  laser.  Figure  10  shows  the  mapping  of  different  characteristic  operation  regions  of 
the  laser  under  injection  cis  a  function  of  these  two  parameters.  In  this  figure,  as  well  as  in 
the  two  following  figures,  the  thick  solid  curves  are  obtained  from  the  linearized  analytical 
solution  whereas  the  open  symbols  represent  data  obtained  from  numerical  simulation  using 
the  full  nonlinear  model.  The  thick  broken  curves  are  merely  guides  of  the  eye. 

The  locking  region  shown  in  Figure  10  is  asymmetric  with  respect  to  the  zero  frequency 
detuning  and  opens  up  significantly  at  high  injection  levels.  Curve  a  is  the  Hopf  bifurcation 
boundary  above  which  the  laser  becomes  unstable  though  it  may  be  locked.  Curve  e  is  the 
boundary  below  which  the  laser  cannot  be  locked.  The  region  between  curves  c  and  e  is 
the  locking-unlocking  bistability  region.  Note  that  the  locking-unlocking  bistability  does  not 
occur  at  an  injection  level  lower  than  ^  =  0.06  or  at  a  detuning  less  negative  than  —8  GHz. 
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Figure  10:  Mapping  of  different  characteristic  operation  regions  of  the  laser  under  injection 
as  a  function  of  the  injection  parameter  and  the  frequency  detuning  of  the  injection  signal. 
The  thick  solid  curves  are  obtained  from  the  linearized  analytical  solution  while  the  open 
symbols  are  data  obtained  from  numerical  simulation  using  the  full  nonlinear  model.  The 
thick  broken  curves  are  merely  guides  of  the  eye.  Curve  a  is  the  Hopf  bifurcation  boundary 
below  which  the  laser  is  stably  locked.  Noise  reduction  under  injection  locking  is  found  in 
the  region  below  curve  b.  A  locking-unlocking  bistability  is  observed  in  the  region  bounded 
by  curves  c  and  e.  Bandwidth  enhancement  under  injection  locking  occurs  in  the  region 
above  curve  d.  The  laser  cannot  be  locked  in  the  region  below  curve  e. 
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In  this  bistability  region,  the  laser  can  be  either  locked  or  unlocked,  depending  on  the  initial 
condition.  If  a  laser  is  initially  injection-locked  with  a  set  of  parameters  in  the  stable  locking 
region  between  curves  a  and  c,  it  remains  locked  while  the  injection  signal  is  slowly  tuned  in 
frequency  or  in  strength  to  cross  the  boundary  marked  by  curve  c  into  the  bistability  region. 
Conversely,  if  it  is  initially  unlocked,  it  remains  unlocked  while  the  injection  signal  is  tuned 
within  the  bistability  region  until  curve  c  is  crossed  to  enter  the  stable  locking  region  when 
the  laser  suddenly  becomes  locked.  Therefore,  unconditionally  stable  locking  occurs  only  in 
the  region  between  curves  a  and  c. 

2.3.1  Resonance  Frequency  Enhancement  and  Noise  Reduction 


Returning  to  the  analysis  of  the  SDL  laser.  Figure  11  shows  the  resonance  frequency  of  the 
laser  under  injection  locking  as  a  function  of  the  injection  parameter.  Line  d  marks  the 
resonance  frequency  fr  =  2.93  GHz  for  the  laser  in  the  free-running  condition.  Each  of 
the  thick  curves  corresponds  directly  to  the  curve  that  has  the  same  label  in  Figure  10.  For 
example,  curve  a  in  Figure  11  shows  how  the  resonance  frequency  of  the  laser  under  injection 
locking  varies  along  the  Hopf  bifurcation  boundary  and  curve  c  shows  how  it  varies  along  the 
bistability  boundary.  The  thin  dotted  curves  in  Figure  11  are  contour  curves  for  fixed  values 
of  frequency  detuning.  They  show  how  the  resonance  frequency  varies  with  the  injection 
parameter  when  the  frequency  detuning  is  fixed  at  a  certain  value  in  gigahertz  that  labels  a 
curve. 

While  not  a  precise  correlation,  the  resonance  frequency  is  often  a  good  indicator  of  the  mod¬ 
ulation  bandwidth.  From  Figure  11,  we  find  that  bandwidth  enhancement  due  to  injection 
locking  occurs  within  almost  the  entire  locking  region  of  interest  shown  in  Figure  10.  Only 
in  a  very  narrow  range  between  curves  d  and  e  within  the  bistability  region  does  injection 
locking  not  lead  to  bandwidth  enhancement.  Most  significant  is  the  fact  that  bandwidth 
enhancement  does  occur  in  the  entire  stable  locking  region  between  curves  a  and  c.  In  fact, 
it  can  be  seen  from  Figure  11  that  substantial  bandwidth  enhancement  can  be  realized  any¬ 
where  in  the  stable  locking  region  for  >  0.1.  Along  the  lower  boundary  of  the  locking  region 
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Figure  11;  Resonance  frequency  of  the  laser  under  injection  locking  as  a  function  of  the 
injection  parameter.  Line  d  marks  the  resonance  frequency  fr  =  2.93  GHz  for  the  laser 
in  the  free-running  condition.  Each  of  the  thick  curves  corresponds  directly  to  the  curve 
that  has  the  same  label  in  Figure  10.  The  open  symbols  are  data  obtained  from  numerical 
simulation  using  the  full  nonlinear  model.  The  closed  triangles  are  experimental  data  taken 
under  the  conditions  on  curve  e  of  Figure  4.  The  thin  dotted  curves  are  contour  curves  for 
fixed  values  of  frequency  detuning  in  gigahertz. 
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marked  by  curve  e,  the  experimentally  observed  resonance  frequencies,  which  are  confirmed 
by  the  simulation  data,  disagree  significantly  with  the  analytically  calculated  values. 

The  data  shown  in  Figure  11  seem  to  indicate  that  the  largest  benefit  of  strong  injection 
locking  can  be  realized  by  operating  the  injection-locked  laser  along  or  near  the  Hopf  bifurca¬ 
tion  boundary  as  the  largest  bandwidth  enhancement  occurs  there.  However,  this  conclusion 
does  not  consider  the  noise  of  the  laser.  Figure  12  shows  the  total  power  of  the  amplitude 
noise  of  the  laser  under  injection  locking,  normalized  to  that  of  the  laser  in  the  free-running 
condition,  as  a  function  of  the  injection  parameter.  Line  b  marks  the  noise  level  of  the  laser 
in  the  free-running  condition.  Again,  each  thick  curve  corresponds  directly  to  the  curves 
in  Figures  10  and  11  that  have  the  same  label.  Curves  d  and  e  cannot  be  distinguished  in 
Figure  12  because  they  virtually  overlap.  Each  of  the  thin  dotted  contour  curves  in  Fig¬ 
ure  12  shows  how  the  total  noise  of  the  laser  varies  with  the  injection  parameter  when  the 
frequency  detuning  of  the  injection  signal  is  fixed  at  the  value  in  gigahertz  that  labels  a 
particular  curve. 

It  can  be  seen  that  noise  reduction  by  injection  locking  occurs  only  in  the  region  between 
curves  b  and  e.  In  the  region  above  curve  6,  the  laser  noise  increases  dramatically  as  the 
operation  condition  approaches  the  Hopf  bifurcation  boundary  marked  by  curve  a.  This  is 
expected  because  dynamic  instability  of  the  injected  laser  starts  to  appear  in  the  region  above 
the  Hopf  bifurcation  boundary.  The  significant  effect  of  dynamic  instability  on  the  noise  of 
the  laser  can  be  very  clearly  seen  by  following  the  contour  line  of  zero  frequency  detuning. 
For  injection  at  the  free-running  frequency  of  the  laser  with  zero  frequency  detuning,  the 
laser  undergoes  several  different  nonlinear  dynamic  states,  including  period  doubling  and 
chaos,  when  the  injection  parameter  varies  between  |  =  0  and  (  «  0.18.  The  high  noise 
levels  represented  by  the  zero-detuning  contour  curve  within  this  range  are  a  clear  indication 
of  the  instability  associated  with  such  nonlinear  dynamics. 

Though  the  region  without  noise  reduction  covers  a  large  range  between  curves  a  and  b 
in  Figure  12,  it  represents  only  a  small  range  of  operation  conditions  shown  in  Figure  10 
and  a  still  smaller  range  of  variation  in  the  resonance  frequency  shown  in  Figure  11.  In 
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Figure  12:  Total  noise  power  of  the  Iciser  under  injection  locking,  normalized  to  that  of  the 
laser  in  the  free-running  condition,  as  a  function  of  the  injection  parameter.  Line  b  marks 
the  noise  level  of  the  laser  in  the  free-running  condition.  Each  thick  curve  corresponds 
directly  to  the  curves  in  Figures  10  and  11  that  have  the  same  label.  Curves  d  and  e  cannot 
be  distinguished  in  this  figure  because  they  virtually  overlap.  The  thin  dotted  curves  are 
contour  curves  for  fixed  values  of  frequency  detuning  in  gigahertz. 
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contrast,  the  bistability  region  bounded  by  curves  c  and  e  occupies  a  relatively  large  range 
of  operation  conditions  shown  in  Figure  10  and  an  even  larger  range  of  variation  in  the 
resonance  frequency  shown  in  Figure  11,  but  it  covers  only  a  small  range  of  maximum  noise 
reduction  shown  in  Figure  12.  In  the  large  region  bounded  by  curves  b  and  c,  bandwidth 
enhancement,  noise  reduction,  and  stable  locking  can  be  simultaneously  fulfilled. 

By  correlating  the  data  shown  in  Figures  11  and  12,  we  find  that  when  the  injection  parameter 
is  fixed  at  a  certain  value,  the  frequency  detuning  for  maximum  bandwidth  enhancement 
differs  from  that  for  maximum  noise  reduction.  However,  at  a  given  value  of  frequency 
detuning,  both  bandwidth  enhancement  and  noise  reduction  are  improved  as  the  injection 
parameter  is  increased  over  the  range  of  parameters  considered.  At  injection  levels  above 
those  considered  here,  the  resonance  frequency  can  be  sufficiently  enhanced  that  there  is  a 
significant  dip  in  the  modulation  response  at  frequencies  below  the  resonance  frequency  [12]. 
Significant  bandwidth  enhancement  coupled  with  noise  reduction  under  a  stable  locking 
condition  can  be  realized  with  a  strong  injection  signal  over  a  properly  chosen  frequency 
detuning  range.  For  example,  for  the  laser  biased  at  J  =  2/3,  stable  locking  of  the  laser 
with  an  injection  signal  of  ^  =  0.3  at  a  frequency  detuning  of  —10  GHz  triples  the  resonance 
frequency  to  9  GHz  while  simultaneously  reducing  the  total  laser  noise  by  5  dB.  Greater 
bandwidth  enhancement  with  equivalent  noise  reduction  requires  that  the  laser  be  subjected 
to  a  stronger  injection  field  and/or  be  biased  to  a  larger  value  of  J. 

2.3.2  Dynamic  Range 

In  addition  to  the  bandwidth  and  noise  characteristics,  one  is  also  interested  in  the  dynamic 
range  of  a  semiconductor  laser.  Because  semiconductor  lasers  are  a  nonlinear  medium,  a 
large  modulation  signal  will  eventually  lead  to  saturation  and  to  the  generation  of  additional 
frequencies  in  the  modulation  spectum.  Generation  of  additional  frequencies  also  occurs 
because  of  chirp,  the  tendency  of  the  laser  to  be  both  frequency  and  amplitude  modulated 
by  a  change  in  the  carrier  density.  Previous  work  has  shown  that  operation  under  strong 
injection  locking  leads  to  a  reduction  of  the  frequency  chirp  [36]  and  recent  calculations  using 
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a  linearized  model  similar  to  the  one  described  above  has  shown  that  chirp  reduction  and 
modulation  distortion  can  both  be  achieved  when  an  injection-locked  laser  is  operated  in  the 
region  of  bandwidth  enhancement  [37,  38]. 

As  the  current  program  was  ending,  we  have  started  some  initial  investigations  of  dynamic 
range  and  nonlinear  distortion  to  strong  signal  modulation.  Our  results  are  quite  preliminary. 
As  an  example.  Figure  13  compares  the  relative  gain  for  a  modulation  current  which  is  added 
to  the  steady  (dc)  bias  current.  The  operating  point  discussed  above,  J  =  2/3,  ^  =  0.2  and 
an  offset  frequency  of  -10  GHz,  is  used.  When  the  modulation  index,  m,  the  ratio  of  the 
current  modulation  amplitude  to  the  constant  bias  current  level,  is  approximately  0.1  or 
less,  then  the  current-to-light  gain  closely  follows  the  small  signal  curve.  Compared  to  the 
free-running  case,  the  magnitude  of  the  modulation  is  generally  greater  at  the  modualtion 
frequencies  of  interest.  However,  because  the  injected  optical  signal  simultaneously  increases 
the  average  optical  output,  the  fractional  modulation  is  actually  reduced  relative  to  the  free- 
running  case.  At  higher  modulation  levels,  there  is  some  deviation  from  the  small-signal 
modulation  characteristic,  as  expected.  For  instance,  there  is  a  clear  drop  in  the  gain  for 
modulation  frequencies  near  the  resonance  frequency  for  m  =  1,  as  shown  by  comparing  the 
two  sets  of  calculated  data  in  Figure  13.  Even  here,  however,  the  gain  is  reduced  by  <  3  dB. 
Our  full  nonlinear  calculations  confirm  the  chirp  reduction  previously  calculated  using  the 
linearized  model. 

High  modulation  levels  also  generate  extra  frequency  components  in  the  nonlinear  medium. 
Also  shown  in  Figure  13  is  a  calculated  spectrum  showing  the  amplitude  modulation  signals 
and  noise  spectrum  that  would  be  measured  when  the  laser  is  modulated  by  two  components, 
at  4  GHz  and  7  GHz,  each  with  m  =  0.06.  The  amplitude  of  the  noise  assumes  that  the 
frequency  resolution  for  the  measurement  is  approximately  120  MHz,  and  the  noise  spectrum 
reflects  the  dominant  role  of  the  spontaneous  emission  noise.  Appearing  above  the  noise 
background  are  nonlinear  mixing  peaks  at  3,  8,  11,  and  14  GHz.  These  components  would 
obscure  any  real  signals  at  these  frequencies  with  m  <  0.0006.  This  dynamic  range  compares 
favorably  with  a  free-running  laser  diode.  The  lower  limit  of  the  dynamic  range  is  ultimately 
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Figure  13:  Saturation  and  nonlinear  frequency  generation  due  to  a  large  modulation  current 
which  is  added  to  the  constant  bias  current  of  an  injection-locked  laser.  The  parameters  of  the 
injection-locked  laser  are  discussed  in  the  text.  The  relative  amplification  at  the  modulation 
frequency  is  nearly  constant  for  a  modulation  index  of  m  =  0.1,  triangles,  or  less,  but 
saturates  near  the  resonance  frequency  for  m  =  1,  diamonds.  The  amplitude  modulation 
spectrum,  crosses,  for  the  laser  undergoing  simultaneous  modulation  at  4  GHz  and  7  GHz, 
each  with  a  modulation  index  of  m  =  0.06  shows  harmonic,  sum  and  difference  fequencies 
just  visible  above  the  noise  background.  The  noise  spectrum,  dominated  by  spontaneous 
emission  noise,  assumes  a  frequency  resolution  of  ~  120  MHz. 
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set  by  the  noise  level.  The  background  noise  level  in  Figure  13,  set  by  the  bandwidth  of 
our  calculation,  is  quite  high.  With  a  smaller  frequency  resolution  bandwith,  the  dynamic 
range  can  be  improved  when  the  modulation  index  of  the  peak  frequency  component  is 
reduced.  Further  work  is  underway  to  quantify  modulation  distortion  and  dynamic  range 
characteristics. 
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3.0  EXPERIMENTAL  APPARATUS 


In  this  program,  we  have  carefully  investigated  the  modulation  and  noise  characteristics  and 
nonlinear  dynamics  of  two  very  different  semiconductor  laser  types.  The  first  is  the  5300/5400 
series  of  laser  diodes  from  SDL  when  subject  to  external  optical  injection.  For  our  purposes 
they  are  attractive  because  of  their  nearly  single-mode  operation,  though  as  we  describe 
later,  mode  hopping  sometimes  occured  under  external  optical  injection,  and  the  fact  that  in 
past  experimental  work  we  have  determined  the  key  dynamic  characteristics  of  these  lasers 
[13].  The  second  laser  type  is  a  VCSEL  whose  properties  have  been  published  elsewhere 
[39,  40].  Due  to  its  very  short  cavity,  wavelength-specific  mirrors  and  gain  bandwidth,  the 
VCSEL  oscillates  in  only  one  longitudinal  mode.  As  we  shall  describe  later  in  this  report,  this 
feature  weis  useful  for  experimentally  verifying  specific  features  of  the  nonlinear  dynamics  in 
semiconductor  lasers  under  optical  injection. 

Experimental  observation  of  the  opical  output  of  a  semiconductor  laser  under  external  op¬ 
tical  injection  is  a  two-laser  experiment  while  optical  probing  of  the  laser  under  external 
optical  injection  requires  three  lasers.  The  two-laser  configuration  is  also  useful  for  mea¬ 
suring  the  noise  characteristics  and  dynamic  parameters  of  the  free-running  laser  and  for 
determining  the  characteristic  resonance  frequencies  of  the  injection-locked  laser.  This  is 
because  the  spontaneous  emission  noise  source  effectively  acts  as  a  weak  optical  probe  of 
the  injection-locked  laser.  Therefore,  we  will  first  describe  the  two-laser  configuration  used 
for  the  nonlinear  dynamics  measurements  and  then  continue  with  a  description  of  the  three- 
laser  apparatus  used  for  the  measurement  of  the  optical  modulation  characteristics  of  an 
injection-locked  semiconductor  laser. 

Figure  14(a)  shows  the  two-laser  experimental  setup.  The  output  of  one  laser  diode  (LDl), 
the  master  laser,  at  a  frequency  i/j  waspassed  through  an  optical  isolator  and  variable  at¬ 
tenuators  and  directed  into  a  second  laser  (LD2),  the  slave  laser,  oscillating  at  j/q  =  u;o/27r. 
Thus  the  frequency  offset  /  =  n/27r  is  given  hy  uo  +  f  =  I'l  and  can  be  varied  by  varying 
t/i.  Both  lasers  were  temperature  and  current  stabilized.  To  bring  them  to  near-degenerate 
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Figure  14(b).  Three-layer  configuration. 


_ LEGEND _ 

TLD:  Tunable  laser  diode  A:  Variable  attenuator 

LD:  Laser  diode  BS:  Beamspiitter 

L:  Lens  D:  Fast  photodiode 

‘  Optical  isolator  M.S.A:  R/F*niicrowave  spectrum  analyzer 

A.O.:  Acoustooptic  modulator  O.S.A.:  Optical  spectrum  analyzer 
M:  Mirror 


Figure  14:  Schematic  diagram  of  the  experimental  apparatus,  (a)  two-laser  configuration 
and,  (b)  three-laser  configuration. 
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operating  frequencies,  the  temperature  was  adjusted  for  gross  frequency  changes  and  then 
the  current  was  adjusted  for  fine  tuning  of  the  frequency  offset  f.  Optical  isolators  were 
used  to  avoid  mutual  injection  and  to  reject  back-reflected  light  from  any  components  in 
the  optical  path.  The  output  of  the  injected  laser  was  monitored  using  both  optical  and 
microwave/radio-frequency  spectrum  analyzers  to  measure  both  the  optical  field  and  power 
spectra.  In  this  two-laser  configuration,  the  optical  spectrum  was  monitored  with  a  Newport 
SR-240C  scanning  Fabry-Perot  which  has  a  free  spectral  range  of  2  THz  and  a  finesse  of 
approximately  50,000,  giving  an  optical  frequency  resolution  of  about  40  MHz.  The  power 
spectrum  was  obtained  by  detecting  the  output  of  the  injected  laser  with  a  fast  photodiode 
of  7  GHz  frequency  response  and  displaying  the  photodiode  signal  on  the  spectrum  analyzer. 
Both  spectra  are  necessary  for  the  determination  of  key  dynamic  parameters  required  for  the 
comparison  of  experimental  data  with  numerical  data  [13].  Measuring  both  spectra  simulta¬ 
neously  allowed  us  to  distinguish  between  amplitude  and  phase  modulations  of  the  output 
of  the  injected  laser. 

The  first  master /slave  combination  in  the  two-laser  configuration  consisted  of  two  Spectra 
Diode  Labs  SDL-5301-G1  single-longitudinal  and  transverse  mode  GaAs/AlGaAs  quantum- 
well  lasers  operating  at  827.6  nm.  These  lasers  are  conventional  Fabry-Perot  edge-emitting 
lasers  with  one  facet  coated  for  high  reflection  and  the  output  facet  coated  for  a  reflection 
of  a  few  percent.  The  threshold  current  for  the  injected  laser  was  approximately  24  mA 
and  the  data  presented  here  were  taken  at  an  injection  current  level  of  40  mA,  J  =  2/3. 
Other  dynamic  properties  for  this  laser  are  summarized  in  Table  1  on  page  11.  The  master 
laser  was  operated  at  a  higher  current  level,  approximately  55-65  mA,  to  have  lower  noise 
characteristics.  The  intrinsic  spontaneous  emission  noise  of  the  laser  oscillator  acts  like  a 
weak  optical  probe  of  the  injection-locked  slave  laser. 

During  this  program,  after  several  years  of  consistently  good  performance,  the  slave  laser 
underwent  a  serious  performance  drop  and  could  no  longer  be  overlapped  in  frequency  with 
the  master  laser.  Overlapping  two  of  the  SDL  lasers  is  not  a  trivial  matter  because  they  do 
not  have  a  continuous  tuning  spectrum  under  temperature  and  current  changes.  The  lasers 
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undergo  longitudinal  mode  hops  at  points  which  vary  randomly  from  laser  to  laser.  We  tried 
a  variety  of  laser  combinations  using  new  SDL-5411-G1  lasers.  However,  while  the  original 
lasers  were  made  with  a  500  /xm  cavity,  the  newer  lasers  all  had  a  750  ^m  cavity.  The  longer 
cavity  meant  a  reduced  value  of  the  laser-cavity  photon  decay  rate.  This  translates  into  a 
smaller  modulation  bandwidth  for  both  the  free-running  and  injection-locked  lasers,  all  other 
parameters  being  equal.  The  longer  cavities  did  allow  the  lasers  to  be  operated  at  higher 
pump  currents,  yielding  larger  output  powers.  After  some  trial  and  error,  we  found  that  the 
best  combination  of  available  lasers  was  to  use  the  original  master  laser  as  the  slave  lasers 
and  one  of  the  newer,  higher  power  lasers  as  the  pump  laser.  We  were  able  to  achieve  good 
frequency  overlap  at  approximately  826.8  nm. 

With  an  operating  master  and  slave  laser,  we  could  proceed  with  the  modifications  re¬ 
quired  for  measuring  optical  modulation  characteristics  of  the  injection-locked  slave  laser. 
A  schematic  of  the  new  experimental  configuration  is  shown  in  Figure  14b.  An  optically 
isolated  master  laser,  LDl,  operating  at  the  optical  frequency  i>i,  injection  locks,  in  a  dy¬ 
namically  stable  state,  a  slave  laser,  LD2,  initially  operating  at  the  free-running  frequency 
i/Q.  The  output  from  LD2,  after  passing  through  other  optical  isolators,  is  detected  by  either 
a  scanning  Fabry-Perot  optical  spectrum  analyzer  or  a  fast  photodiode,  microwave  spec¬ 
trum  analyzer  combination.  The  output  of  an  optically  isolated,  tunable,  external  cavity 
semiconductor  laser,  TLD,  operating  at  the  frequency  z/i,  is  passed  through  an  acoustooptic 
modulator  where  part  of  the  beam  is  split  off  and  frequency  shifted  by  80  MHz.  This  shifted 
output  is  used  to  probe  the  modulation  characteristics  of  the  slave  laser  while  the  unshifted 
part  of  the  beam  is  used  as  a  local  oscillator.  When  the  microwave  spectrum  analyzer  is 
set  to  a  narrow  frequency  window  at  80  MHz,  then,  as  the  probe  laser  is  scanned  about 
!/£,,  the  microwave  spectrum  analyzer  measures  the  spectrum  of  the  regeneratively  amplified 
probe  signal.  When  the  local  oscillator  beam  is  blocked  so  that  only  the  slave  laser  output  is 
captured  by  the  fast  photodiode,  then  the  microwave  spectrum  analyzer  measures  the  am¬ 
plitude  modulation  spectrum  of  the  laser.  An  upgraded  photodiode  extended  the  frequency 
response  of  the  amplitude  modulation  measurements  to  15  GHz. 
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The  slave  has  very  similar  dynamic  parameters  to  the  same  type  of  laser  that  we  characterized 
earlier  [13].  Specifically,  7c  =  2.5  x  10“  s"',  7,  =  1.5  x  10^  s-^  7n  =  2Jx  10^  s'^  and 
7p  =  3.6  jx  10®  s“\  where  J  is  the  differential  injection  current  above  threshold,  normalized 
to  the  threshold  current.  A  different  value  was  found  for  the  linewidth  enhancement  factor, 
5  =  5  ±  1.  For  the  measurements  reported  here,  J  =  2.13.  The  master  laser  was  an  SDL 
5411-Gl  laser.  Both  lasers  were  temperature  and  current  stabilized.  The  probe  laser,  a  New 
Focus  Model  6226  external  cavity  semiconductor  laser,  could  be  scanned  over  a  calibrated 
frequency  range  of  approximately  40  GHz  by  dithering  a  mirror  mounted  on  a  piezoelectric 
crystal. 

Feedback  to  the  slave  laser  from  its  collimating  lens  was  discovered  to  play  a  role  in  the 
understanding  of  details  of  some  of  the  regenerative  amplification  and  amplitude  modulation 
spectra  that  we  measured.  The  effect  of  feedback  into  the  oscillating  mode  of  the  injection- 
locked  slave  laser  spectrum  will  be  discussed  later.  Here,  we  will  briefly  discuss  the  effects  of 
feedback  on  the  optical  probe.  Because  the  probe  laser  follows  the  same  optical  train  as  the 
master  laser,  its  coupling  into  the  slave  laser  could  not  be  simultaneously  optimized  with  the 
master  laser.  Even  when  feedback  effects  on  the  master  laser  and  injection-locked  slave  laser 
were  minimized,  it  was  found  that  the  probe  laser  was  influenced  by  etaloning  effects  due  to 
the  collimating  lens.  An  example  is  shown  in  Figure  15.  The  figure  shows  the  raw  data  of 
two  regenerative  amplification  spectra  and  the  noise  level  of  the  measurement  system.  One 
of  the  spectra,  with  the  strong  peak,  shows  the  measurement  of  the  free-running  laser  biased 
at  J  =  2.13.  The  frequency  axis  of  the  figure  is  referenced  to  the  free- running  oscillation 
frequency.  The  other  spectrum  is  taken  at  the  same  optical  frequency  but  with  the  bias 
reduced  to  J  =  0.5.  The  oscillation  frequency  has  now  shifted  so  that  the  laser  resonance 
is  not  recovered  in  the  heterodyne  spectrum.  However,  there  are  clearly  repeating  maxima 
and  minima  offset  by  approximately  10  GHz.  Inspection  of  the  J  =  2.13  spectra  reveals  that 
the  same  periodic  feature  is  evident.  The  frequency  spacing  identifies  a  cavity  formed  by 
the  slave  laser  and  the  collimating  lens.  The  fact  that  the  feature  is  approximately  the  same 
in  strength  in  both  the  on-resonance  and  off-resonance  spectra  shows  that  the  etaloning  is 
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Frequency  Offset  (GHz) 


Figure  15:  Spectra  of  the  regenerative  amplification  sideband  due  to  a  weak  optical  probe 
showing  the  effects  of  a  weak  optical  cavity  for  the  probe  laser  formed  by  the  slave  laser  and 
its  collimating  lens.  The  cavity  does  not  use  the  oscillating  mode  of  the  slave  laser  as  is  seen 
by  the  nearly  equal  strength  etaloning  in  a  resonant  spectrum  taken  at  J  =  2.13,  solid  curve, 
and  a  nonresonant  spectrum  taken  at  J  =  0.5,  dashed  curve.  .Also  shown  is  the  noise  level 
of  the  measurement  system,  dotted  curve. 
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not  due  to  coupling  into  the  oscillating  mode.  The  weak  optical  probe  has  a  frequency- 
dependent  amplitude  due  to  the  frequency-dependent  feedback  of  the  weak  optical  cavity 
formed  for  the  probe  beam.  The  effects  of  this  cavity  can  be  taken  into  account  by  simply 
correcting  for  the  amplitude  variation.  Because  the  spectrum  analyzer  gives  the  signal  in 
dB,  one  simply  cancels  out  the  etalon  effects  by  adding  an  equal  etalon  component  to  the 
spectrum  which  is  180  degrees  out  of  phase.  All  subsequent  modulation  data  presented  here 
have  been  corrected  for  this  linear  effect  on  the  raw  data. 

Data  from  one  other  master-slave  combination  are  also  presented  in  this  report.  A  New 
Focus  Model  6226  tunable  external  cavity  laser  was  the  master  laser  and  a  VCSEL  was  the 
injected  laser.  Unfortunately,  from  the  perspective  of  getting  useful,  quantitative  results, 
the  VCSEL  had  the  habit  of  spontaneously  changing  its  output  characteristics,  such  as 
threshold  current,  output  wavelength,  and  dynamic  parameters,  during  the  course  of  this 
research  program.  The  source  of  this  feature  is  not  clear.  Such  behavior  has  not  been 
reported  in  the  literature,  but,  in  private  conversations,  other  experimentalists  investigating 
VCSELs  have  observed  such  behavior.  Throughout  much  of  the  program,  and  for  most  of  the 
data  presented  in  this  report,  the  VCSEL  had  the  characteristics  described  below.  When 
the  VCSEL  changed,  it  could  be  recharacterized  by  our  model  with  appropriate  dynamic 
parameters.  Initially,  the  VCSEL  had  a  threshold  for  laser  oscillation  of  approximately  3.9 
mA  and  emission  wavelength  of  approximately  850  nm  [39, 40, 41].  The  total  output  power  as 
a  function  of  bias  current  is  shown  in  Figure  16.  At  currents  under  approximately  6  mA,  the 
VCSEL  oscillates  in  the  lowest  order  spatial  mode  but,  due  to  the  symmetry  of  the  device,  can 
emit  in  two  orthogonally  polarized  modes.  Also  shown  in  Figure  16  is  the  polarized  output 
power  as  a  function  of  the  injection  current.  The  higher  threshold  polarization  component 
eventually  dominates  the  output  and,  when  the  output  starts  to  be  divided  among  several 
spatial  modes,  all  show  this  polarization.  The  laser  operated  with  these  characteristics  during 
the  period  that  the  threshold  data  and  the  polarization  switching  data  that  are  presented 
later  were  taken.  The  dynamic  characteristics  of  the  laser  were  reported  in  Reference  [41]. 
More  recently,  the  laser  operated  with  a  threshold  of  approximately  3.2  mA.  Among  the 
dynamic  parameters  the  most  notable  change  was  the  linewidth  enhancement  factor  which 
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Figure  16:  Output  power  of  the  VCSEL  under  study  at  bias  currents  where  the  output  power 
is  in  one  spatial  mode  with  two  orthogonal  polarization  modes. 
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Wcis  reduced  from  approximately  5-6  to  2-3.  The  data  describing  operation  above  threshold 
and  showing  the  locking/unlocking  bistability  were  taken  with  these  characteristics. 
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4.0  SEMICONDUCTOR  LASER  ABOVE  THRESHOLD 


When  a  semiconductor  laser  is  subjected  to  injection  by  a  near-resonant  external  optical 
field,  the  laser  output  will  exhibit  a  variety  of  optical  spectra.  By  comparison  between  the 
experimentally  observed  spectra  and  calculated  spectra,  we  have  been  able  to  categorize  the 
dynamics.  Calculated  spectra  are  obtained  by  numerically  integrating  the  normalized,  cou¬ 
pled  equations,  then  performing  a  Fourier  transform  on  the  resulting  time  series  to  calculate 
the  spectra. 

As  we  have  already  discussed,  some  key  features  can  be  predicted  from  the  perturbation 
analysis  with  the  set  of  coupled  equations  which  assumes  that  the  changes  to  the  amplitude, 
phase  and  carrier  density  are  sufficiently  small  that  the  three  equations  can  be  linearized. 
Figure  17  shows  a  typical  mapping  of  the  predicted  locking  and  stability  range  based  on  a 
linearized  analysis  [12].  The  actual  values  used  to  calculate  the  map  correspond  to  the  in¬ 
jected  SDL  laser  at  the  bias  current  of  40  mA,  J  =  2/3,  but  the  qualitative  characteristics  are 
more  general.  The  analysis  predicts  an  asymmetric  locking  range  [19]  with  a  bounded  range 
of  unstable  operation  within  the  locking  range.  The  circulating  field  amplitude  increases 
with  increasing  injection  parameter  and/or  decreasing  offset  frequency  within  the  locked 
region.  A  larger  value  of  b  increases  the  low  frequency  boundary  of  the  locked  operating 
range  and  increases  the  range  of  unstable  dynamics.  A  larger  value  of  J  shifts  the  unstable 
region  to  larger  values  of  ^  with  the  low  ^  limit  shifting  approximately  in  proportion  to  J 
or  Aq  and  the  upper  limit  shifting  approximately  with  Aq  [12].  The  analysis  also  predicts 
that  the  eigenvalues  of  the  coupled  equations  are  strong  functions  of  the  operating  point. 
Figure  18  plots  the  eigenvalues,  or  resonance  frequencies  of  the  field-free  carrier  coupling,  as 
a  function  of  ^  for  two  different  cases.  The  first  case  is  for  external  optical  injection  at  the 
free-running  frequency  of  the  injected  laser.  This  corresponds  to  a  phase  offset  of  arctan  b 
between  the  master  and  injected  laser  optical  fields  in  the  linearized  analysis.  For  this  case 
the  resonance  frequencies  increase  monotonically  with  In  contrast,  when  the  injection 
frequency  is  changed  as  ^  is  increased  so  that  there  is  no  phase  offset  between  the  master 
and  injected  laser  optical  fields,  a  case  which  approximately  corresponds  to  following  the  low 
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Figure  17:  Mapping  of  the  locking  and  stability  characteristics  of  the  conventional,  edge- 
emitting  semiconductor  laser  with  J  =  2/3  as  a  function  of  the  frequency  offset  between  the 
locking  and  free- running  laser  fields  and  the  injection  parameter,  The  thick  solid  line  is 
the  boundary  of  the  injection-locked  region.  Shown  in  the  interior  of  the  locking  region  is 
the  Hopf  bifurcation  line  (thick  dashed  line)  separating  the  stable  and  unstable  operating 
regimes.  The  positive  frequency  boundary  of  the  locked  region  is  a  line  of  constant  output 
field  amplitude  at  the  free-running  level,  Aq.  Also  shown  are  the  lines  for  1.1  Aq  (dotted) 
and  1.2.4o  (dot-dashed). 
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Figure  18:  Calculated  eigenfrequency  and  measured  frequency  of  the  optical  spectrum  side¬ 
band  peak  for  the  conventional,  edge-emitting  semiconductor  laser  with  J  =  2/3  as  a  func¬ 
tion  of  the  injection  parameter,  Experimental  data  are  shown  as  solid  symbols,  calculated 
frequencies  using  the  full  nonlinear  equations  are  shown  as  open  symbols  and  calculated 
eigenfrequencies  from  the  linearized  approximation  are  shown  as  curves.  Two  cases  are  plot¬ 
ted:  Dashed  Line  and  Squares  -  injection  at  the  free-running  optical  frequency.  Dotted  Line 
and  Triangles  -  injection  corresponding  to  =  0.  The  former  case  shows  unstable  operation 
for  0.007  <  ^  <  0.24  while  the  latter  case  is  always  stable. 
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frequency  boundary  of  the  locked  operating  region,  the  linearized  analysis  predicts  a  decreas¬ 
ing  resonance  frequency.  Also  shown  in  Figure  18  are  the  resonance  frequencies  calculated 
by  numerically  integrating  the  full  nonlinear  equations  and  Fourier  transforming  the  result¬ 
ing  time  series  and  the  experimentally  measured  resonance  frequencies.  The  experimental 
data  show  excellent  agreement  with  the  full  nonlinear  calculation,  and  good  agreement  with 
the  linearized  analysis.  The  major  discrepancies  with  the  linearized  analysis  occurs  within 
the  region  where  the  full  nonlinear  analysis  indicates  chaotic  dynamics  at  the  <f>i  =  arctan  b 
phase  offset  or  bistable  operation  at  the  phase  offset  of  =  0. 

4.1  NONLINEAR  DYNAMICS 


Using  equations  74  to  76  with  experimentally  determined  parameters,  we  have  been  able  to 
produce  calculated  spectra  quite  similar  to  the  observed  spectra.  When  the  noise  source  is 
eliminated  from  the  calculated  spectra,  the  spectra  can  be  directly  related  to  the  underlying 
dynamics  induced  by  the  optical  injection  [7,  8,  9].  As  an  example.  Figure  19  shows  the 
progression  of  observed  spectra  from  the  conventional,  edge-emitting  laser  when  subjected 
to  resonant,  fl  =  0,  optical  injection  as  the  injection  parameter,  is  increased.  To  determine 
the  injection  level,  we  calibrated  the  system  by  measuring  the  four-wave  mixing  spectrum  in 
the  weak  injection  limit  and  used  our  previously  developed  techniques  to  compare  the  gener¬ 
ated  sideband  signal  with  the  central  peak  [13].  The  injected  power  was  varied  by  changing 
the  attenuation  in  the  path  of  the  master  laser  and  the  power  level  was  monitored.  In  the 
spectra,  the  frequency  is  relative  to  the  free-running  frequency  or,  equivalently,  the  injection 
frequency  and  the  scale  of  the  signal  level  is  consistent  throughout  the  series.  Expanded 
versions  of  the  spectra  showing  weak  details  are  referenced  to  the  right  vertical  axis  scale. 
Figure  19(a)  is  the  free-running  spectrum  consisting  of  a  single  peak  and  weak  relaxation 
resonance  sidebands  which  are  visible  in  the  expanded  scale  version  of  the  spectrum.  The 
asymmetry  is  due  to  6.  In  Figure  19(b)  =  1  x  10"^,  and  the  spectrum  consists  of  relatively 

narrow  peaks  separated  by  a  frequency  spacing  of  2.9  GHz,  which  is  the  relaxation  resonance 
frequency  of  the  free-running  laser.  These  are  typical  features  of  highly  unstable  injection 
locking  [19].  In  Figure  19(c)  =  3  x  10“^,  broad  period-doubling  features  appear  in  the 
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Figure  19:  Measured  optical  spectra  of  the  conventional,  edge-emitting  semiconductor  laser 
with  J  =  2/3  under  optical  injection  at  the  free- running  optical  frequency:  (a)  free- running 
operation  with  expanded  scale  showing  the  relaxation  resonance  sidebands,  (b)  unstable 
injection  locking,  limit  cycle,  at  =  1  x  10“'*  with  expanded  scale  spectrum,  (c)  period 
doubling  at  =  3  x  10““*  with  features  shown  clearly  in  the  expanded  scale,  (d)  chaotic 
dynamics  at  =  1.3  x  10“^,  (e)  limit  cycle  at  =  3.6  x  10“^,  (f)  period  doubling  at 
=  7.2  X  10“^,  (g)  transition  back  to  limit  cycle  at  =  1.7  x  10“^,  and  (h)  stable 
operation  at  =  8.5  x  10"^  along  with  expanded  scale  spectrum  and  the  expanded  scale 
sideband  spectrum  of  the  meister  Iciser. 
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spectrum  between  the  narrow  oscillation  peaks  [8].  They  are  stronger  at  negative  detunings, 
a  reflection  of  the  positive  value  of  b.  Comparing  the  expanded  versions  of  (b)  and  (c)  clearly 
shows  the  qualitative  diffference  between  the  spectra  before  and  after  the  period-doubling 
features  emerge.  The  field  noise  source  causes  the  broadening  of  the  period-doubling  fea¬ 
tures.  When  it  is  removed  in  the  calculated  spectra,  sharp  period-doubling  peaks  result  [8]. 
In  Figure  19(d),  =  1.3  x  10“^,  the  spectrum  becomes  dominated  by  a  broad  pedestal  and 

many  secondary  peaks  develop.  Relatively  little  energy  remains  in  the  narrow  injection  spike. 
At  this  stage  chaos  has  fully  developed  and  elimination  of  the  noise  source  in  the  calculated 
spectra  does  not  change  the  qualitative  features  [8].  Within  the  region  of  chaotic  dynamics, 
a  fraction  of  the  oscillating  power,  increasing  to  up  to  35%  as  the  injection  level  is  increased, 
is  shifted  from  the  principal  oscillating  mode  into  several  of  the  weak  side  modes.  Over 
a  narrow  injection  range,  the  broadened  spectrum  collapses  again  into  an  equally  spaced 
set  of  narrow  features,  as  shown  in  Figure  19(e)  where  =  3.6  x  10“^,  and  the  principal 
mode  regains  its  full  power  [7].  Now,  however,  the  separation  of  the  peaks  has  increased 
and  the  spectrum  is  much  more  strongly  shifted  to  the  negative  components.  At  still  higher 
injection  levels,  Figure  19(f)  where  =  7.2  x  10“^,  a  new,  clear  period  doubling  is  observed 
with  a  further  increase  of  the  peak  separation  and  relative  strengthening  of  the  negative 
frequency  components.  The  strongest  spectral  feature  is  clearly  frequency-shifted,  pushed 
by  the  resonant  injection.  The  period  doubling  peaks  then  steadily  decrease  in  amplitude 
as  the  resonance  peak  monotonically  increases  while  the  peak  separation  continues  to  in¬ 
crease.  The  spectra  in  Figure  19(g)  is  observed  when  =  1.7  x  10"^  and  the  injection  peak 
again  begins  to  dominate  the  spectrum.  Eventually,  the  side  peaks  decrease  in  magnitude  so 
that  only  a  single  peak  at  the  injection  frequency  dominates  the  spectrum.  In  Figure  19(h) 
=  3,5  X  10“^  and  stable  dynamics  has  been  reestablished.  The  peak  separation  has  now 
increased  to  10.5  GHz  and  the  side  peak  is  visible  in  the  expanded  spectrum.  Also  visible 
in  the  expanded  spectrum  are  replicas  of  the  relaxation  resonance  features  of  the  master 
laser.  The  identiflcation  of  these  features  is  confirmed  by  comparison  with  sideband  spec¬ 
trum  of  the  master  laser,  also  shown.  At  high  injection  levels,  the  noise  sidebands  of  the 
master  laser  appear  in  the  spectrum  of  the  injection-locked  slave.  These  noise  features  do 
not  reflect  the  dynamics  of  the  injected  laser  system.  The  dynamics  are  reflected  in  the 
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shifted  resonance  feature.  Note  that  the  strong  asymmetry  of  the  resonance  peaks  continues 
when  stable  dynamics  has  been  reestablished.  The  weak  asymmetry  in  the  resonance  peaks 
of  the  free-running  laser  has  become  a  strong  asymmetry  in  the  injection-locked  laser.  The 
measured  shift  in  the  resonance  peaks  is  compared  with  the  calculated  eigenfrequencies  of 
the  coupled  equations  in  Figure  18  and  shows  very  good  agreement  with  the  full  nonlinear 
calculation  and  reasonable  agreement  with  the  linearized  treatment. 

To  understand  the  deterministic  dynamics,  the  full  nonlinear  coupled  equations  are  first 
solved  with  the  noise  source  terms  set  to  zero.  Figure  20  depicts  the  numerically  obtained 
bifurcation  diagram  of  the  values  of  the  extrema  of  the  amplitude,  a{t),  versus  the  injection 
parameter,  As  the  injection  level  is  increased,  the  steady  state  is  destabilized  and  the 
relaxation  frequency  is  undamped  (Hopf  bifurcation),  in  agreement  with  the  linear  stabil¬ 
ity  analysis.  Further  increasing  the  injection  leads  to  a  period-doubling  bifurcation  route  to 
chaos  and  then  a  similar,  but  reversed,  route  out  of  chaos.  To  confirm  the  chaotic  nature,  the 
Lyapounov  exponents  and  the  corresponding  Kaplan- Yorke  dimension  of  the  attractor  have 
been  calculated  [42].  For  ^  =  0.03,  the  three  Lyapounov  exponents  are  0.01457c,  0.0,  and 
-0.02987c,  and  the  Kaplan- Yorke  dimension  is  2.48.  The  positive  Lyapounov  exponent  is  a 
clear  measure  that  the  system  is  chaotic  in  this  parameter  range.  At  even  higher  injection 
levels,  the  bifurcation  diagram  shows  a  second  period-doubling.  Finally,  at  the  highest  injec¬ 
tion  levels,  the  laser  diode  reverts  to  stable  operation  (reverse  Hopf  bifurcation).  Therefore, 
the  deterministic  dynamics  is  seen  to  largely  determine  the  experimentally  observed  optical 
optical  spectra  and  to  reproduce  all  of  the  key  changes  in  the  spectral  characteristics  which 
we  observe. 

When  both  coherent  injection  and  spontaneous  emission  are  present,  there  is,  effectively,  a 
fluctuating  injected  field.  This  leads  to  a  blurring  of  the  period-doubling  cascades  into  and 
out  of  chaos.  In  the  forward  cascade,  only  the  first  period-doubling  is  not  obscured,  and 
the  new  frequency  components  are  severely  broadened.  The  reverse  cascade  is  completely 
obscured.  This  effect  can  be  numerically  recovered  by  including  spontaneous  emission  noise 
sources  in  the  amplitude  and  phase  equations.  Figure  21  compares  calculated  spectra  with 
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Figure  20:  Numerically  calculated  bifurcation  diagram  of  the  extrema  of  the  normalized 
optical  field  amplitude,  a{t),  versus  the  normalized  injection  level,  Injection  is  at  the 

free-running  frequency  of  the  slave  laser.  Experimentally  determined  parameters  of  the 
semiconductor  laser  are  used  in  the  calculation. 


and  without  the  noise  source  terms  at  three  injection  levels.  Without  noise,  a  region  of  period 
doubling,  Figure  21(a),  and  period  quadrupling,  Figure  21(b),  are  clearly  distinguishable  in 
the  forward  cascade.  With  noise,  however,  the  calculations  show  only  broadened  period¬ 
doubling  features.  Figures  21(d)  and  (e),  which  are  similar  to  the  experimental  data  shown 
in  Figure  19(c).  In  the  period-doubling  bubble  at  higher  injection  levels,  however,  the  noise 
does  not  severely  broaden  the  new  spectral  peaks,  as  shown  in  Figures  21(c)  and  (f).  Likewise, 
the  measured  spectrum  in  Figure  19(f)  has  sharp  period-doubling  features. 

Before  proceeding  to  a  more  general  picture  of  the  dynamics  induced  under  external  optical 
injection,  we  distinguish  a  special  type  of  period-doubling-like  spectra,  the  region  of  sub¬ 
harmonic  resonance  [10].  Here,  the  half-offset  frequency  feature  can  totally  dominate  the 
frequency  component  oflFset  from  the  central  peak  by  the  resonance  frequency.  Figure  22 
shows  a  set  of  spectra  which  illustrate  the  subharmonic  resonance.  Again,  the  offset  fre¬ 
quency  is  referenced  to  the  frequency  of  the  free-running  laser.  In  these  spectra  the  master 
laser  is  offset  by  «  5.5  GHz  and  the  injection  level  is  varied.  A  complementary  set  of  spectra, 
where  the  injection  level  is  held  constant  while  the  offset  frequency  is  varied  from  below  to 
above  twice  the  relaxation  resonance  frequency  of  the  free-running  laser,  has  been  presented 
elsewhere  [10].  At  low  injection  levels  the  spectrum  consists  of  narrow,  weak  sidebands, 
equally  and  oppositely  offset  from  a  main  peak  with  its  relaxation  resonance  sidebands,  as 
shown  in  Figure  22(a).  These  narrow  features  are  due  to  regenerative  amplification  of  the 
injected  signal  and  four- wave  mixing  between  the  injected  signal  and  the  central  peak  [13]. 
As  the  injected  signal  is  increased,  the  spectral  features  near  the  relaxation  resonance  fre¬ 
quency  narrow  and  increase.  Figure  22(b)  where  the  injection  level  is  =  2.1  x  10“^.  Note, 
however,  that  these  features  are  halfway  between  the  injection  feature,  or  its  oppositely 
shifted  partner,  and  the  central  peak.  This  frequency  is  slightly  less  than  the  relaxation 
resonance  frequency  of  the  free-running  laser.  There  is  also  a  small  frequency  pushing  of 
the  central  peak.  At  still  higher  injection  levels.  Figure  22(c)  with  =  4  x  10"^,  the  sub¬ 
harmonic  resonance  nature  of  the  dynamics  becomes  clear.  The  sideband  at  the  injection 
frequency  is  now  dominated  by  the  feature  at  half  the  offset.  Also,  the  frequency  pushing 
has  stopped  and  the  central  peak  is  slightly  pulled  toward  the  injected  frequency  [10].  As 
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Figure  21:  Computed  optical  spectra  showing  the  effects  of  the  spontaneous-emission  noise 
at  three  levels  of  injection,  (a)  and  (d)  ^  =  0.0145,  (b)  and  (e)  (  =  0.0152.  (c)  and  (f) 
^  =  0.064.  (d)-(f)  include  the  spontaneous-emission  noise  source  term  while  the  others  do 
not.  The  strongest  spectral  feature  in  (a),  (b),  and  (d),  and  the  two  strongest  in  (c)  and  (f) 
are  clipped  to  emphasize  the  detailed  features. 
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Figure  22:  Measured  optical  spectra  of  the  conventional,  edge-emitting  semiconductor  laser 
with  J  =  2/3  under  optical  injection  at  a  frequency  offset  «  5.5  GHz:  (a)  weak  injection 
with  regenerative  amplification  and  four-wave  mixing  sidebands,  (b)  enhanced  subharmonic 
signal  at  =  2.1  X  10“^,  (c)  subharmonic  resonance  at  =  4  x  10““*.  (d)  weak  period 
quadrupling  at  =  2.6  x  10~^,  (e)  period  doubling  at  =  6  x  10“^,  and  (f)  limit  cycle  at 
^2  =  1.2  X  10-2. 
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the  injection  level  is  increased  to  =  2.6  x  10~^,  Figure  22(d),  the  pushing  of  the  central 
peak  and  the  dominance  of  the  resonance  frequency,  now  significantly  larger  than  the  initial 
offset  freqency,  are  reestablished.  This  is  actually  a  period-quadrupling  spectrum,  with  weak, 
broadened  features  just  discernible  above  the  noise  level  in  between  the  some  of  the  sharp 
spectral  features.  At  higher  injection  levels  the  spectrum  follows  the  progression  to  period 
doubling.  Figure  22(e)  with  =  6  x  10“^,  and  to  limit  cycle  oscillations,  Figure  22(f)  with 
=  1.2  X  10“^.  The  frequency  pushing  continues  to  increase  and  the  pushed  feature  gets 
relatively  weaker.  At  very  high  injection  levels,  >  3  -  4  x  10“^,  the  laser  output  begins 
to  smoothly  shift  to  other  longitudinal  modes. 

Using  spectra  like  that  of  Figures  19  and  22,  we  have  constructed  a  mapping  of  the  dynamics 
for  the  laser  at  this  pump  level  as  the  injection  parameter  and  the  frequency  offset  are 
varied.  In  the  mapping  of  Figure  23,  the  frequency  axis  is  relative  to  the  free-running 
frequency  of  the  injected  laser.  The  map  symbols  are:  4  -  perturbation  spectrum  with  weak 
regenerative  amplification  and  four-wave  mixing  sidebands;  S  -  stable  injection  locking;  SR 
-  subharmonic  resonance;  PI  -  limit  cycle  oscillation;  P2  -  period  doubling;  P4  -  period 
quadrupling;  Chaos  -  deterministic  chaos;  M’  -  multiwave  mixing  with  most  output  on 
another  longitudinal  mode;  hatched  regions  -  principal  output  on  another  longitudinal  mode; 
thin  lines  -  smooth  transition  between  dynamic  regions;  thick  dotted  lines  -  abrupt  mode 
hop  transitions  with  minor  hysteresis;  thick  dashed  lines  with  arrow  -  one  way  mode  hops 
out  of  mode;  thick  solid  lines  -  abrupt  transition  to/from  a  region  of  chaos  or  multiwave 
mixing  where  there  is  significant  power  in  another  longitudinal  mode,  from/to  a  region  with 
power  primarily  in  the  principal  mode.  The  smooth  transitions  represented  by  the  thin  lines 
are  approximate.  For  instance,  a  peak-to-sideband  ratio  of  10:1  was  used  for  the  transition 
line  from  stable  to  unstable  dynamics.  For  small  values  of  the  optical  injection  acts  as  a 
perturbation.  Weak  sidebands  are  generated  at  the  offset  frequency  of  the  external  optical 
injection  and  equally  and  oppositely  shifted  from  the  central  peak.  These  are  the  regenerative 
amplification  and  four-wave  mixing  sidebands,  respectively.  As  the  injection  parameter  is 
increased,  various  dynamic  instabilities  develop.  At  negative  detunings  the  locked  region 
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Figure  23:  Mapping  of  the  experimentally  observed  dynamic  regions  for  a  conventional  edge- 
emitting  semiconductor  laser  with  J  =  2/3.  The  dynamics  are  determined  by  comparison 
of  observed  spectra  with  calculated  spectra  with  and  without  noise.  Symbols  are  defined  in 
the  text. 
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resembles  the  predictions  of  the  mapping  from  the  linearized  analysis.  There  is  a  narrow 
region  of  stable  injection  locking  at  lower  values  of  (  which  opens  up  for  ^  ^  0.1. 

Most  of  the  dynamic  regions  and  transitions  are  directly  recovered  in  the  single-mode  model 
of  semiconductor  laser  operation.  The  value  of  the  linewidth  enhancement  factor,  6,  is  a 
critical  parameter  in  determining  the  limits  of  the  range  of  nonlinear  dynamics  and  whether 
that  range  will  include  regions  of  chaotic  dynamics.  From  the  mapping  it  can  be  seen  that 
the  second  region  of  period-doubling  dynamics  along  the  fi  =  0  line  is  associated  with 
a  second  region  of  chaotic  dynamics  at  positive  offset  frequencies.  The  positive  detuning 
region  is  more  complicated  because  boundaries  between  different  types  of  nonlinear  dynamics 
are  difficult  to  draw.  For  instance,  when  the  offset  is  approximately  equal  to  the  free- 
running  relaxation  resonance  frequency,  the  distinction  between  the  perturbation,  four-wave- 
mixing-type,  spectrum  and  the  limit-cycle  spectrum  blurs  into  one  of  qualitative  difference. 
There  are  no  clear  features  in  the  spectra  which  distinguish  between  locked  and  unlocked 
operation.  The  various  unstable  dynamics  show  no  qualitative  distinction  as  they  cross 
the  line  where  the  linearized  analysis  divides  locked  and  unlocked  operation.  Within  the 
nominally  locked  operating  regime,  the  various  spectral  features  of  the  optical  spectrum 
clearly  show  broadening  with  respect  to  the  narrow  feature  at  the  injection  frequency  and 
are,  therefore,  not  locked  to  the  injection  source. 

Also  associated  with  injection  at  positive  offset  frequencies  is  the  general  trend  of  frequency 
pushing  of  the  originally  free-running  spectral  feature.  One  anomally  of  the  subharmonic 
resonance  is  the  slight  pulling  of  the  original  oscillation  peak  in  that  region.  For  optical 
injection  at  positive  offset  frequencies,  the  excitation  has  a  tendency  to  push  the  original 
oscillation  peak  to  negative  offset  frequencies  with  the  degree  of  pushing  increasing  as  the 
offset  frequency  of  the  excitation  is  decreased  or  the  injection  parameter  is  increased.  Figure 
24  plots  the  frequency  pushing  for  the  injection  offset  at  5.5  GHz.  More  details  of  the  weak 
frequency-pulling  effect  around  the  subharmonic  resonance  are  given  elsewhere  [10].  Here, 
we  concentrate  on  the  general  frequency-pushing  trend.  The  pushing  effect  highlights  a  key 
point  which  is  central  to  the  asymmetry  of  the  mapping.  The  external  optical  injection 
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Figure  24:  Frequency  pushing  of  the  initially  free-running  spectral  feature  as  a  function  of 
injection  with  the  master  laser  frequency  offset  by  5.5  GHz,  experimental  data  -  squares  and 
dotted  line,  full  nonlinear  equations  -  solid  line,  eigenvalues  of  linearized  equations  -  dashed 
line.  Frequency  pulling  in  the  subharmonic  resonance  region  causes  the  abrupt  dip  in  the 
experimental  data  and  nonlinear  calculation  curves. 
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shifts  the  Fabry-Perot  oscillation  frequency  of  the  laser  cavity  by  changing  the  steady-state 
carrier  density.  The  excitation  follows  the  new  Fabry-Perot  frequency  along  the  zero  phase 
offset  line  in  the  stable  operating  region.  When  the  optical  excitation  is  off  of  the  Fabry- 
Perot  resonance,  there  is  a  new  frequency  in  the  coupled  field-carrier  system.  The  shifting  of 
the  resonance  oscillation  frequency  with  excitation  level  is  due  to  the  competition  between 
the  original  relaxation  resonance,  and  its  associated  damping,  and  the  frequency  difference 
between  the  injected  signal  and  this  new  resonance,  and  its  associated  phase  offset,  to  control 
the  response.  The  system  restabilizes  within  the  injection-locked  operating  region  when  the 
Fabry-Perot  frequency  shifts  and  associated  damping  are  large  compared  to  the  original 
system  resonances  [12].  In  this  region  of  stable  operation  under  strong  injection,  where 
linearized  dynamics  again  prevail,  the  modulation  characteristics  of  the  semiconductor  laser 
are  strongly  modified. 

4.2  ENHANCED  MODULATION  BANDWIDTH  AND  REDUCED  NOISE  IN 
INJECTION-LOCKED  SEMICONDUCTOR  LASERS 


Recently,  several  studies  using  conventional  coupled-equation  models  have  predicted  that  the 
modulation  bandwidth  of  strongly  injection-locked  semiconductor  lasers  can  be  significantly 
improved  relative  to  the  free-running  case  [11,  12,  43,  37].  However,  experimental  work  to 
date  has  only  confirmed  that  the  resonance  frequencies  of  an  injection-locked  laser  scale  in 
accordance  with  the  predictions  of  the  model  [12].  In  a  free- running  semiconductor  laser,  the 
maximum  current  modulation  bandwidth  is  set  by  the  /i-factor  [44],  K  =  (27r)^(7„-b7p)/flr, 
and  fmax  =  Here,  7„  and  7p  are  the  decay  rates  due  to  stimulated  emission  and  gain 

saturation,  respectively,  and  flr/2jr  is  the  relaxation  resonance  frequency,  =  7c7„  -f  75 7p, 
where  7c  is  the  photon  decay  rate  from  the  laser  cavity  and  7,  is  the  free-carrier  spontaneous 
decay  rate.  Because  7„  and  7p  are  proportional  to  the  circulating  power  in  the  laser  cavity, 
fmax  gives  the  limiting  modulation  bandwidth  at  high  power  levels.  Here,  we  demonstrate 
that  at  accessible  operating  conditions  an  injection-locked  laser  can  reach  higher  modulation 
bandwidths  than  fmax  for  the  laser. 
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Semiconductor  lasers  have  an  asymmetric  frequency  range  for  injection-locked  operation 
[19].  This  is  because  a  change  in  the  free  carrier  density  modifies  both  the  gain  and  the 
refractive  index.  The  change  in  the  refractive  index  is  related  to  the  change  in  the  gain  by 
the  linewidth  enhancement  factor,  b.  When  the  injected  optical  power  is  within  a  bounded 
range,  typically  from  10"^  to  lO'^  of  the  output  power,  a  sizeable  fraction  of  the  locked 
operating  region  displays  unstable  dynamics.  However,  for  strong  injection  locking,  at  a 
level  higher  than  10”^,  stable  dynamics  again  dominates  the  locked  operating  region  and 
bandwidth  enhancement  has  been  calculated  [12].  In  regions  of  stable  dynamics,  a  small- 
signal  analysis  is  expected  to  be  appropriate  for  calculating  the  modulation  characteristics 
due  to  a  weak  optical  probe  or  current  modulation.  A  weak  optical  probe  at  a  frequency 
offset  from  the  locked  operating  frequency  generates  a  regenerative  amplification  sideband  at 
the  probe  frequency.  When  there  is  no  optical  feedback  into  the  slave  laser,  the  small-signal 
regenerative  amplification  spectrum  is  given  by  [12]: 

+  (K  +  h,L/2)fi  +  ny2  +  (7,  +  T'nL)2£//2p 

|A.f  1L>|2 

-|-[— n(7rL  -  7pI-/2)  +  V{l3  +  Inl)  +  b^\l2f].  ("9) 

Here,  Ar  and  A,-  are  the  amplitudes  of  the  regenerative  amplification  sideband  and  the 
injected  probe,  respectively,  t]  is  the  coupling  parameter  for  external  fields  into  the  laser, 
and  H/2ir  is  the  offset  between  the  probe  frequency  and  the  operating  frequency  injection- 
locked  laser.  Injection  locking  shifts  the  operating  point  of  the  complex  gain  as  described 
by  7„L  =  7nai,  7pL  =  7poi,  U  =  (7/7c)aL  cos(^l)  and  V  =  {v/lc)aLsm{(t>L),  where  = 
]Ax,|/|Ao|  is  the  ratio  of  the  locked  to  free-running  field  amplitudes  and  is  the  phase  offset 
between  the  circulating  and  injected  fields.  Also,  =  (7c  ~  2C/)7nL  +  '^slpL 

D  =  {-in  +  U)[-Q‘^  +  0.1  +  {7s  +  InLW  -  i^7rL] 

+{nl  -  io^pLW  -f  {-in  -h  7.  +  7nL)v\  (80) 

where  7rL  =  Is  +  lnL  +  lpL  +  U.  Note  that  we  write  2f//2  for  one  term  in  Equation  79  because 
the  factor  of  2  in  the  denominator  was  incorrectly  included  in  Equation  47  of  Reference  [12]. 
Similarly,  the  amplitude  modulation  spectrum  due  to  a  weak  current  modulation  that  would 
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be  measured  by  direct  detection  of  the  injection-locked  laser  output  is  described  by  [12,  37, 
43]: 


\^^  +  iU  +  bVnhc-2U)WnL  (81) 

\(^\  «  pp 

where  d  is  the  complex  amplitude  modulation  coefficient.  Fitting  the  regenerative  ampli¬ 
fication  spectrum  determines  all  of  the  parameters  needed  to  calculate  the  parasitic-free 
amplitude  modulation  characteristic  due  to  a  weak  current  modulation. 

The  strength  of  the  signal  from  the  master  laser  is  given  by  the  injection  parameter,  (  = 
{‘ql'fc)aL,  which  can  be  experimentally  determined  [8].  The  injected  power  is  proportional 
to  Figure  25  shows  the  regenerative  amplification  spectra  for  the  free- running  laser  and 
for  three  different  operating  points  when  the  laser  is  injection  locked  at  the  highest  injection 
level  we  could  achieve,  ^  =  0.4  i  0.1.  The  frequency  axis  is  with  respect  to  the  free- running 
frequency  of  the  slave  laser.  The  signal  axis  is  relative  to  the  system  noise  level.  Changing 
the  offset  frequency  of  the  master  laser  with  respect  to  the  free-running  slave  laser  modifies 
the  regenerative  amplification  characteristic  of  the  injection- locked  laser.  Detunings  of  the 
master  laser  of  0,  -4.8  and  -9.7  GHz  are  shown.  Each  of  the  injection-locked  laser  spectra 
contains  a  narrow  peak  at  the  injection  frequency  corresponding  to  the  injection-locked  peak. 
It  is  not  part  of  the  true  modulation  spectrum.  Along  with  the  experimental  data,  calculated 
spectra  are  shown  using  the  experimentally  determined  parameters.  For  the  calculations, 
the  overall  signal  level  is  a  free  parameter.  It  was  adjusted  to  match  the  free-running  data. 
The  injection  parameter  was  fixed  at  ^  =  0.4.  The  value  of  the  linewidth  enhancement 
factor  was  then  adjusted,  within  its  allowable  range,  to  give  a  good  fit  for  the  frequency 
offset  of  the  resonance  peak  in  the  zero  detuning  data  yielding  b  =  5.4.  The  amplitude  of 
the  signal  is  no  longer  a  free  parameter  for  these  data.  There  are  then  no  free  parameters 
for  the  fitting  of  the  detuned  injection  spectra.  Note  that  the  injection-locked  laser  spectra 
are  strongly  asymmetric,  reflecting  the  shift  of  the  Fabry-Perot  cavity  resonance  frequency 
under  the  optical  injection  [12]. 
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Figure  25:  Spectra  of  the  regenerative  amplification  sideband  due  to  a  weak  optical  probe 
under  four  different  operating  conditions.  The  frequency  axis  is  relative  to  the  free-running 
frequency  of  the  slave  laser  and  the  signal  axis  is  relative  to  the  noise  level  of  the  measurement 
system.  Solid  symbols  are  experimental  data  and  curves  are  calculated  spectra  using  the 
experimentally  determined  dynamic  parameters  and  Equation  52.  Solid  Line  and  Squares 
-  Free-running  laser.  Dashed  Line  and  Diamonds  -  Injection  locking  at  the  free-running 
frequency,  Dotted  Line  and  Triangles  -  Injection  locking  offset  at  -4.8  GHz,  Dot-dashed 
Line  and  Ovals  -  Injection  locking  offset  at  -9.7  GHz. 
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The  good  agreement  between  the  measured  and  calculated  spectra  allows  us  to  calculate 
the  amplitude  modulation  spectra  due  to  a  small  oscillating  current  when  the  laser  is  at 
the  operating  point  corresponding  to  the  four  regenerative  amplification  spectra.  The  cal¬ 
culated  spectra,  normalized  to  the  zero-frequency  signal  level,  are  plotted  in  Figure  26.  The 
free-running  laser  shows  the  characteristic  relaxation  resonance  peak  in  the  spectrum.  The 
modulation  spectra  of  the  injection-locked  laser  all  show  a  slight  dip  before  a  shifted  reso¬ 
nance  peak.  Based  on  the  deviation  between  the  calculated  and  measured  spectra  for  the 
detuned  injection  locking,  the  calculated  curves  slightly  overestimate  the  resonance  frequency 
and  underestimate  the  width  of  the  resonance  peak  of  these  two  spectra.  For  the  laser  un¬ 
der  investigation,  fmax  ~  19.3  GHz,  a  value  higher  than  can  be  practically  achieved  due  to 
constraints  on  the  injection  current  level.  All  three  injection-locked  modulation  curves  have 
3-dB  rolloff  frequencies  higher  than  fmax- 

The  dip  in  the  current  modulation  characteristic  highlights  the  fact  that,  in  the  strong 
injection  limit,  the  spectrum  is  made  up  of  two  resonances  [12].  The  zero-frequency  resonance 
has  a  damping  corresponding  primarily  to  the  field-enhanced  free-carrier  decay  rate,  75+7ni- 
The  circulating  power  increases  as  the  detuning  becomes  more  negative  within  the  injection- 
locked  operating  region,  decreasing  the  size  of  the  dip.  The  non-zero  resonance  corresponds 
to  the  offset  of  the  shifted  Fabry-Perot  resonance  frequency  with  respect  to  the  injection- 
locked  peak.  In  the  large  injection  limit,  the  offset  is  largely  determined  by  V ,  with  associated 
damping  determined  by  U.  As  the  detuning  of  the  master  laser  becomes  more  negative,  4>l 
decreases  towards  zero  in  the  injection-locked  region.  This  explains  the  combined  decrease 
in  the  resonance  frequency  and  increase  in  the  damping  rate. 

In  the  design  of  high-speed  semiconductor  lasers,  great  attention  is  paid  to  decreasing  the 
values  of  b  and  7p  [26].  In  the  injection-locked  case,  this  is  not  as  important.  The  value 
of  b  needs  to  be  kept  low  to  minimize  the  frequency  chirp.  However,  under  strong  injec¬ 
tion  locking,  it  has  recently  been  shown  that  the  frequency  chirp  is  greatly  reduced  in  the 
bandwidth-enhanced  operating  region  [37].  Also,  increasing  b  increases  the  frequency  shift  of 
the  Fabry-Perot  resonance,  all  other  parameters  being  unchanged.  Decreasing  7p  yields  an 
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Figure  26:  Calculated  amplitude  modulation  spectra  due  to  a  weak  current  modulation 
corresponding  to  the  regenerative  amplification  spectra  plotted  in  Figure  25.  All  curves  are 
normalized  to  the  signal  at  zero  frequency  offset.  Solid  Line  -  Free-running,  Dashed  Line 
-  Injection  locking  at  the  free-running  frequency,  Dotted  Line  -  Injection  locking  offset  at 
—4.8  GHz.  Dot-dashed  Line  -  Injection  locking  offset  at  —9.7  GHz. 
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increased  fmax  through  the  /i'-factor.  However,  the  /f-factor  does  not  apply  to  the  injection- 
locked  laser.  Under  injection  locking,  increasing  7p  has  the  principal  effect  of  increasing  the 
damping  rates  of  the  charactersitic  resonances.  Therefore,  it  can  actually  enhance  the  useful 
operating  range  of  the  system. 

We  have  shown  that  parasitic-free  modulation  bandwidths  in  excess  of  20  GHz,  and  in  excess 
of  the  A'-factor  limit,  can  be  achieved  in  relatively  low  cost,  commercially  available  lasers  by 
using  strong  injection  locking.  The  bandwidth  enhancement  relative  to  the  free-running  case 
was  approximately  a  factor  of  three.  Factors  of  2-3  improvement  in  bandwidth  are  calculated 
to  be  achievable  in  a  wide  range  of  laser  diodes,  and  parasitic-free  bandwidths  in  excess  of 
60  GHz  appear  to  be  feasible  with  high-speed  lasers. 

We  have  previously  described  how  the  noise  spectra  of  a  laser  diode  can  be  understood  in 
terms  of  the  external  injection  model  in  the  weak  injection  limit  [15].  This  analysis  was  then 
extended  to  the  case  of  the  injection- locked  laser  to  understand  changes  in  the  amplitude 
noise  spectrum  [11].  Specifically,  data  was  taken  of  the  injection-locked  laser  when  the 
offset  frequency  is  changed  to  maintain  the  <f>L  —  0  condition.  In  the  previous  publication, 
the  injection-locked  data  showed  only  qualitative  agreement  with  calculated  spectra.  The 
discrepancy  was  largely  due  to  improperly  accounting  for  the  increased  circulating  power 
under  injection  locking  in  the  previously  calculated  spectra.  Figure  27  replots  the  previously 
published  data  of  the  noise  amplitude  spectra.  It  compares  the  data  with  properly  scaled 
calculated  spectra  using  the  linearized  analysis  and  full  nonlinear  calculations.  There  is  now 
good  agreement  between  calculated  and  measured  spectra.  As  long  as  the  laser  is  not  biased 
too  close  to  the  threshold  for  laser  oscillation,  J  =  0,  then  the  noise  can  be  considered  a  linear 
perturbation  and  the  noise  improvements  calculated  above  are  recovered  experimentally. 

4.3  FEEDBACK  EFFECTS  ON  MODULATION  SPECTRA  IN 
INJECTION-LOCKED  SEMICONDUCTOR  LASERS 

During  the  course  of  our  experimental  investigations  on  bandwidth  enhancement  in  injection- 
locked  laser  diodes,  we  found  that  we  observed  spectra  which  could  not  be  fully  explained 
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Figure  27:  Chajiges  in  the  amplitude  noise  spectrum  of  a  quantum-well  laser  diode,  where 
the  dominant  noise  source  is  spontaneous  emission,  due  to  stable  injection  locking  by  an 
external  laser.  The  experimental  data,  solid  symbols,  are  for  the  detuning  that  minimized 
the  amplitude  of  the  noise  spectrum.  Two  different  calculated  spectra  are  shown.  The  open 
symbols  are  the  calculated  spectra  using  the  full  nonlinear  coupled  equations  that  minimized 
the  noise  spectrum.  The  curves  are  the  calculated  spectra  at  the  different  injection  levels 
for  the  offset  frequency  where  oi  =  0.  Dotted  Line  and  Diamonds  -  ^  =  0,  Solid  Line  and 
Squares  -  ^  =  .06,  Dashed  Line  and  Triangles  -  ^  =  .2. 
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without  including  feedback  into  the  slave  laser  cavity  from  its  collimating  lens.  We  have 
already  described  how  the  lens  and  laser  surfaces  formed  a  weak  cavity  for  the  linear  probe 
laser.  Here,  we  discuss  the  effects  of  the  feedback  on  our  experimentally  observed  spectra 
using  the  expanded  model  for  laser  diode  operation,  with  the  effects  of  feedback  included, 
that  was  described  above.  These  formulas  are  appropriate  when  the  external  optical  injection 
is  sufficiently  strong  to  induce  stable,  injection-locked  operation  and  the  feedback  is  not  so 
strong  that  it  destabilizes  the  injection-locked  operating  point. 

When  there  is  feedback  into  an  injection-locked  laser,  there  is  a  shift  in  the  phase,  cI>l, 
and  normalized  amplitude,  fl£,,  of  the  output  field  as  described  by  Equations  14-16.  These 
changes  will  be  reflected  in  the  offset  frequency  and  damping  characteristics  of  the  resonance 
associated  with  the  shifted  cavity  Fabry-Perot  resonance  frequency.  Depending  upon  the 
phase  of  the  feedback  signal  due  to  the  round-trip  delay,  the  damping  can  be  either  enhanced 
or  reduced.  The  effects  of  feedback  are  further  reflected  in  the  modulation  characteristics 
with  the  additional  etaloning  associated  with  the  feedback  cavity. 

The  effects  of  feedback  can  be  clearly  seen  in  the  regenerative  amplification  spectra  of  Fig¬ 
ure  28  and  amplitude  modulation  spectra  of  Figure  29.  The  experimental  data  was  taken 
with  the  laser  subject  to  an  injection  level  of  ^  =  0.2  ±  0.05  and  the  master  laser  was  offset 
from  the  free-running  frequency  of  the  slave  leiser  by  —22  GHz.  Two  calculated  spectra, 
with  and  without  feedback,  are  shown  in  each  figure.  In  Figure  28  they  are  calculated  from 
Equation  52  and  in  Figure  29  from  Equation  54.  The  calculated  spectra  used  (  =  0.225  and 
the  strength  and  phase  shift  of  the  feedback  was  adjusted  to  give  a  good  fit  to  the  data. 
The  calculated  spectra  with  feedback  used  r  =  0.055  and  arccos(u;LT)  =  2.945,  nearly  out 
of  phase  with  the  externally  injected  signal.  The  other  parameters  were  the  same  as  for 
the  data  above  which  did  not  show  the  feedback.  The  optical  power  which  is  fed  back  into 
the  cavity  is  approximately  an  order  of  magnitude  less  than  the  injected  signal  and,  from  a 
realistic  applications  point  of  view,  represents  fairly  strong  undesired  feedback. 
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Figure  28:  Spectra  of  the  regenerative  amplification  sideband  due  to  a  weak  optical  probe 
showing  the  effects  of  feedback  on  an  injection-locked  laser.  The  frequency  axis  is  relative  to 
the  injection-locked  operating  frequency  which  is  shifted  from  the  free-running  frequency  of 
the  slave  laser  by  —22  GHz.  The  signal  axis  is  relative  to  the  noise  level  of  the  measurement 
system.  Open  symbols  are  experimental  data  and  curves  are  calculated  spectra  using  the 
experimentally  determined  dynamic  parameters  and  Equation  52.  The  injection  level  is 
^  =  0.225  with,  Dashed  Curve-  feedback  amplitude  is  r  =  0,  and  Dotted  Curve  -  feedback 
amplitude  is  r  =  0.055  and  phase  is  arccos(u;£,T)  =  2.945. 
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Figure  29:  Spectra  of  the  amplitude  modulation  due  to  a  weak  optical  probe  showing  the 
effects  of  feedback  on  an  injection-locked  laser.  The  frequency  axis  is  relative  to  the  injection- 
locked  operating  frequency  which  is  shifted  from  the  free-running  frequency  of  the  slave  laser 
by  —22  GHz.  The  signal  axis  is  relative  to  the  noise  level  of  the  measurement  system.  Open 
symbols  are  experimental  data  and  curves  are  calculated  spectra  using  the  experimentally 
determined  dynamic  parameters  and  Equation  52.  The  injection  level  is  ^  =  0.225  with, 
Dashed  Curve-  feedback  amplitude  is  r  =  0,  and  Dotted  Curve  -  feedback  amplitude  is 
r  =  0.055  and  phase  is  &tccos{u>lt)  =  2.945. 
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As  further  examples  of  the  effects  feedback  can  have  on  the  modulation  spectra,  we  show 
calculated  regenerative  amplification  and  current  modulation  induced  amplitude  modula¬ 
tion  spectra.  The  calculations  compare  spectra  which  use  parameters  relevant  to  the  data 
described  above  with  spectra  where  either  the  amplitude  or  the  phase  of  the  feedback  has 
been  changed.  Figures  30  and  31  show  that  when  the  phase  or  amplitude  of  the  feedback 
are  adjusted,  one  can  observe  both  enhanced  and  reduced  stability  of  the  shifted  resonance 
peak.  The  curve  with  the  phase  shifted  by  tt  corresponds  to  feedback  close  to  being  in 
phase  with  the  externally  injected  signal.  Such  feedback  can  diminish  the  magnitude  of  the 
resonance  peak  rather  than  enhance  it.  When  the  feedback  is  approximately  out  of  phase 
with  the  injected  signal,  as  was  the  case  with  the  experimental  data,  an  increase  in  the 
feedback  enhances  the  magnitude  of  the  resonance  peak. 

To  a  reasonable  approximation,  the  effects  of  feedback  at  this  level  are  the  same  as  changing 
the  strength  and  frequency  offset  of  the  injected  optical  signal.  Indeed,  the  data  could  be 
adequately  fit  using  the  feedback-free  formulas  of  Reference  [12]  with  injection  parameters 
and  frequency  offsets  that  did  not  match  the  experimental  conditions.  If  one  operates  too 
close  to  the  stability/instability  boundary  of  injection- locked  operation,  feedback  can  cause 
a  nominally  stable  system  to  undergo  amplitude  and  phase  oscillations  at  the  resonance 
frequency.  Therefore,  as  long  as  one  operates  well  within  the  unconditionally  stable  operating 
region  and  the  feedback  level  is  not  too  high,  feedback  will  not  preclude  the  bandwidth 
enhancement  and  improved  modulation  characteristics. 

4.4  VCSEL  ABOVE  THRESHOLD 


The  effects  of  optical  injection  on  the  optical  spectrum  of  a  VCSEL  are  shown  in  Figure  32. 
For  these  spectra,  the  VCSEL  was  operated  at  a  bias  current  of  6  mA.  Again,  zero  frequency 
corresponds  to  the  free-running  optical  frequency  of  the  VCSEL.  The  master  laser  was  also 
tuned  to  this  frequency  for  this  series  of  spectra.  At  this  injection  current  level,  there  is  some 
output  power  in  higher  order  spatial  modes  which  are  spectrally  offset  from  the  region  shown. 
Optical  injection  at  the  free-running  frequency  of  the  dominant,  lowest  order  spatial  mode, 
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Figure  30:  Spectra  of  the  regenerative  amplification  sideband  characteristic  due  to  a  weak 
optical  modulation,  showing  the  effects  of  optical  feedback.  The  frequency  offset  is  with 
respect  to  the  free-running  optical  frequency.  The  dotted  and  solid  lines  correspond  to  the 
calculated  curves  of  Figure  28  with  and  without  feedback,  respectively.  The  dashed  curve  is 
for  the  feedback  amplitude  reduced  to  r  =  0.0275  and  the  phase  of  the  feedback  unchanged 
at  arccos(a;i,T)  =  2.945  while  the  dot-dashed  curve  is  for  the  amplitude  of  the  feedback 
unchanged  at  r  =  0.055  and  the  phase  shifted  by  tt  to  arccos(a.’£,T)  =  6.087. 
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Figure  31:  Spectra  of  the  amplitude  modulation  characteristic  due  to  a  weak  current  mod¬ 
ulation.  The  dotted  and  solid  lines  use  the  same  parameters  as  the  calculated  curves  of 
Figures  28  and  29  with  and  without  feedback,  respectively.  The  dashed  curve  is  for  the 
feedback  amplitude  reduced  to  r  =  0.0275  and  the  phase  of  the  feedback  unchanged  at 
arccos(u;£,T)  =  2.945  while  the  dot-dashed  curve  is  for  the  amplitude  of  the  feedback  un¬ 
changed  at  r  =  0.055  and  the  phase  shifted  by  tt  to  arccos(a.‘£,r)  =  6.087. 
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Figure  32:  Measured  optical  spectra  of  the  VCSEL  with  J  =  0.9  under  optical  injection  at 
the  free-running  frequency.  The  frequency  axis  is  referenced  to  the  free-running  frequency 
of  the  VCSEL.  The  injection  conditions  are:  (a)  free-running  and  relative  injection  power 
levels  of,  in  arbitrary  units  (b)  0.045,  (c)  0.234,  (d)  0.441,  (e)  0.907,  and  (f)  1.536  where  the 
laser  has  reestablished  stable  operation. 
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the  case  shown  here,  leads  to  a  decreeise  in  the  power  in  the  higher  order  modes.  However, 
it  did  not  fully  quench  this  emission.  The  set  of  spectra  in  Figure  32  show  an  initially  free- 
running  laser,  32(a),  which  undergoes  a  Hopf  bifurcation  to  a  dynamically  unstable  state  as 
the  injection  level  is  increased,  32(b).  The  resonance  frequency  monotonically  increases  as 
the  injection  level  increases,  32(c)-32(e).  Finally,  the  laser  reestablishes  stable  dynamics  at 
very  high  injection  levels,  32(f).  The  progression  is  similar  to  the  one  previously  observed. 
The  laser  did  not  show  period  doubling  and  a  transition  to  chaotic  dynamics  because  the 
linewidth  enhancement  factor  had  a  relatively  low  value,  6  =  2-3.  Earlier  data  taken  when 
the  laser  displayed  a  larger  linewidth  enhancement  factor  clearly  displayed  the  transition  to 
chaos.  All  of  these  features  are  consistent  with  the  dynamics  expected  from  the  coupled 
equation  model. 

The  increase  in  the  resonance  frequency  as  a  function  of  the  injection  level  is  shown  in 
Figure  33.  Here,  the  resonance  frequency  again  shows  the  nearly  linear  dependence  with 
the  injection  parameter  seen  in  the  earlier  data  on  the  conventional,  edge-emitting  laser. 
We  have  previously  shown  that  this  dependence  falls  to  two-thirds  power,  cube  root  of  the 
injected  power,  when  the  injected  power  becomes  large  compared  to  the  free-running  output 
power  [41].  Both  of  these  results  follow  from  the  coupled  equation  model  [12].  The  spectra 
in  Figures  32(c)-(e)  show  that  there  are  two  dominant  components  of  the  optical  spectrum 
as  the  resonance  frequency  varies  over  a  significant  range  of  microwave  frequencies.  The 
beating  of  these  two  optical  components  leads  to  a  strong  microwave  signal  upon  conventional 
detection  of  the  optical  beam.  This  tunability  is  achieved  by  changing  the  injected  power 
from  a  constant  output  power  optical  beam  without  the  requirement  for  any  high-speed 
electrical  imput.  The  tuning  range  shown  in  Figure  33  should  not  be  taken  as  a  performance 
limit  for  VCSELs  in  general  or  even  the  VCSEL  under  study.  Others  have  observed  the 
changing  resonance  frequency  in  VCSELs  to  be  as  large  as  60  GHz  [45]. 

Referring  back  to  the  mapping  of  the  conventional,  edge-emitting  laser.  Figure  23,  some 
features  of  the  mapping  reflect  specific  characteristics  of  the  laser  under  investigation,  a 
conventional,  edge-emitting,  Fabry-Perot  laser  diode,  where  the  oscillating  mode  is  always 
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Figure  33:  Variation  of  the  resonance  frequency  of  the  VCSEL  under  external  optical  injec¬ 
tion  as  a  function  of  the  square  root  of  the  injection  power.  The  squares  are  the  experimental 
data  and  the  solid  line  is  a  visual  aid. 


in  competition  with  other  longitudinal  modes  which  are  close  to  threshold.  Because  of 
this,  the  external  optical  signal  can  induce  a  mode  hop  and  certain  regions  of  the  detuning- 
injection  plane  become  inaccessible.  There  is  an  abrupt  mode  hop  near  the  locking-unlocking 
boundary  which  has  a  small  hysteresis,  not  shown,  associated  with  it.  For  larger  offsets,  some 
of  the  mode  hops  are  one  way  and  the  laser  does  not  reestablish  operation  on  the  original 
longitudinal  mode  at  this  boundary.  At  negative  offsets  there  is  a  bounded  range  of  unstable 
dynamics  which  follows  a  period-doubling  progression,  largely  obscured  by  noise,  to  chaotic 
dynamics.  Within  the  range  of  chaotic  dynamics  a  varying  fraction  of  the  optical  output 
begins  to  shift  to  other  modes.  One  edge  of  the  region  of  chaotic  dynamics  is  an  abrupt 
transition  where  power  returns  to  the  original  longitudinal  mode  as  the  laser  returns  to 
oscillatory  dynamics.  Again  there  is  a  small  hysteresis,  not  shown,  associated  with  this 
transition. 

Because  the  longitudinal  modes  in  a  VCSEL  have  such  a  large  frequency  spacing,  only  one  is 
resonant  with  the  high-reflectivity  mirrors  that  define  the  laser  cavity.  Therefore,  the  VCSEL 
should  not  show  these  mode  hops  and  we  did  not  observe  them.  In  addition,  it  should  be 
able  to  show  the  locking/unlocking  bistability  that  was  apparently  masked  in  the  other  laser 
by  the  mode  hop.  The  bistability  could  be  observed.  Figure  34  shows  a  set  of  spectra  taken 
with  the  bias  current  to  the  VCSEL  set  at  6  mA  and  the  detuning  of  the  master  laser  varied 
while  the  injection  level  was  held  constant.  When  the  master  laser  frequency  offset  is  shifted 
to  negative  values  away  from  the  free-running  slave  frequency,  the  laser  remains  locked  until 
the  detuning  is  beyond  —12  GHz.  However,  when  the  detuning  shift  is  then  reversed,  locking 
is  not  reestablished  until  —10.8  GHz.  Between  —12  and  —10.8  GHz  the  laser  displays  the 
locking/unlocking  bistability. 

4.5  VCSEL  POLARIZATION  AND  NOISE  CHARACTERISTICS 


The  fact  that  the  VCSEL  has  two  nearly  degenerate,  orthogonally  polarized  modes  can 
complicate  the  interpretation  of  the  laser  spectra.  Even  when  the  output  power  of  one 
polarization  dominates  the  other  by  a  factor  of  100,  the  influence  of  the  weak  polarization 
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Figure  34:  Measured  optical  spectra  of  the  VCSEL  with  J  =  0.9  as  the  injection  frequency 
offset  is  shifted.  The  frequency  axis  is  referenced  to  the  free-running  frequency  of  the  VCSEL. 
Locked  injection  as  the  injection  frequency  is  shifted  away  from  resonance;  (a)  0  GHz  offset, 
(b)  -4.8  GHz.  (c)  -12  GHz.  Unlocked  injection  as  the  injection  frequency  is  shifted  towards 
resonance:  (d)  —15  GHz,  (e)  —12.6  GHz,  and  (f)  — lO.S  GHz. 
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component  is  easily  observed.  First,  it  causes  a  broadening  of  the  laser  linewidth  by  more 
than  a  factor  of  10.  Second,  it  causes  a  resonance  feature  in  the  amplitude  (power)  spectrum. 
Figure  35  shows  a  key  characteristic  of  this  resonance.  In  the  figure  the  resonance  is  centered 
at  approximately  600  MHz,  0  dB  corresponds  to  the  noise  level  of  the  measurement  apparatus 
and  —5  dB  is  the  resolution  limit.  Three  different  power  spectra  are  shown.  When  the 
measured  polarization  component  is  directed  to  be  along  either  of  the  two  orthogonally 
polarized  components,  the  resonance  feature  is  at  a  relative  maximum.  However,  when  the 
measured  polarization  direction  is  rotated  to  be  45  degrees  between  the  two  polarization 
components,  equivalent  to  measuring  the  total  spectrum  without  polarization  selection,  the 
resonance  feature  essentially  disappears.  This  shows  that  the  two  polarization  components 
display  anti-correlated  dynamics.  Power  in  one  component  comes  explicitly  at  the  expense 
of  the  other  so  that  the  total  power  is  left  unchanged. 

External  optical  injection  into  either  of  the  two  components  can  cause  an  enhancement  or  a 
quenching  of  the  resonance  depending  on  the  specific  injection  frequency  and  power.  This  is 
similar  to  the  previously  discussed  case  where  the  optical  injection  can  cause  unstable  dyan- 
mics  or  stable  injection  locking  in  the  single  polarization  case.  Optical  spectra  under  optical 
injection  are  shown  in  Figure  36.  Here,  the  laser  Weis  operated  at  J  =  0.31,  zero  frequency 
corresponds  to  the  free-running  frequency  of  the  dominant  peak  and  the  orthogonal  polar¬ 
ization  component  is  offset  by  approximately  8.5  GHz.  In  both  cases,  injection  on  the  strong 
peak.  Figure  36(a),  and  the  weak  peak,  Figure  36(b),  the  injection  matches  the  polarization 
direction  of  the  injected  component.  In  the  weak  component  spectra,  there  is  residual  signal 
from  the  strong  component  passed  by  our  polarizers.  In  this  spectra,  therefore,  we  see  the 
effect  of  the  injection  on  the  stronger  orthogonal  component.  The  orthogonal  component 
shows  resonance  sidebands  due  to  the  injection.  Little  effect  was  observed  when  the  po¬ 
larization  of  the  injected  signal  was  orthogonal  to  the  injected  component.  The  injection 
strength  was  not  sufficient  to  induce  unstable  dynamics  directly  in  the  strong  component  at 
the  offset  of  the  weak  component.  While  the  central  peak  of  both  polarization  components 
is  narrowed  by  the  injection,  indicating  a  quenching  of  the  resonance,  the  dominant  polar¬ 
ization  component  is  clearly  destabilized.  Even  more  complicated  spectra  due  to  optical 
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Figure  35;  Low-frequency  power  spectrum  showing  the  resonance  due  to  the  antimode  cou¬ 
pling  between  the  two  orthogonally  polarized  modes.  The  VCSEL  was  operated  at  J  =  0.31 
where  the  dominant  polarization  is  approximately  two  orders  of  magnitude  stronger  than 
the  orthogonal  component.  The  resonance  disappears  when  the  output  polarized  45  degrees 
between  the  two  orthogonal  polarizations  is  passed  to  the  photodetector. 
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injection  are  observed  when  the  two  polarization  components  are  of  nearly  equal  strength  in 
the  free-running  case.  Figure  37  shows  the  case  when  J  =  0.38  and  the  measured  polarization 
direction  is  45  degrees  between  the  two  polarization  components.  Injection  on  either  peak 
induces  a  complicated  structure.  The  details  of  this  structure  are  not  well  understood.  How¬ 
ever,  they  clearly  represent  additional  noise-like  background  in  any  modulation  spectrum. 
For  this  reason,  it  is  necesary  to  avoid  VCSEL  operation  when  there  are  two  polarization 
components  present  for  the  types  of  applications  that  we  have  proposed  in  this  program. 

A  more  genercil  problem  was  the  strong  influence  of  the  intrinsic  spontaneous  emission  noise 
on  the  spectral  features.  Many  potential  applications  require  low  noise  signals  of  good  spec¬ 
tral  purity.  However,  when  the  laser  was  operated  under  conditions  where  there  would  be  a 
strong  microwave  modulation  frequency,  the  resonance  feature  was  quite  broad.  Figure  38 
shows  a  typical  example.  The  master  laser,  and  thus  the  optical  component  at  the  locking 
frequency,  had  an  linewidth  of  under  5  MHz.  Therefore,  the  shifted  Fabry-Perot  resonance 
component  was  not  locked  but  had  a  much  broader  linewidth,  yielding  the  order  of  magni¬ 
tude  larger  linewidth  of  the  microwave  modulation  component  in  the  power  spectrum.  The 
spontaneous  emission  noise  broadens  the  shifted  resonance  feature  and  breaks  the  locking  of 
the  master  laser.  We  have  calculated  the  effects  of  a  current  modulation  at  the  resonance 
freqeuncy,  in  addition  to  the  injected  optical  signal.  The  modulation  current  has  the  effect  of 
locking  the  two  optical  components  together  and  sharply  reducing  the  linewidth.  However, 
the  modulation  component  continued  to  retain  fairly  strong  noise  sidebands.  Clearly,  the 
noise  characteristics  of  the  VCSEL  under  investigation  are  not  acceptable  for  the  types  of 
analog  systems  requiring  good  spectral  purity. 

Due  to  the  effects  of  the  polarization  instability  and  the  excitation  of  higher  order  transverse 
modes  at  higher  bias  currents,  the  VCSEL  under  investigation  had  a  relatively  narrow  range 
of  operating  conditions  where  it  had  the  potential  to  operate  as  a  single-mode  device.  Single¬ 
mode  semiconductor  laser  noise  characteristics  generally  can  be  improved  significantly  by 
pushing  the  bicis  current  far  above  threshold.  However,  the  ability  to  do  this  with  the 
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Figure  36:  Optical  spectra  showing  the  effect  of  external  optical  injection  on  the  two  polariza¬ 
tion  components.  J  =  0.31,  (a)  injection  and  measurement  parallel  to  the  strong  polarization 
component,  and  (b)  injection  and  measurement  parallel  to  the  weak  polarization  component. 
The  free-running  spectra  are  shown  offset  for  clarity. 
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Figure  37:  Optical  spectra  showing  the  effects  of  optical  injection  when  the  two  polarization 
components  are  of  similar  magnitude,  J  =  0.38.  Top,  injection  on  the  left  peak.  Middle, 
injection  on  the  right  peak.  Bottom,  free-running.  The  spectra  are  offset  for  clarity. 
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Figure  38:  Power  spectrum  at  the  beat  frequency  between  the  injection- locked  optical  signal 
component  and  the  shifted  Fabry-Perot  cavity  resonance  feature  for  the  laser  under  optical 
injection  where  the  optical  spectrum  showed  two  strong  optical  components.  The  bias  current 
was  4.8  mA  and  the  master  laser  had  a  linewidth  of  under  5  MHz. 
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VCSEL  was  limited  by  the  appearance  of  the  higher  order  transverse  modes.  This  problem 
is  not  intrinsic  to  VCSELs  and  better  devices  are  being  reported  in  the  literature. 
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5.0  VCSEL  NEAR  THRESHOLD 


The  threshold  for  laser  oscillation  is  a  prototype  for  the  disorder-to-order  transition  that 
occurs  in  a  wide  range  of  physical  systems  [46].  It  has  a  strong  analogy  with  phase  transitions 
of  systems  in  thermal  equilibrium  [47].  At  such  transition  points,  the  fluctuations  induced  by 
noise  sources  play  a  critical  role.  Noise  in  lasers  near  threshold  has  been  studied  extensively, 
but  most  analyses  have  concentrated  on  laser  systems  where  the  refractive  index  of  the 
gain  medium  is  relatively  unaffected  by  changes  in  the  optical  gain  [48].  Here,  we  present 
data  starting  from  a  free-running  condition  which  shows  a  striking  change  in  the  lineshape 
of  the  semiconductor  laser.  While  the  lineshape  is  Lorentzian  below  and  well  above  the 
threshold  region  of  laser  oscillation,  it  is  near  Gaussian  at  and  just  above  the  threshold. 
The  Gaussian  lineshape  reflects  the  fluctuations  of  the  cavity  resonance  frequency  due  to 
the  carrier  density  fluctuations.  Injection  of  an  external  optical  field  induces  spectral  holes 
and  spikes  as  well  as  an  overall  shift  in  the  Gaussian  lineshape.  The  typical  optical  spectra 
of  a  semiconductor  laser  undergoing  deterministic  chaos  is  observed  as  well  as  a  transition 
back  to  linear  dynamics  and  a  Lorentzian  lineshape  of  the  field-noise-induced  spectrum  at 
high  injection  levels.  We  observed  these  features  in  a  VCSEL.  However,  the  key  point  is  not 
that  the  laser  structure  is  very  different  from  the  conventional  edge-emitting  laser  described 
above.  It  is  that  this  laser  had  significantly  larger  values  of  b  and  7c.  The  larger  value  of  7c 
means  that  the  laser  is  noisier.  All  novel  features  observed  in  the  spectra  could  be  recovered 
using  the  coupled  equation  model. 


5.1  FREE-RUNNING  OPERATION 


Well  below  and  well  above  the  threshold  for  laser  oscillation,  the  lineshape  of  a  semiconductor 
has  been  well  studied  and  is  well  understood  [3].  At  an  injection  current  of  3.8  mA,  below  the 
oscillation  threshold,  the  laser  exhibits  a  Lorentzian  lineshape  determined  by  the  difference 
between  the  gain  and  7c,  as  shown  in  Figure  39(a).  As  hais  been  observed  previously,  the 
linewidth  of  a  semiconductor  laser  then  begins  to  broaden  as  the  injection  current  is  increased 
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Figure  39;  Measured  optical  spectra  of  a  VCSEL  as  the  bias  current  is  varied:  (a)  3.8  mA, 
just  below  the  oscillation  threshold,  (b)  4.1  mA,  above  the  oscillation  threshold,  and  (c)  4.4 
mA,  approximately  10%  above  the  oscillation  threshold.  Also  shown  are  Lorentzian  (solid 
line)  and  Gaussian  (dashed  line)  curves  fit  to  match  the  spectra  full  width  at  half  maximum. 
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[49]  due  to  the  linewidth-enhancing  effects  of  b  [3].  We  also  observe  in  the  VCSEL,  and 
this  is  a  key  experimental  feature  that  we  wish  to  emphasize  here,  that  this  broadening 
is  accompanied  by  a  change  in  the  lineshape  from  Lorentzian  to  near  Gaussian.  This  is 
shown  in  Figure  39(b)  where  the  injection  current  to  the  VCSEL  is  4.1  mA.  As  the  injection 
current  is  increased  further,  the  Lorentzian  lineshape,  accompanied  by  relaxation  resonance 
sidebands,  reappears.  This  is  the  typical  optical  spectrum  of  a  semiconductor  laser  biased 
above  threshold,  and  the  spectrum  for  an  injection  current  of  4.4  mA  is  shown  in  Figure  39(c). 
In  the  wings  the  laser  lineshape  falls  off  more  rapidly  than  the  Lorentzian  curve.  This  feature, 
along  with  the  relaxation  resonances,  has  been  previously  observed  in  semiconductor  lasers 
and  explained  by  pointing  out  that  the  linewidth  enhancement  factor,  acting  through  carrier 
density  fluctuations,  is  not  an  instantaneous  effect  [3].  Ultimately,  for  large  offset  frequencies, 
a  Lorentzian  lineshape  reflecting  only  field-noise  sources  acting  directly  on  the  optical  phase 
should  reappear.  It  has  been  previously  established  in  lasers  where  changes  in  the  gain  are 
not  accompanied  by  significant  changes  in  the  refractive  index  that  the  lineshape  always 
remains  Lorentzian  [48].  Past  work  on  the  lineshape  of  semiconductor  lasers  near  threshold 
has  shown  non-Lorentzian  lineshapes,  though  this  point  has  not  been  emphasized,  but  these 
lineshapes  were  not  so  clearly  Gaussian  as  was  observed  in  our  laser  [49]. 

We  have  calculated  the  optical  spectra  for  the  VCSEL  using  the  coupled  equation  model. 
To  emphasize  the  generic  nature  of  the  observed  phenomena,  we  use  a  simplified  form  of 
the  model.  Near  threshold,  where  the  circulating  optical  field  is  small,  it  is  more  convenient 
to  use  as  the  reference  point  the  injection  current,  Jo,  and  the  carrier  density,  Nq,  where 
Tg{No)  =  7c  when  no  optical  field  is  present.  Coupled  equations  for  the  two  quadrature 
components  of  the  optical  field  then  replace  the  equations  for  amplitude  and  phase  used  for 
the  above-threshold  case.  Near  threshold,  we  can  ignore  the  contributions  of  7„  and  7p  to 
the  damping  due  to  the  small  value  of  the  field.  We  retain  the  important  contributions  of 
the  field  dependence  to  the  cavity  resonance  frequency  shift  and  to  the  field-carrier  resonant 
coupling.  If  we  make  the  assumptions  that  7,  =  7n  J  when  J  =  1  and  that  cavity  losses  are 
due  to  output  coupling  only,  then  the  coupled  equations  take  on  a  particularly  simple  form. 

—37-  =  ^{Ei  +  hE2)n  -f  Ai  cos(flt)  -|-  Fi,  (82) 

at  L 
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(83) 


^  =  -^(6E,-£i)n- Aism(n()  +  fi, 

at  Z 

§  =  +  (84) 

Here,  Ei  and  E2  and  Fi  and  F2  are  the  quadrature  field  and  field-noise  components,  re¬ 
spectively.  The  experimentally  accessible  control  parameter  is  the  injection  current  which  is 
normalized  for  numerical  calculation  in  the  parameter  A  =  (J  —  Jo)fed/^sNo.  A  convenient 
normalization  for  the  equations  can  be  made  by  defining  Rap  =  "^cFq  and  normalizing  the 
field  equations  with  respect  to  Fq.  Values  for  the  parameters  are  derived  from  previously 
reported  data  [41]  with  6  =  6,  7c  =  5  x  10^^  s“^,  7a  =  1.1  x  10“^7c,  and 

-r — 9nFI  =  10~®7c.  (85) 

hujQ 

Using  the  simplified  equations,  we  can  recover  the  transition  from  Lorentzian  to  near  Gaus¬ 
sian  to  Lorentzian  as  A,  or  the  injection  current,  is  increased.  Calculated  optical  spectra 
of  the  free-running  laser  are  plotted  in  Figure  40.  Due  to  the  more  simplified  model,  the 
agreement  is  relatively  less  quantitative  than  for  the  edge-emitting  laser  described  above. 
Lorentzian  wings  are  calculated  for  the  near-Gaussian  lineshape  at  offset  frequencies  beyond 
25  GHz.  If,  however,  the  linewidth  enhancement  factor  is  reduced  to  6  =  3,  or  the  strength 
of  the  noise  power,  is  reduced  by  a  factor  of  2,  the  Gaussian  lineshape  at  A  =  0.1 
disappears,  as  shown  in  Figure  41.  In  either  of  these  cases  the  lineshape  remains  essen¬ 
tially  Lorentzian,  with  the  relaxation  resonance  sidebands  above  threshold,  throughout.  We 
should  point  out  that  a  change  in  the  noise  power  is  also  equivalent  to  a  simultaneous  scaling 
of  73  and  gN  because  of  the  intimate  relationship  between  Rgp  and  7c  [50]. 

Unlike  the  optical  spectrum,  the  amplitude  (intensity)  spectrum  of  the  laser  field  and  the 
spectrum  of  carrier  density  fluctuations  are  not  influenced  by  a  change  in  the  linewidth 
enhancement  factor.  Similarly,  for  the  range  of  values  considered  here,  they  show  only  minor 
changes  in  shape,  accompanied  by  an  overall  scaling  in  strength,  due  to  a  change  in  the 
amplitude  of  the  noise  source.  The  calculated  spectra  for  the  intensity  and  carrier  density 
follow  the  expected  profile  from  a  linearized  analysis  [15],  except  at  low  frequencies  due  to 
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Figure  40:  Calculated  optical  spectra  (solid  line)  and  histograms  of  the  Fabry-Perot  reso¬ 
nance  frequency  variation  based  on  fluctuations  in  the  carrier  density  (squares)  using  the 
experimentally  determined  VCSEL  parmeters  for  different  values  of  the  normalized  bias 
current  parameter:  (a)  A  =  -0.02,  (b)  A  =  0.1,  and  (c)  A  =  0.2. 
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Figure  41;  Calculated  optical  spectra  (solid  line)  and  histograms  of  the  Fabry-Perot  reso¬ 
nance  frequency  variation  based  on  fluctuations  in  the  carrier  density  (squares)  using  modi¬ 
fied  values  of  the  experimentally  determined  VCSEL  parmeters  for  the  case  of  Figure  16(b). 
In  (a)  the  linewidth  enhancement  factor,  6=3,  not  6.  and  in  (b)  the  field  noise  power  is 
reduced  bv  a  factor  of  2.  All  other  parameters  remain  the  same. 
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the  incomplete  amplitude  pinning  near  threshold.  The  amplitude  of  the  laser  field  and  the 
carrier  density  are  mutually  coupled  but  are  independent  of  the  phase  of  the  field.  While  the 
phase  fluctuations  do  not  couple  back  to  these  two  quantities,  they  are  strongly  influenced 
by  them  through  the  linewidth  enhancement  factor  and  the  dependence  of  the  gain  on  the 
carrier  density.  However,  the  near-Gaussian  lineshape  is  not  a  result  of  nonlinear  dynamics. 
All  three  equations  can  be  linearized  and  the  same  lineshape  will  appear.  In  a  semiconductor 
laser,  the  Lorentzian  lineshape  occurs  when  the  relaxation  resonances  are  well  outside  of  the 
central  linewidth  so  that  the  resonant  peaks  separated  by  the  relaxation  resonance  frequency 
overlap  weakly  [3,  51].  Near  threshold  this  approximation  will  fail.  In  cavities  characterized 
by  a  fast  photon  decay  rate,  with  correspondingly  large  field-noise  source  term,  and  a  large 
linewidth  enhancement  factor,  like  the  VCSEL  studied  here,  strong  deviations  from  the 
Lorentzian  lineshape  can  be  expected. 

To  determine  the  source  of  the  near-Gaussian  lineshape  it  is  necessary  to  analyze  the  calcu¬ 
lated  time  series  for  the  carrier  density  equation.  The  carrier  density  fluctuations  are  pro¬ 
portional  to  index  of  refraction  changes  through  the  dependence  of  the  gain  on  the  carrier 
density.  These  fluctuations  lead  to  instantaneous  variations  in  the  cavity  resonant  frequency. 
A  histogram  of  the  frequency  variations  is  also  shown  in  Figures  40(a)  and  41.  It  is  gen¬ 
erated  by  counting  the  number  of  times  the  frequency  (carrier  density)  value  falls  within  a 
frequency  bin  interval  and  plotting  the  count  as  a  function  of  frequency.  This  generates  a 
curve  with  values  proportional  to  the  fraction  of  time  that  the  particular  frequency  bin  is  the 
instantaneous  resonant  frequency.  The  curve  of  resonant  frequencies  is  dissimilar  to  the  laser 
lineshape  in  Figures  40(a)  and  (c)  and  in  Figure  41,  but  follows  the  same  Gaussian  profile  in 
Figure  40(b).  At  this  pump  level,  it  is  the  range  of  carrier  density  fluctuations  that  deter¬ 
mine  the  laser  lineshape  and  not  the  usual  phase  diffusion  fluctuations.  The  instantaneous 
cavity  resonant  frequency  fluctuations  are  caused  by  the  same  field-noise  source  terms  as  the 
phase  diffusion  fluctuations,  but  usually  they  are  not  a  significant  factor  in  determining  the 
lineshape  because  the  linewidth  of  each  of  the  relaxation  resonance  features  that  forms  the 
overall  lineshape  is  less  than  the  resonance  frequency.  Therefore,  the  normal  phcise  diffusion 
term  reasserts  its  dominance  above  threshold  even  though  the  carrier  density  fluctuations 
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persist.  Only  when  the  different  relaxation  resonance  components  show  substantial  overlap 
is  there  strong  deviation  from  a  central  Lorentzian  lineshape. 


5.2  OPERATION  UNDER  EXTERNAL  OPTICAL  INJECTION 


Having  understood  the  free-running  characteristics  of  the  VCSEL  near  threshold,  we  then 
investigated  the  effects  of  external  optical  injection.  We  first  investigated  the  spectral  changes 
when  the  VCSEL  was  biased  at  3.9  mA  where  the  free-running  spectrum  shows  only  modest 
deviation  from  a  Lorentzian  spectrum.  This  study  was  motivated,  in  part,  by  the  publication 
of  an  analysis  of  near-threshold  VCSEL  operation  which  used  an  amplifier  model,  as  opposed 
to  the  lumped-element  oscillator  model  we  use,  to  describe  some  novel  features  in  the  optical 
spectra  of  a  VCSEL  under  strong  optical  injection  [24].  The  model  predicted  that  stimulated 
emission  and  absorption  due  to  the  coherent  transfer  of  energy  significantly  enhanced  the 
semiconductor  response  and  produced  new  resonances  in  the  optical  spectrum  which  are 
distinct  from  the  relaxation  resonances  typically  observed  in  semiconductor  lasers.  Here,  we 
show  that  the  laser  cavity  plays  a  profound  role  in  the  generation  of  the  new  resonances  and 
in  the  enhancement  of  the  modulation  bandwidth.  Therefore,  the  lumped-element  model  is 
more  appropriate  than  the  travelling  wave  amplifier  model. 

When  the  free-running  laser  is  biased  near  threshold  and  is  subjected  to  a  near-resonant 
locking  field,  the  limit  \AllAl\  <  1,  and  7e  >  2U  and  7,,  then  U  ^  812  +  iclnLl2-is.  In 
this  limit.  Equations  52  and  53  can  be  solved  to  show  that  the  regen  term  dominates  the 
fwm  term.  If  the  gain  defect,  8,  is  large  compared  to  7,  and  V ,  the  central  peak  of  the  regen 
spectrum  is  shifted  from  the  free-running  spectrum.  The  difference  between  the  two  in  the 
weak  locking,  large  8  limit  has  maxima  and  minima  shifted  from  the  free-running  frequency 
by  -7,[1/6±  (1  +  1/6^)^/^],  the  same  as  predicted  by  the  amplifier  analysis  [24].  However, 
Equation  52  predicts  that  the  amplifier  analysis  is  inadequate  when  the  free-running  laser  is 
near  or  above  threshold,  and/or  when  V  becomes  comparable  to  8.  As  the  injecting  field  is 
increased  so  that  U  and  V  dominate  the  other  rates,  except  for  7c,  the  regen  term  contains 
a  strong  resonance  at  H  w  -V  =  wo  -  wl  -  }nclnLl‘2ls-  The  corresponding  resonance  at 


positive  frequencies  is  much  weaker  due  to  a  cancelling  term,  fl  —  V,  in  the  numerator.  The 
magnitude  of  the  shift  is  proportional  to  the  circulating  locked  power,  due  to  the  coherent 
field-induced  carrier  decay  rate,  7nL-  Unlike  the  amplifier  analysis,  there  is  a  multiplicative 
factor  of  approximately  for  the  resonance  shift.  This  factor  is  typically  on  the  order 

of  100  -  1000  in  semiconductor  lasers.  When  jc  ^  7si  ^  large  shift  in  the  resonance  can  be 
observed  even  if  7,  >  7nL. 

To  verify  these  predictions,  we  investigated  the  VCSEL  biased  near  threshold.  Many  of  the 
characteristics  of  this  VCSEL  have  been  described  previously  [39,  40].  Figure  42  shows  the 
output  spectrum  and  the  regen  spectrum  when  the  laser  is  biased  at  4.8  mA,  well  above 
the  threshold  value  of  3.9  mA.  Total  output  power  at  this  injection  current  is  w  0.35  mW 
and  the  coherent  output  power  is  ~  0.3  mW.  Here,  we  can  determine  the  key  dynamic 
parameters  of  the  laser  from  its  spectra.  Figure  42  shows  the  good  agreement  between  data 
and  model  that  is  achieved  using  the  determined  parameters,  6  =  6,  7c  =  5.5  x  10^^  s“\ 
73  =  5  X  10^  s"\  7„  =  3.5P  X 10^  s"^  and  7p  =  5.3P  x  10®  s"U  P  in  the  formulas  refers  to  the 
coherent  output  power  in  milliwatts  and  the  uncertainty  for  the  parameters  is  ±20%.  The 
enhancement  factor  for  the  frequency  shifts,  relative  to  the  amplifier  case,  is  approximately 
330. 

For  the  injection  locking  measurements,  the  bias  current  was  set  to  3.9  mA  and  the  output 
from  a  tunable,  narrow  linewidth,  external  cavity  laser  (New  Focus  Model  6126)  was  injected 
into  the  VCSEL.  Output  power  from  the  free-running  VCSEL  was  w  0.02  mW.  Near  thresh¬ 
old  operation  insured  that  |Ao/A^l  C  1.  Here,  we  concentrate  on  optical  injection  at  the 
free-running  frequency.  Consistent  results  were  obtained  for  detuned  injection.  Figure  43 
shows  free-running  and  injection-locked  optical  spectra  taken  with  a  high  finesse  optical 
spectrum  analyzer  (Newport  Model  SR-240C)  with  w  60  MHz  resolution.  The  resonance 
feature  shifts  to  lower  frequencies  as  the  injection  power  increases. 

Features  from  these  spectra  can  be  compared  with  the  predictions  of  the  model.  The 
injection-locked  output  power  can  be  determined  from  the  spectra  and  the  measured  output 
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Figure  42:  (a)  The  optical  spectrum  of  the  free-running  VCSEL  biased  well  above  threshold 
and  (b)  the  regeneratively  amplified  spectrum  due  to  the  injection  of  a  weak  optical  probe. 
The  experimental  data  is  compared  with  the  calculated  spectra  using  the  parameters  listed 
in  the  text. 
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Figure  43:  Representative  optical  spectra  of  the  free-running  VCSEL  biased  near  threshold 
(a),  and  under  increasing  injected  power  from  a  narrowband  la^er  tuned  to  the  free  running 
oscillation  frequency  (b)-(e).  The  injected  power  ratio  of  (b):(c):(d):(e)  is  1:10:100:1000. 
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power.  Figure  44  shows  the  dependence  of  the  injection-locked  coherent  power  and  the  fre¬ 
quency  shift  of  the  resonance  on  the  power  injected  into  the  VCSEL.  Because  we  are  unable 
to  independently  determine  the  coupling  parameter,  Tj,  we  can  only  make  a  relative  measure¬ 
ment  of  the  injection  power  and  the  relative  uncertainty  is  ±50%.  Using  the  experimentally 
determined  parameters,  model  calculations  for  both  the  injection-locked  coherent  power  and 
the  frequency  shift  are  in  good  agreement  with  the  data.  At  high  injection  power  levels,  the 
offset  frequency  scales  linearly  with  the  injection- locked  coherent  power,  and  both  scale  with 
the  cube  root  of  the  injected  power.  Even  at  the  highest  injected  power  measured,  where  the 
injection-locked  coherent  power  is  ~  0.07  mW,  the  field-induced  enhancements  of  the  decay 
rate  are  7„l  «  0.057,  and  jpi  «  0.0757,.  These  values  coincide  with  the  expected  values  at 
the  lower  operating  power. 

The  lumped-circuit  analysis  of  laser  oscillation  assumes  that  spatial  effects  can  be  averaged 
over  the  mode  profiles.  It  hcis  given  excellent  quantitative  agreement  with  a  wide  variety 
of  single-mode  semiconductor  lasers,  including  Fabry-Perot  edge  emitting  lasers  with  large 
output  coupling  [7,  13].  The  linearized  treatment,  like  that  given  here,  fails  to  accurately 
reproduce  the  central  linewidth  of  the  free-running  laser  and,  more  generally,  the  dynamics 
whenever  a  resonance  becomes  unstable,  but  it  shows  good  accuracy  in  the  prediction  of 
the  positions  and  shifts  of  the  resonance  frequencies.  The  laser  cavity  strongly  enhances  the 
frequency  shifts  induced  by  the  injection  field  beyond  what  is  expected  from  the  increased 
stimulated  emission  due  to  the  stronger  oscillating  field.  Phase  and  amplitude  characteristics 
must  be  analyzed  for  a  detailed  quantitative  understanding  of  the  spectral  features. 

When  the  free-running  Gaussian-lineshape  laser  is  subjected  to  external  optical  injection, 
several  spectral  changes  occur.  As  before,  we  first  concentrate  on  resonant  injection.  The 
first  observable  change  as  the  injection  level  is  increased  is  a  shifting  of  the  spectrum  to 
lower  optical  frequencies.  This  reflects  the  drop  in  the  overall  carrier  density  due  to  the 
higher  circulating  field  and  enhanced  stimulated  emission.  As  the  injection  level  is  increased 
further,  the  spectrum  shows  a  new  spike  at  negative  offset  frequencies  and  a  dip  at  positive 
offset  frequencies  in  addition  to  the  overall  shift  and  a  spike  at  the  injection  frequency,  as 


112 


Injected  Power  (arb.  units) 


Figure  44:  The  dependence  of  the  injection-locked  power,  (Dj-data  and  ( —  )-modeI,  and 

the  shift  of  the  resonance  from  the  free-running  frequency,  (A)-data  and  ( - )-model,  for 

the  VCSEL  biased  near  threshold  when  subject  to  an  injected  signal  at  the  free-running 
oscillation  frequency. 
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shown  in  Figure  45(a).  The  free- running  spectrum  is  also  shown  for  reference.  The  spike  and 
dip  are  oppositely  offset  by  1  GHz,  approximately  the  relaxation  resonance  frequency.  As 
the  excitation  level  is  increased  further,  spikes  and  dips  appear  at  multiples  of  the  original 
offset,  corresponding  to  the  different  resonance  features.  At  higher  excitation  levels,  the 
spectra  are  similar  to  the  chaotic  dynamics  spectra  discussed  above.  Figure  45(b).  Finally, 
at  very  high  excitation  levels,  the  spectral  structure  has  shifted  away  from  the  excitation 
peak  and  assumed  a  more  Lorentzian  lineshape.  Figure  39(c).  Calculated  spectra  show  the 
progression  from  shift  with  spike  and  dip.  Figure  46(a),  to  chaotic-like.  Figure  46(b),  to 
shifted  Lorentzian,  Figure  46(c). 

We  have  previously  presented  the  shift  of  the  Fabry-Perot  resonance  feature  due  to  optical 
injection  when  this  laser  was  initially  biased  at  3.9  mA,  at  the  laser  threshold  [41].  There,  the 
initial  deviation  from  Lorentzian  was  not  large  but  the  spectra  clearly  showed  a  narrowing 
and  a  progression  to  a  more  Lorentzian  lineshape  as  the  injected  power  was  increased.  This 
effect  can  also  be  seen  as  the  laser  frequency  is  tuned  across  the  resonance  at  strong  injection 
levels  at  the  higher  operating  current  described  by  the  spectra  presented  here.  Figure  47 
plots  this  trend  under  strong  external  injection.  For  large  positive  detunings  the  Gaussian 
lineshape  is  retained  but  there  is  a  shift  of  the  feature.  As  the  detuning  is  decreased  the 
feature  narrows,  with  increased  amplitude.  For  injection  close  to  the  edge  of  the  locking 
range,  near  —24  GHz,  the  emission  feature  becomes  weak,  indicating  stable,  injection-locked 
operation.  The  smooth  progression  of  the  shift  of  the  Fabry-Perot  resonance  feature  is 
summarized  in  Figure  48.  For  positive  detunings,  the  shift  decreases  with  increasing  offset 
as  the  injection  undergoes  less  amplification  in  the  cavity.  For  negative  detunings  the  shift 
is  greater  than  but  approaching  the  magnitude  of  the  injection  offset  as  the  offset  decreases 
to  more  negative  values.  At  the  edge  of  the  locking  range  this  smooth  progression  abruptly 
halts  and  the  laser  displays  a  weaker,  more  complex  spectrum.  The  shift  of  the  cavity 
resonance  frequency  and  the  associated  damping  are  caused  by  the  reduced  carrier  density 
under  strong  external  injection.  The  dynamics  are  strongly  influenced  by  the  magnitude  of 
the  shift  and  damping. 
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Figure  45:  Measured  optical  spectra  of  the  VCSEL  biased  at  4.1  mA  and  subjected  to 
external  optical  injection.  The  ratio  of  injected  power  is:  (a)  1,  (b)  10,  (c)  940.  The  free- 
running  spectrum  is  also  shown  for  reference. 
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Figure  46:  Calculated  optical  spectra  of  the  VCSEL  biased  at  4.1  mA  and  subjected  to 
external  optical  injection.  The  ratio  of  injected  optical  power  is:  (a)  1,  (b)  4,  (c)  100. 
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Figure  47;  Spectra  of  the  VCSEL  under  strong  optical  injection  as  the  detuning  is  varied: 
12  GHz  (squares),  —12  GHz  (diamonds),  and  —24  GHz  (triangles).  Each  spectrum  has 
two  principal  features,  the  regeneratively  amplified  injected  signal  and  the  regeneratively 
amplified  spontaneous  emission  at  the  shifted  Fabry-Perot  resonance  peak  of  —8.5  GHz, 
—19.5  GHz,  and  —29.5  GHz,  respectively.  The  —12  GHz  offset  excitation  spectrum  also 
shows  a  weak  four-wave  mixing  peak  at  —4.5  GHz. 
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Figure  48:  Measured  pushing  of  the  Fabry-Perot  resonance  as  a  function  of  the  detuning 
of  the  injected  signal.  The  diamonds  are  for  an  injected  power  approximately  an  order  of 
magnitude  stronger  than  the  squares.  The  solid  line  is  a  slope  of  one  line  at  negative  offsets 
and  slope  of  zero  line  at  positive  offsets  to  aid  the  eye. 
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6.0  CONCLUSIONS 


The  work  done  under  this  program  has  identified  the  key  output  characteristics  and  nonlinear 
dynamics  of  a  semiconductor  laser  subject  to  external  optical  injection.  It  has  confirmed  that 
the  coupled  equation  model  for  a  semiconductor  laser  under  optical  injection  predicts  and/or 
reproduces  the  spectral  characteristics  that  we  observe.  This  has  allowed  us  to  identify  the 
underlying  deterministic  dynamics  that,  to  a  large  extent,  determine  the  output  character¬ 
istics  of  the  semiconductor  laser  under  optical  injection.  We  have  verified  that  bandwidth 
enhancement  and  noise  reduction  can  be  simultaneously  achieved  in  an  injection- locked  laser 
diode.  Certainly,  there  is  no  sacrifice  in  the  dynamic  range  of  an  appropriately  injection- 
locked  semiconductor  laser  and  there  can  be  considerable  enhancement  as  well  as  reduction 
of  the  frequency  chirp  which  limits  the  effectiveness  of  laser  diodes  in  certain  applications.  In 
the  region  of  unstable  dynamics,  we  have  been  able  to  generate  a  broadly  tunable  modulation 
signal  when  the  master  laser  is  appropriately  tuned  to  the  region  of  limit  cycle,  period  one, 
oscillations  under  strong  injection.  However,  the  microwave  signal  generated  by  the  beating 
of  optical  components  is  quite  noisy  when  the  leiser  is  excited  to  display  unstable  dynamics 
due  to  the  strong  effects  of  spontaneous  emission  noise.  More  work  is  necessary  to  improve 
the  operating  characteristics  in  this  region.  We  have  observed  that  feedback  can  modify  the 
modulation  characteristics.  Perhaps  it  can  be  used  to  improve  the  stability  of  the  unstable 
modulation  feature  as  well.  Several  novel  spectral  characteristics  have  been  observed  in  a 
VCSEL.  These  include  a  Gaussian  lineshape  near  threshold,  strong  frequency  pushing  of 
the  Fabry-Perot  resonance  signal,  and  a  locking/unlocking  bistability.  All  observed  features 
have  been  understood  by  comparing  measured  spectra  with  the  predictions  of  the  model. 
Under  strong  injection  locking,  semiconductor  lasers  are  calculated  to  be  able  to  achieve 
bandwidths  broadly  tunable  in  the  millimeter-wave  region,  with  tunable  modulation  char¬ 
acteristics  and  reduced  noise  relative  to  the  free-running  case.  This  technique  may  be  very 
useful  for  augmenting  the  performance  of  available  laser  diodes  so  that  they  can  be  useful 
in  a  variety  of  applications. 
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